OFFSET
0,2
COMMENTS
a(n) ~ 2^n*c where c=0.683268303731665578137197731736227240453607...
PROG
(PARI) {a(n)=local(B=[1], t); if(n==0, 1, for(k=1, n, t=2*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