Estas son todas las soluciones con $$ -2 \leq x < y < z < 1700. $$
Por alguna razón, hay sólo un par de soluciones con algunos positivos entradas y algunos negativos, que probablemente tiene una breve prueba. De todos modos, una vez que los tres son positivos y distintos, es fácil ver que podemos exigir $x<y<z.$
x^2 + y^2 + z^2 x y z x^2 + y^2 + z^2 FACTORED
6 -2 -1 1 6 = 2 3 gcd(x,y,z) 1
2 -1 0 1 2 = 2 gcd(x,y,z) 1
6 -1 1 2 6 = 2 3 gcd(x,y,z) 1
686 7 14 21 686 = 2 7^3 gcd(x,y,z) 7
20250 55 85 100 20250 = 2 3^4 5^3 gcd(x,y,z) 5
20250 60 75 105 20250 = 2 3^4 5^3 gcd(x,y,z) 15
31250 75 100 125 31250 = 2 5^6 gcd(x,y,z) 25
73002 115 161 184 73002 = 2 3 23^3 gcd(x,y,z) 23
140250 175 205 260 140250 = 2 3 5^3 11 17 gcd(x,y,z) 5
156750 155 250 265 156750 = 2 3 5^3 11 19 gcd(x,y,z) 5
332750 275 330 385 332750 = 2 5^3 11^3 gcd(x,y,z) 55
384846 259 266 497 384846 = 2 3 7^3 11 17 gcd(x,y,z) 7
1647750 650 715 845 1647750 = 2 3 5^3 13^3 gcd(x,y,z) 65
1825346 679 776 873 1825346 = 2 97^3 gcd(x,y,z) 97
2409066 736 943 989 2409066 = 2 3^2 11 23^3 gcd(x,y,z) 23
3188646 891 1053 1134 3188646 = 2 3^13 gcd(x,y,z) 81
3188646 918 999 1161 3188646 = 2 3^13 gcd(x,y,z) 27
2676750 775 790 1205 2676750 = 2 3 5^3 43 83 gcd(x,y,z) 5
3668250 860 1195 1225 3668250 = 2 3 5^3 67 73 gcd(x,y,z) 5
4276866 1063 1184 1321 4276866 = 2 3 11^2 43 137 gcd(x,y,z) 1
5343750 1150 1175 1625 5343750 = 2 3^2 5^6 19 gcd(x,y,z) 25
6885902 1359 1510 1661 6885902 = 2 151^3 gcd(x,y,z) 151
x^2 + y^2 + z^2 x y z x^2 + y^2 + z^2 FACTORED