The first term y(2z+x)x(2x−y)+1=y(2z+x)2(x2+yz). The similar transformation of two other terms yields:
f(x,y,z)=y(2z+x)x2+yz+z(2x+y)y2+zx+x(2y+z)z2+xy≥2
By Cauchy-Schwarz inequality for positive x1,…,xn
(x1+⋯+xn)(x1a12+⋯+xnan2)≥(a1+⋯+an)2(1)
Therefore,
g(x,y,z)=y(2z+x)x2+z(2x+y)y2+x(2y+z)z2≥3(xy+yz+zx)(x+y+z)2
h(x,y,z)=2z+xz+2x+yx+2y+zy≥2(x2+y2+z2)+xy+yz+zx(x+y+z)2
Therefore,
f=g+h≥(x+y+z)2(3(xy+yz+zx)1+2(x2+y2+z2)+xy+yz+zx1)=3(xy+yz+zx)(2(x2+y2+z2)+xy+yz+zx)2(x+y+z)4=3(xy+yz+zx)(2(x+y+z)2−3(xy+yz+zx))2(x+y+z)4≥(23(xy+yz+zx)+2(x+y+z)2−3(xy+yz+zx))22(x+y+z)4=2
(in the last inequality we applied AM-GM inequality). Done.