|
|
A114050
|
|
x-values in the solution to x^2 - 22*y^2 = 1.
|
|
3
|
|
|
1, 197, 77617, 30580901, 12048797377, 4747195585637, 1870383011943601, 736926159510193157, 290347036464004160257, 114395995440658128948101, 45071731856582838801391537, 17758147955498197829619317477
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
a(n) = 394*a(n-1) - a(n-2); a(1)=1, a(2)=197. - Benoit Cloitre, Feb 03 2006
|
|
MATHEMATICA
|
LinearRecurrence[{394, -1}, {1, 197}, 20] (* Harvey P. Dale, Nov 03 2011 *)
|
|
PROG
|
(PARI) g(n, k=22) = for(y=0, n, x=k*y^2+1; if(issquare(x), print1(sqrtint(x), ", ")))
(PARI) a=vector(12); a[1]=1; a[2]=197; for(i=3, #a, a[i]=394*a[i-1]-a[i-2]); a \\ Benoit Cloitre, Feb 03 2006
(Magma) I:=[1, 197]; [n le 2 select I[n] else 394*Self(n-1)-Self(n-2): n (Magma) in [1..40]]; // Vincenzo Librandi, Dec 21 2011
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Offset changed from 0 to 1, g.f. and formula adapted by Vincenzo Librandi, Dec 21 2011
|
|
STATUS
|
approved
|
|
|
|