Para un pequeño proyecto en el que estoy trabajando, tengo el deseo de encontrar las soluciones para $$x^2+y^2=z(4z+1)$$ en números naturales $x,y,z$. Quiero automatizar la búsqueda de soluciones para $z$ hasta un valor máximo igual de eficiente como sea posible, corriendo a través de las $z$ valores de 1 y, a continuación, tratar de encontrar posibles $x$$y$.
Una cosa que he encontrado es que me puede hacer caso omiso de $z \equiv 3,6,7 \mod 8$, como la suma de dos cuadrados, sólo puede ser $\equiv 0,1,2,4,5 \mod 8$. Pero me pregunto qué otros criterios que puede utilizar para excluir $z$ valores. También me pregunto si para determinado$z$, $x$ I puede ignorar, así que yo sólo la comprobación de las $x$ valores que puedan dar una solución.