Dada la ecuación $$ 2 n^2 - 1 = ( m - n )^2. $$
A partir de esto podemos escribir $$ 2 m n = m^2 - n^2 + 1. $$
Dejemos que $$ x = m^2 - n^2, y = 2 m n, z = m^2 + n^2 \Rightarrow x^2 + y^2 = z^2. $$
Pero $$ y = x + 1, $$
por lo que obtenemos $$ x^2 + ( x + 1 )^2 = z^2, $$
que puede escribirse como $$ 2 z^2 - 1 = ( 2 x + 1 )^2, $$
o $$ 2 \big( \underbrace{m^2 + n^2}_{n'} \big)^2 - 1 = \big( \underbrace{ 2 \big[ m^2 - n^2 \big] + 1 }_{m'-n'} \big)^2. $$
Por lo tanto, a partir de $$ \big( m, n \big) $$
podemos "generar" el siguiente par $(m',n')$ utilizando $$ \big( m', n' \big) = \big( 3 m^2 - n^2 + 1, m^2 + n^2 \big). $$
Tenga en cuenta que $m=2$ y $n=1$ rinde $$ 2 \cdot 1^2 - 1 = ( 2 - 1 )^2 $$
Ahora podemos "generar" los siguientes pares...
$$ \big( m', n' \big) = \big( 3 \cdot 2^2 - 1^2 + 1, 2^2 + 1^2 \big) = \big( 12 , 5 \big). $$
$$ \big( m'', n'' \big) = \big( 3 \cdot 12^2 - 5^2 + 1, 12^2 + 5^2 \big) = \big( 408 , 169 \big). $$
$$ \big( m''', n''' \big) = \big( 3 \cdot 408^2 - 169^2 + 1, 408^2 + 169^2 \big) = \big( 470832, 195025 \big). $$
El número que buscamos son los números $n$ la primera dada por $$ n = 1\\ n = 5\\ n = 169\\ n = 195025 $$
Esto no son todos los números, ya que el es un segundo método para generar pares $(m',n')$ .
Todavía estoy trabajando en el segundo método.
El método general viene dado por
$$ 2 x_n^2 - 1 = y_n^2, $$
donde ambos $x_n$ y $y_n$ son números enteros, entonces $$ x_n = \frac{7 + 5 \sqrt{2}}{2 \sqrt{2}} \Big( 3 + 2 \sqrt{2} \Big)^{n-1} -\frac{7 - 5 \sqrt{2}}{2 \sqrt{2}} \Big( 3 - 2 \sqrt{2} \Big)^{n-1}. $$
Los resultados son $$ \begin{array}{l|l} n & x_n\\ \hline 0 & 1\\ 1 & 5\\ 2 & 29\\ 3 & 169\\ 4 & 985\\ 5 & 5741\\ 6 & 33461\\ 7 & 195025\\ 8 & 1136689\\ 9 & 6625109\\ 10 & 38613965 \end{array} $$
(necesito publicar cómo derivarlo, pero es largo)
La idea básica para derivarlo:
Dado $$ 2 x_n^2 - 1 = y_n^2. $$
Tenga en cuenta que $$ \left( \begin{array}{c} x_{n+1} \\ y_{n+1} \end{array} \right) = \left( \begin{array}{cc} 3 & 2 \\ 4 & 3 \end{array} \right) \left( \begin{array}{c} x_n \\ y_n \end{array} \right) $$
y $$ \left( \begin{array}{c} x_0 \\ y_0 \end{array} \right) = \left( \begin{array}{c} 1 \\ 1 \end{array} \right). $$
Así que $$ \left( \begin{array}{c} x_n \\ y_n \end{array} \right) = \left( \begin{array}{cc} 3 & 2 \\ 4 & 3 \end{array} \right)^n \left( \begin{array}{c} 1 \\ 1 \end{array} \right). $$
Dejemos que $\chi$ sea la traza de la matriz y $\Delta$ el determinante, entonces los valores propios vienen dados por $\lambda_\pm = \chi/2 \pm \sqrt{\chi^2/4 - \Delta}$ .
Entonces $$ \left( \begin{array}{cc} 3 & 2 \\ 4 & 3 \end{array} \right)^n = \frac{\lambda_+^n - \lambda_-^n}{\lambda_+ - \lambda_-} \left( \begin{array}{cc} 3 & 2 \\ 4 & 3 \end{array} \right) - \lambda_+ \lambda_- \frac{\lambda_+^{n-1} - \lambda_-^{n-1}}{\lambda_+ - \lambda_-} \left( \begin{array}{cc} 1 & 0 \\ 0 & 1 \end{array} \right) $$
De esto se deduce que $$ x_n = 5 \frac{\lambda_+^n - \lambda_-^n}{\lambda_+ - \lambda_-} - \lambda_+ \lambda_- \frac{\lambda_+^{n-1} - \lambda_-^{n-1}}{\lambda_+ - \lambda_-} $$
y $$ y_n = 7 \frac{\lambda_+^n - \lambda_-^n}{\lambda_+ - \lambda_-} - \lambda_+ \lambda_- \frac{\lambda_+^{n-1} - \lambda_-^{n-1}}{\lambda_+ - \lambda_-}. $$
Si lo resolvemos, obtenemos $$ x_n = \frac{7 + 5 \sqrt{2}}{2 \sqrt{2}} \Big( 3 + 2 \sqrt{2} \Big)^{n-1} -\frac{7 - 5 \sqrt{2}}{2 \sqrt{2}} \Big( 3 - 2 \sqrt{2} \Big)^{n-1}. $$
(He escrito este post deprisa, así que perdonadme por algunas erratas :)