La pregunta sería, naturalmente, ser muy sencillo si hubo un $2xy$ en lugar de un $4xy$. Entonces es simplemente una cuestión de hacer:
$$ x^2+2xy+y^2=1\\ (x+y)^2=1\\ \sqrt{(x+y)^2}=\sqrt{1}\\ |x+y|=1 $$
Claramente, existen infinitos pares de números enteros que se diferencian por una, así que hay un número infinito de entero soluciones a $x^2+2xy+y^2=1$. Por desgracia, el mismo principio no se aplica a $x^2+4xy+y^2=1$ donde si puedo intento de construir una similar a prueba todo lo que puedo hacer es:
$$ x^2+4xy+y^2=1\\ (x+y)^2+2xy=1\\ 2xy=1-(x+y)^2\\ 2xy=(1+x+y)(1-x-y) $$
Y, a partir de ahí, no tengo idea de cómo proceder para completar la prueba de que hay un número infinito de entero de soluciones. Me pregunto si me estoy acercando a la pregunta enteramente en el camino equivocado, o si simplemente estoy perdiendo algo. Un empujón en la dirección correcta sería muy apreciada!