OFFSET
1,2
COMMENTS
This equation can also be written as (2*a(n) - b(n))^2 - D(n)*b(n)^2 = +4 or -4 with D(n) := A077425(n) = 1 + 4*G(n).
This is from Perron's table (see reference p. 108, for n = 1..28) which gives the minimal x,y values which solve the above mentioned Diophantine equations.
REFERENCES
O. Perron, "Die Lehre von den Kettenbruechen, Bd.I", Teubner, 1954, 1957 (Sec. 30, Satz 3.35, p. 109 and table p. 108).
LINKS
FORMULA
MATHEMATICA
g[n_] := Ceiling[Sqrt[n]] + n - 1; r[n_] := Reduce[an > 0 && bn > 0 && (an ^2 - an*bn - g[n]*bn^2 == 1 || an^2 - an*bn - g[n]*bn^2 == - 1), {an, bn}, Integers] /. C -> c; ab[n_] := DeleteCases[ Flatten[ Table[{an, bn} /. {ToRules[r[n]]} // Simplify, {c[1], 0, 1}], 1], an | bn]; a[n_] := a[n] = Min[ab[n][[All, 1]]]; Table[Print[{n, a[n]}]; a[n], {n, 1, 62}] (* Jean-François Alcover, Oct 04 2012 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Wolfdieter Lang, Nov 29 2002
EXTENSIONS
More terms from Max Alekseyev, Feb 06 2010
a(9), a(33), a(54) corrected (after notice by Jean-François Alcover); a(58) through a(62) added. - Wolfdieter Lang, Oct 04 2012
STATUS
approved