OFFSET
1,3
COMMENTS
The fundamental solution of the Pell equation x^2 - 5*(n*y)^2 = 1, is the smallest solution of x^2 - 5*y^2 = 1 satisfying y = 0 (mod n).
REFERENCES
Michael J. Jacobson, Jr. and Hugh C. Williams, Solving the Pell Equation, Springer, 2009, pages 1-17.
LINKS
A.H.M. Smeets, Table of n, a(n) for n = 1..20000
H. W. Lenstra Jr., Solving the Pell Equation, Notices of the AMS, Vol.49, No.2, Feb. 2002, pp. 182-192.
FORMULA
MATHEMATICA
b[n_] := b[n] = Switch[n, 0, 0, 1, 4, _, 18 b[n - 1] - b[n - 2]];
a[n_] := For[k = 1, True, k++, If[Mod[b[k], n] == 0, Return[k]]];
a /@ Range[100] (* Jean-François Alcover, Nov 16 2019 *)
PROG
(Python)
xf, yf = 9, 4
x, n = 2*xf, 0
while n < 20000:
n = n+1
y1, y0, i = 0, yf, 1
while y0%n != 0:
y1, y0, i = y0, x*y0-y1, i+1
print(n, i)
CROSSREFS
KEYWORD
nonn
AUTHOR
A.H.M. Smeets, Jan 15 2018
STATUS
approved