Putting m=a2, n=b2, p=c2 then m,n,p>0. We have the system
{mx+ny+pz=1xy+yz+zx=1.
From this system, we get mx=1−ny−pz and m(xy+yz+zx)=m, therefore
ny2+y(nz+pz−mz−1)+m−z+pz2=0(∗)
Equation (∗) has only one solution when (nz+pz−mz−1)2−4n(m−z+pz2)=0.
This equation is equivalent to ((m−n−p)2−4np)z2+2(m+n−p)z+1−4mn=0. (∗∗)
Equation (∗∗) has only one solution when
(m+n−p)2=(1−4mn)((m−n−p)2−4np).
or by simple algebra m2+n2+p2+1=2(mn+np+pm).
Conversely, with this condition, putting back into (∗) and (∗∗), we see that the system has only one solution: (x,y,z)=(n+p−m,p+m−n,m+n−p).
Thus, m2+n2+p2+1=2(mn+np+pm)⇔a4+b4+c4+1=2(a2b2+b2c2+c2a2) is the necessary and sufficient condition for the system to have only one solution.
Moreover, from this condition, we have
2(a2b2+b2c2+c2a2)>a4+b4+c4⇔(a+b+c)(a+b−c)(b+c−a)(c+a−b)>0.
From which follows that all of the numbers a+b−c, b+c−a, c+a−b are positive, so a,b,c are the sides of certain triangle. (Q.E.D)