Suppose $f(x) = ax^3 + bx^2 + cx + d$ is a cubic equation with roots $\alpha, \beta, \gamma.$ Then we have:
$\alpha + \beta + \gamma= -\frac{b}{a}\quad (1)$
$\alpha\beta + \beta\gamma + \gamma\alpha = \frac{c}{a}\quad (2)$
$\alpha\beta\gamma = -\frac{d}{a}\quad (3)$
We can find $\alpha^2\beta + \beta^2\gamma + \gamma^2\alpha + \alpha^2\gamma + \gamma^2\beta + \beta^2\alpha$ in terms of $a,b,c,d$ with the formula:
$$ \alpha^2\beta + \beta^2\gamma + \gamma^2\alpha + \alpha^2\gamma + \gamma^2\beta + \beta^2\alpha = (\alpha+\beta+\gamma)(\alpha\beta+\alpha\gamma+\beta\gamma) - 3\alpha\beta\gamma $$ $$=\left(\frac{-b}{a}\right) \left(\frac{c}{a}\right) - 3\left(-\frac{d}{a}\right).$$
But I was wondering if there was some way to find $ \alpha^2\beta + \beta^2\gamma + \gamma^2\alpha\ $ and therefore also $\ \alpha^2\gamma + \gamma^2\beta + \beta^2\alpha\ $ in terms of $a,b,c,d,\ $ with some algebraic manipulation, i.e. without finding the roots with a cubic formula?
Notice that there are two possible values of $\alpha^2\beta+\beta^2\gamma+\gamma^2\alpha,$ namely $\alpha^2\beta+\beta^2\gamma+\gamma^2\alpha = \beta^2\gamma+\gamma^2\alpha+\alpha^2\beta = \gamma^2\alpha+\alpha^2\beta+\beta^2\gamma$ and $\alpha^2\gamma+\gamma^2\beta+\beta^2\alpha = \gamma^2\beta+\beta^2\alpha+\alpha^2\gamma = \beta^2\alpha + \alpha^2\gamma+\gamma^2\beta.$
You know the elementary symmetric polynomials evaluated at the roots (by the Vieta relations): $$ \begin{aligned}s_1&=\alpha+\beta+\gamma=-b/a,\\ s_2&=\alpha\beta+\beta\gamma+\gamma\alpha=c/a,\\ s_3&=\alpha\beta\gamma=-d/a. \end{aligned} $$ The fundamental theorem of symmetric polynomials says that every symmetric polynomial can be written in terms of the elementary ones.
You are interested in finding $$ \begin{aligned} u&=\alpha^2\beta+\beta^2\gamma+\gamma^2\alpha,&\text{or}\\ v&=\alpha^2\gamma+\gamma^2\beta+\beta^2\alpha. \end{aligned} $$ The problem, as explained by the others, is that you cannot tell which is which because $u$ and $v$ only follow cyclic symmetry.
However, the combinations $u+v$ and $uv$ are fully symmetric. A banal (but a bit tedious) calculation shows that $$ \begin{aligned} u+v&=s_1s_2-3s_3,\\ uv&=s_1^3s_3-6s_1s_2s_3+s_2^3+s_3^2. \end{aligned} $$ This means that we know the coefficients of the quadratic polynomial $$ p(x)=(x-u)(x-v)=x^2-[u+v]x+uv $$ that has $u$ and $v$ as its roots. All you need to do is plug in the known values of $s_1,s_2,s_3$ into the formulas above, and solve the quadratic $p(x)=0$. The roots are the two choices.