OFFSET
1,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000 (terms 1..153 from Sean A. Irvine)
PROG
(PARI) lista(n)={my(v=vector(n), M=Map()); v[2]=9; mapput(M, 3, 1); for(i=3, n, my(s=v[i-1]+v[i-2]); fordiv(s, d, my(t=(d-s/d)/2); if(t>0&&!frac(t)&&!mapisdefined(M, t), mapput(M, t, 1); v[i]=t^2; break)); if(!v[i], v=v[1..i-1]; break)); v} \\ Andrew Howroyd, Feb 17 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Sep 20 2002
EXTENSIONS
More terms from David Wasserman, Jan 16 2005
STATUS
approved