El siguiente método, que utiliza el teorema de Cayley-Hamilton, puede considerarse como el uso secreto de las identidades de Newton, pero en mi opinión es más transparente:
Sea $$ X = \left(\begin{matrix} x & 0 & 0 \\ 0 & y & 0 \\ 0 & 0 & z \\ \end{matrix} \right) \,.$$
Nos interesan los rastros de $X^n$ . El polinomio característico de $X$ es $$p(\lambda) = \det(\lambda - X) = \lambda^3 - (x+y+z) \lambda^2 + (xy + yz + zx)\lambda - xyz\,.$$ El segundo coeficiente es $-12$ el tercer coeficiente se obtiene calculando $(x+y+z)^2 - (x^2 + y^2 + z^2)$ , dando $66$
Para encontrar $xyz$ se puede hacer álgebra como en otras respuestas, o utilizar el teorema de Cayley-Hamilton, que dice que $p(X) = 0$ . Tomando la traza, se obtiene $3xyz = 12 - 12\cdot 12+12\cdot 66$ o $xyz = 220$ .
Así (de Cayley-Hamilton) $$ p(X) = X^3 - 12X^2 + 66X - 220 = 0\,.$$ A partir de esta identidad, es fácil obtener la relación de recursión $$S_n = x^n + y^n + z^n = \mathrm{tr} [X^n] = \mathrm{tr}[X^{n-3} (12 X^2 - 66X + 220) ] = 12 S_{n-1} - 66 S_{n-2} + 220 S_{n-3}$$ dando $$S_4 = (12-66+220)\cdot 12 = 1992$$