$$8y^2=4x^2+4x=(2x+1)^2-1\implies z^2-8y^2=1\text{ where } z=2x+1$$
Utilizando este ,
$$z_k\pm2\sqrt2y_k=(z_1\pm2\sqrt2\cdot y_1)^k$$
Por observación, $z_1=\pm3,y_1=\pm1$ como $z^2-8y^2=1=3^2-8\cdot1$
Si tomamos $z_1=3,y_1=1$ $$z_k+2\sqrt2y_k=(3+2\sqrt2)^k \text{ and } z_k-2\sqrt2y_k=(3-2\sqrt2)^k$$
Resuelva ahora $z_k,y_k$
Del mismo modo, para las demás combinaciones
$z_1=-3,y_1=-1$
$z_1=3,y_1=-1$
et $z_1=-3,y_1=1$
A partir de la relación de recurrencia en Wikipedia si $x_1,y_1,x_k,y_k$ son números enteros, también lo serán $x_{k+1},y_{k+1}$