After the change of variables (a,b,c)=(x+y,y+z,z+x) it suffices to prove that
a2+b2+c2+1+∣a−b∣2ab+1+∣b−c∣2bc+1+∣c−a∣2ac≥0
for all real numbers a, b, c. Let (A,B,C)=(∣b−c∣,∣c−a∣,∣a−b∣), since (ab)⋅(bc)⋅(ca)≥0 we may assume that ab≥0. Since A,B,C≥0 and A+B=∣b−c∣+∣c−a∣≥∣a−b∣=C, we have
1+Cab≥(1+A)(1+B)ab
(a2+b2+1+C2ab−(1+Ba+1+Ab)2)+(c+1+Ba+1+Ab)2≥0
Second Solution. Since we can change (a,b,c) by (−a,−b,−c) we can assume that a≥b≥0≥c. Then, if we put C=∣a−b∣, A=∣b−c∣, we can find ∣a−c∣=A+C. Our inequality would indeed be
a2+b2+c2+1+C2ab+1+A2bc+1+C+A2ac≥0
It is a quadratic expression with respect to c and its discriminant is
(1+A2b+1+C+A2a)2−4(a2+b2+1+C2ab)
Then, (1+A2b)2≤4b2, (1+C+A2a)2≤4a2, and 2⋅1+A2b⋅1+C+A2a≤1+C8ab. Hence, the discriminant is non-positive and we are done. ■