OFFSET
0,3
FORMULA
a(n) = [x^n] Product_{i=0..n-1} (1 + x^(2^i))^(n-i) for n>=0.
MAPLE
A131824 := proc(n) local f, i ; f := 1 ; if n = 1 then RETURN(1) ; fi ; for i from 0 to n-1 do f := taylor(f*(1+x^(2^i))^(n-i), x=0, n+1) ; od; coeftayl(f, x=0, n) ; end: seq(A131824(n), n=0..27) ; # R. J. Mathar, Jul 24 2007
PROG
(PARI) a(n)=if(n==0, 1, polcoeff(prod(i=0, n-1, (1 + x^(2^i) +x*O(x^n))^(n-i)), n))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 19 2007
EXTENSIONS
More terms from R. J. Mathar, Jul 24 2007
STATUS
approved