Si $\alpha, \beta, \gamma$ son las raíces de $p(x)=ax^3+bx^2+cx+d=0$ entonces $$p(\alpha)=p(\beta)=p(\gamma)=p(\alpha)+p(\beta)+p(\gamma)=0$$
Ahora $S_n=\alpha^n+\beta^n+\gamma^n$ y la última ecuación nos dice que $$aS_3+bS_2+cS_1+3d=0$$
La suma de cuadrados $S_2$ es fácil de encontrar a partir de los polinomios simétricos, y $S_1$ es conocido. Esta observación evita tener que recordar complicadas fórmulas para las sumas de potencias. (Podríamos poner $S_0=3$ )
En efecto, para las potencias superiores podemos observar
$$\alpha^r p(\alpha)+\beta^r p(\beta)+\gamma^rp(\gamma)=0$$ que da $$aS_{r+3}+bS_{r+2}+cS_{r+1}+dS_r=0$$ que nos permite calcular sucesivamente las sumas de potencias.
Además, podemos calcular la fórmula para $S_2$ utilizando una cuadrática, como si sólo se tratara de dos raíces, de modo que $$aS_2+bS_1+2c=0$$ con $S_1=-\frac ba$ [1], de modo que $$S_2=\frac {b^2}{a^2}-2\left(\frac ca\right)$$ y esto sigue siendo cierto para cualquier número mayor de raíces. Del mismo modo, la fórmula para tres raíces que podemos derivar de la cúbica se aplica también para cuatro o más raíces.
[1] incluso podríamos haberlo deducido de la ecuación lineal $aS_1+b=0$