Let $n$ be a natural number and $a$ be a real number, where $a \neq 0$. Suppose that the integral expression $x ^ { n } + y ^ { n } + z ^ { n } + a ( x y + y z + z x )$ can be expressed as the product of $x + y + z$ and an integral expression $P$ in $x$, $y$ and $z$, i.e.
$$x ^ { n } + y ^ { n } + z ^ { n } + a ( x y + y z + z x ) = ( x + y + z ) P .$$
We are to find the values of $n$ and $a$.
The equality (1) holds for all $x$, $y$ and $z$. So, consider for example, two triples of $x$, $y$ and $z$ that satisfy $x + y + z = 0$:
$$x = y = 1 , \quad z = - \mathbf { A }$$
and
$$x = y = - \frac { \mathbf { B } } { \mathbf { C } } , \quad z = 1 .$$
By substituting each triple in (1), we obtain the following two equations:
$$\left( - \mathbf { A } \right) ^ { n } = \mathbf { D } \text{ a} - \frac{\mathbf{E}}{\mathbf{F}}$$
$$\left( - \frac { \mathbf { B } } { \mathbf { C } } \right) ^ { n } = \frac { \mathbf{E} }{ \mathbf{F} } \text{ a} - \frac { \mathbf { H } } { \mathbf { I } } .$$
From these two equations, we get an equation in $a$. Solving this, we obtain $a = \mathbf { K }$ and hence by the first equation that $n = \mathbf { L }$.
Conversely, when $a = \mathbf{K}$ and $n = \mathbf{L}$, there exists a $P$ such that (1) holds, and hence these values of $a$ and $n$ are the solution.