I'll give you a solution which uses the theory of quadratic forms on the field of $p$-adic numbers $\mathbb Q_p$. If we are able to prove that there are no solutions of the equation $p=x^2-xy+y^2$ in $\mathbb Q_p$ then, clearly, there are no integers $x,y$ satisfying the equation, since $\mathbb Z\subset \mathbb Q_p$. The first thing we have to do is to find an equivalent diagonal expression of the quadratic form $f(x,y)=x^2-xy+y^2$, in order to easily compute the local invariants of $f$. With respect to the basis $\mathcal B=\{e_1=(1,0),e_2=(1,2)\}$ we can see that the expression of $f$ is diagonal. More precisely $f$ is equivalent to $g(x,y)=x^2+3y^2$. The following theorem holds:
Theorem: Given $a\in\frac{\mathbb Q_p^\times}{\mathbb Q_p^{\times,2}}$ and $f$ a non-degenerate quadratic form of rank $n$, discriminant $d$ and Hasse-invariant $\epsilon$, the equation $f(x_1,\ldots,x_n)=a$ has a solution over $\mathbb Q_p$ if and only if:
- $n=1$ and $a=d$
- $n=2$ and $(a,-d)=\epsilon$
- $n=3$ and $a\neq -d$ or $a=-d$ and $(-1,-d)=\epsilon$
- $n\geq 4$ always
The discriminant and the Hasse-invariant are the same for equivalent forms, so we have that: $d_f=3$ and $\epsilon_f=(1,3)$. Since $p\in\frac{\mathbb Q_p^\times}{\mathbb Q_p^{\times,2}}$ and $n=2$, if we prove that $(p,-d_f)\neq \epsilon_f$ we have done. It's not too much difficult to prove that $(p,-3)=\left(\frac{-3}{p}\right)$, and using the hint of your book: $(p,-3)=-1$ for an odd prime $p$ such that $p\equiv 2$ mod $3$. While we have that $\epsilon_f=(1,3)=1$, since $1$ is a square.
So there are no solutions $x,y\in\mathbb Q_p$ to the equation $f(x,y)=p$, and therefore there are no solutions over integers.