Réponse rédigée (en espérant qu'elle soit juste) :
Notons
)
l'exposant de la plus grande puissance de 2 divisant x.
Soient :
)
d'où

)
d'où

a et b impairs
Supposons

.
Alors l'équation nous donne :
}*(ab)^n=2^{2\beta*m}(2^{2(\alpha -\beta)}*a^2+b^2)^m)
or

donc
}a^2+b^2)
est impair.
d'où, en identifiant les puissances de 2 :
=2\beta*m)
donc
 > 2n*\beta)
donc

Or,
^2 \geq 0)
donc, si

,

donc
^n = (x^2+y^2)^m \geq (2xy)^m > (2xy)^n \geq (xy)^n)
Absurde.
donc

(un raisonnement identique donne

impossible)
donc


Donc :
^n=2^m(a^2+b^2)^m)
Notons que

n'a pas de solutions avec a et b impairs, donc
=1)
d'où, en identifiant les valuations 2-adique des deux côtés :

donc

divise m mais par gauss m divise

donc


On obtient alors facilement l'ensemble des solutions :
=(2^{\alpha}, 2^{\alpha}, 2\alpha, 2\alpha +1))