On the Diophantine equation \((x^2-1)(y^2-1) = [(\tfrac 12 (y-x))^2-1]^2\) (Q2394236)

L'identité \[ (x^2-1)(y^2-1)-[(\tfrac 12 (y-x))^2 -1]^2=- \tfrac 1{16}(x+y)^2(x^2- 6xy+y^2 +8) \] permet à l'auteur de se ramener à la résolution (classique) de \(x^2-2z^2=1\), dont les solutions positives correspondent à toutes les solutions positives \((0<x<y)\) de l'équation initiale (sauf \(x=1\), \(y=3\)), en posant \(y=4z+3x\). L'A. remarque que \(x_0=1, x_1=3, x_{n+2}= 6x_{n+1}-x_n\) donnent toutes les solutions impaires de \((x^2-1)(y^2-1)=(z^2-1)^2\), en posant \(z=\tfrac 12 (x_{n+1}-x_n)\). Le problème des solutions non impaires reste à voir, en dehors des solutions de \textit{K. Szymiczek} [Elem. Math. 22, 37--38 (1967; Zbl 0171.00802)] (en particulier \(x=4, y=31\) et \(z=11\)).