Sea $d$ el máximo común divisor de $x$ e $y$ y escribamos $x=a\cdot d$ e $y=b\cdot d$, donde $a$ y $b$ son enteros primos entre sí. Si dividimos ambos miembros de la ecuación entre $d^2$, obtenemos \[a^2-b^2=2abz.\] Ahora vamos a probar que $a=\pm 1$ y $b=\pm1$. Si $a$ fuera distinto de $\pm 1$, entonces existiría un primo $p$ que divide a $a$, luego $p$ también dividiría a $b^2=a^2-2abz=a(a-2bz)$, contradiciendo que $a$ y $b$ son primos entre sí, luego hemos probado por reducción al absurdo que $a=\pm 1$ y, de la misma forma $b=\pm 1$. De aquí que $x=\pm d$ e $y=\pm d$, es decir, $x$ e $y$ son iguales u opuestos. Tenemos entonces las siguientes soluciones:
Como la ecuación es simétrica en $x$ e $y$, deducimos que todas las soluciones son $(p-p^2,p-1)$, $(p-1,p-p^2)$, $(0,0)$, $(p+1,p^2+p)$, $(p^2+p,p+1)$ y $(2p,2p)$.
En definitiva, hemos probado que las únicas ternas de números que son solución del problema son las de la forma $(0,1,p)$, $(1,0,p)$ para $p$ cualquier número primo y $(5,2,5)$, $(2,5,5)$, $(5,26,13)$ y $(26,5,13)$.