OFFSET
0,3
PROG
(PARI) {a(n)=local(B=vector(n+1, k, if(k==1, 1, sum(i=0, k-2, 2^(i+1)*binomial(2*k-2, i)*binomial(k-1, i+1))/(k-1)))); return(polcoeff(Ser(B)^2+x*Ser(B)^4+x*O(x^n), n))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 01 2005
STATUS
approved