OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..400
EXAMPLE
G.f.: A(x) = 1 + 3*x + 12*x^2 + 64*x^3 + 354*x^4 + 2160*x^5 + 13518*x^6 +...
where
A(x) = ((1+x)*(1+3*x^2)*(1+12*x^3)*(1+64*x^4)*(1+354*x^5)*...)^3.
Related expansion:
A(x)^(1/3) = 1 + x + 3*x^2 + 15*x^3 + 76*x^4 + 454*x^5 + 2742*x^6 +...
MAPLE
A:= proc(n) option remember; local i, p, q; if n=0 then 1 else
p, q:= A(n-1), 1; for i from 0 to n-1 do q:= convert(
series(q*(1+coeff(p, x, i)*x^(i+1))^3, x, n+1), polynom)
od: q fi
end:
a:= n-> coeff(A(n), x, n):
seq(a(n), n=0..30); # Alois P. Heinz, Aug 01 2013
MATHEMATICA
a[n_] := a[n] = SeriesCoefficient[Product[(1+a[k] x^(k+1))^3, {k, 0, n-1}], {x, 0, n}];
a /@ Range[0, 30] (* Jean-François Alcover, Nov 20 2020 *)
PROG
(PARI) {a(n) = polcoeff(prod(k=0, n-1, (1+a(k)*x^(k+1)+x*O(x^n)))^3, n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 02 2012
STATUS
approved