Prueba para todos los no cuadráticos $n$
Lema: Ecuación de Pell $x^2-n y^2 = 1$ con $n$ al no ser un cuadrado perfecto, tiene infinitas soluciones tales que $x$ es impar, $y$ es par y $x\equiv1$ (mod $2n$ ).
Prueba: Es un hecho bien conocido que la ecuación de Pell con no cuadrática $n$ tiene un número inifinito de soluciones. Escoge una solución cualquiera $(x_1,y_1)$ . Cabe destacar que $x_1$ y $y_1$ debe ser coprima, así como $x_1$ y $n$ .
Ahora calcula:
$$x_2=x_1^2+ny_1^2,\quad y_2=2x_1y_1$$
Se puede demostrar fácilmente que $(x_2,y_2)$ es también una solución de la misma ecuación de Pell. Evidentemente, $y_2$ está en paz.
Si sustituye $x_1^2=ny_1^2+1$ en la expresión para $x_2$ lo consigues:
$$x_2=1+2ny_1^2\implies x_2\equiv1\space (\text{mod}\space 2n)$$
Esto también demuestra que $x_2$ tiene que ser impar (lo que tiene mucho sentido porque las soluciones de la ecuación de Pell son siempre coprimas y $y_2$ es par).
Se pueden construir más soluciones de la ecuación de Pell de la misma manera y todas satisfacen los criterios del lema. Por lo tanto, no hay una sola solución de este tipo. En realidad hay infinitas.
Fin de la demostración del lema
Volvemos a la ecuación original (mismo enfoque que AQUÍ ) :
$$\frac{a^2+b}{a+b^2}=n$$
se puede reescribir como
$$u^2-nv^2=1-n^3 \tag1$$
donde:
$$u=2nb-1,\quad v=2a-n$$
Tome $x,y$ tal que:
$$x^2-ny^2=1\tag2$$
Puede demostrar fácilmente que $(-x+y n^2)$ y $(-y+nx)$ satisfacer (1):
\begin{align*}(-x+ yn^2)^2-n(-y+nx)^2&=x^2-2xyn^2+y^2n^4-ny^2+2xyn^2-x^2n^3\\&=(x^2-ny^2)+n^3(ny^2-x^2)\\&=1-n^3.\end{align*}
Esto demuestra que:
$$u=-x+yn^2=2nb-1$$
$$v=-y+xn=2a-n$$
...representan una solución de $(1)$ .
Por lo tanto,
$$a=\frac{(x+1)n-y}{2},\space b=\frac{yn^2-(x-1)}{2n}.\tag3$$
Según nuestro lema La ecuación de Pell tiene infinitas soluciones $x,y$ tal que $x$ es impar, $y$ es par y $x\equiv1$ (mod $2n$ ). Reemplaza estas soluciones en (3) y obviamente obtendrás infinitos valores enteros para $a,b$ .
Fin de la prueba para todos los no cuadrados $n$ .
El siguiente script simple de Mathematica generará una sola $a,b$ para un tamaño no cuadrado $n$ muy rápido (sigue la prueba, palabra por palabra):
ABPair[n_] := Module[
{x, y, a, b, a1, b1, a2, b2},
pellSolutions = Solve[x^2 - n y^2 == 1, {x, y}, Integers] /. C[1] -> 1;
pellSolutions = {x, y} /. pellSolutions;
{a1, b1} = First[Select[pellSolutions, #[[1]] > 0 && #[[2]] > 0 &, 1]];
{a2, b2} = If[Mod[a1, 2 n] == 1 && Mod[b1, 2] == 0, {a1, b1}, {a1^2 + n b1^2, 2 a1 b1}];
a = (n (a2 + 1) - b2)/2;
b = (b2 n^2 - a2 + 1)/(2 n);
{a, b, (a^2 + b)/(b^2 + a)}
];
Por ejemplo:
ABPair[5613]
{60584278414870816497213, 808653403020126409200, 5613}
El tercer número es sólo una comprobación de que los números calculados son válidos. En otras palabras:
$$\frac{60584278414870816497213^2+808653403020126409200}{60584278414870816497213+808653403020126409200^2}=5613$$
El script es rápido como un rayo incluso para $n$ con 12 dígitos:
ABPair[561044335534]
Ver la solución de Sil para la cuadrática $n$ . Caso cerrado :)