Without loss of generality, we may assume gcd(a,b,c) = 1.

(otherwise, if d=gcd(a,b,c) the for a'=a/d, b'=b/d, c'=c/d, the equation still holds for a', b', c' and a'b'c' is still a cube if only if abc is a cube).

We multiply equation by abc, we have:

  \(a^2c+b^2a+c^2b=3abc\)(*)

if \(abc=\pm1\), the problem is solved.

Otherwise, let p be a prime divisor of abc. Since gcd(a,b,c)=1, the (*) implies that p divides exactly two of a, b,c. By symetry, we may assume p divides a, b but not c. Suppose that the lagest powers of p dividing a, b are m, n, respecively.

If n < 2m, then \(n+1\le2m\)  and \(p^{n+1}\)| \(a^2c,b^2c,3abc\). Hence \(p^{n+1}\)|\(c^2b\), forcing \(\)\(p\)|\(c\) (a contradiction). If n > 2m, then \(n\ge2m+1\) and \(p^{2m+1}\)|\(c^2b,b^2a,3abc\). Hence \(p^{2m+1}\)|\(a^2c\), forcing \(p\)|\(c\) (a contradicton). Therefore n = 2m and \(abc=\Pi p^{3m}\), \(p\)|\(abc\), is a cube.