Nota. El teorema de Wilson afirma que $(p-1)!\equiv -1\ (\text{mod }p)$ para cualquier primo $p$.
Si $x_k^2\equiv -7\pmod{2^k}$, entonces $x_k^2\equiv -7\pmod{2^{k+1}}$ o bien $x_k^2\equiv -7+2^k\pmod{2^{k+1}}$. En el primer caso podemos tomar directamente $x_{k+1}=x_k$ y, en el segundo, definimos $x_{k+1}=x_k+2^{k-1}$ y tenemos que \[x_{k+1}^2=x_{k}^2+2^kx_k+2^{2k-2}\equiv x_k^2+2^k\equiv (-7+2^k)-2^k\equiv -7\pmod{2^{k+1}},\] donde hemos usado además que $x_k$ debe ser impar. No obstante, esta forma de pasar de $x_k$ a $x_{k+1}$ es cierta únicamente para $k\geq 3$ ya que debe cumplirse que $2k-2\geq k+1$ para que $2^{2k-2}\equiv 0\pmod{2^{k+1}}$. La prueba por inducción se sigue entonces de los casos base $1\leq k\leq 3$, en los que es basta con tomar $x_1=x_2=x_3=1$.
Vamos a ver que esto implica que $a=b=c=0$. En efecto, si alguno de los tres números $a,b,c\in\mathbb{Z}$ es no nulo, entonces podríamos haber comenzado suponiendo que la solución $(a,b,c)$ es tal que la suma $a^2+b^2+c^2\gt 0$ es lo más pequeña posible (de entre todas las soluciones no nulas habrá una que cumpla esto), pero entonces $(x,y,z)$ es una solución con \[x^2+y^2+z^2=\frac{1}{9}(a^2+b^2+c^2)\lt a^2+b^2+c^2,\] en contradicción con el hecho de que $a^2+b^2+c^2$ es mínimo. Esto es lo que se llama técnica del descenso infinito o principio de minimalidad.