OFFSET
0,2
COMMENTS
In general, for k > 0, if g.f. = product_{m>=1} (1+q^m)^k, then a(n) ~ k^(1/4) * exp(Pi * sqrt(k*n/3)) / (2^((k+3)/2) * 3^(1/4) * n^(3/4)) * (1 + (Pi*k^(3/2) / (48*sqrt(3)) - 3^(3/2) / (8*Pi*sqrt(k))) / sqrt(n)). - Vaclav Kotesovec, Mar 05 2015, extended Jan 16 2017
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..5000
FORMULA
a(n) ~ exp(Pi * 4 * sqrt(2*n/3)) / (65536 * 6^(1/4) * n^(3/4)). - Vaclav Kotesovec, Mar 05 2015
MATHEMATICA
nmax=50; CoefficientList[Series[Product[(1+q^m)^32, {m, 1, nmax}], {q, 0, nmax}], q] (* Vaclav Kotesovec, Mar 05 2015 *)
PROG
(PARI) m=50; q='q+O('q^m); Vec(prod(n=1, m, (1+q^n)^32)) \\ G. C. Greubel, Mar 20 2018
(Magma) Coefficients(&*[(1+x^m)^32:m in [1..40]])[1..40] where x is PolynomialRing(Integers()).1; // G. C. Greubel, Mar 20 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Terms a(19) onward added by G. C. Greubel, Mar 20 2018
STATUS
approved