OFFSET
0,2
PROG
(PARI) {a(n)=local(B=[1]); if(n==0, 1, for(k=1, n, t=3*a(k-1); B=concat(B, t); B[ #B]=t+1-denominator(Vec(Ser(B)^(1/2))[ #B]) )); B[n+1]}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Sep 07 2007
STATUS
approved