The inequality is equivalent to
abc(a+1)(b+1)(c+1)c(a+1)+a(b+1)+b(c+1)≥(1+abc)23,
or after simplifications to
(1+abc)2(ab+bc+ca+a+b+c)≥3abc(ab+bc+ca+a+b+c+abc+1).
We put m=a+b+c, n=ab+bc+ca and x3=abc, to obtain:
(m+n)(1+x3)2≥3x3(x3+m+n+1),
or (m+n)(x6−x3+1)≥3x3(x3+1).
From the inequality of arithmetic–geometric mean we have m≥3x and n≥3x2, and so m+n≥3x(x+1). Therefore, it is enough to prove that
⇔⇔⇔3x(x+1)(x6−x3+1)≥3x3(x+1)(x2−x+1)x6−x3+1≥x2(x2−x+1)x6−x4−x2+1≥0⇔x4(x2−1)−(x2−1)≥0(x2−1)(x4−1)≥0⇔(x2+1)(x2−1)2≥0,
which is valid.