Salut,
Perso, vu la tête du truc, j'essayerais bien de copier la méthode classique concernant les équation du type xy=3x+5y+7 qu'on écrit xy-3x-5y=7 puis où on dit que xy-3x-5y, c'est "presque" (x-5)(y-3) sauf qu'il manque le +15 à la fin donc l'équation équivaut à xy-3x-5y+15=7+15 c'est à dire (x-5)(y-3)=22 et y'a plus qu'a chercher les factorisation possible de 22.
Sauf que là, il y a du xy
z avec le z qui fait c..., mais on peut faire comme si c'était un paramètre pour voir ce que ça donne. Comme pour avoir du xy=... faut diviser par z, on commence par regarder le cas particulier z=0 qui donne clairement l'unique solution x=0 ; y=0 (car les entiers doivent être positifs; dans Z il y aurait d'autres solutions).
Pour z non nul on a alors :
\ \Leftrightarrow\ <br />xy\!-\!\dfrac{4}{z}x\!-\dfrac{4}{z}y\!=\!4\ \Leftrightarrow\ <br />\Big(x\!-\!\dfrac{4}{z}\Big)\Big(y\!-\dfrac{4}{z}\Big)\!=\!4\!+\!\dfrac{16}{z^2}\ \Leftrightarrow\ <br />(zx\!-\!4)(zy\!-\!4)\!=\!4z^2\!+\!16)
qui permet facilement de trouver les solutions pour un

fixé (il suffit de chercher les factorisations de

)
Or, une fois trouvé "à la main" les solutions pour

,

,

, c'est fini vu que, par symétrie de l'énoncé, on connait en fait toutes les solutions où une des trois variables est

or, si elles étaient toutes les trois

, on aurait
(zy\!-\!4)\!\geq\!(4z\!-\!4)^2\!=\!16z^2\!-\!32z\!+\!16\!=\!4z^2\!+\!16\!+\!4z(3z\!-\!8)>4z^2\!+\!16)