OFFSET
0,2
COMMENTS
In general, if g.f. = Product_{k>=1} (1+x^k)^(m*k+c), m > 0, then a(n) ~ (m*Zeta(3))^(1/6) * exp(-c^2 * Pi^4 / (1296*m*Zeta(3)) + (c * Pi^2 * n^(1/3)) / (2^(5/3) * 3^(4/3) * (m*Zeta(3))^(1/3)) + 3^(4/3) * (m*Zeta(3))^(1/3) * n^(2/3) / 2^(4/3)) / (2^(m/12 + c/2 + 2/3) * 3^(1/3) * sqrt(Pi) * n^(2/3)). - Vaclav Kotesovec, Mar 08 2015
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..1000
FORMULA
a(n) ~ Zeta(3)^(1/6) * exp(-Pi^4/(972*Zeta(3)) + Pi^2 * n^(1/3) / (2^(2/3) * 3^(5/3) * Zeta(3)^(1/3)) + 3^(5/3)/2^(4/3) * Zeta(3)^(1/3) * n^(2/3)) / (2^(23/12) * 3^(1/6) * sqrt(Pi) * n^(2/3)), where Zeta(3) = A002117.
MATHEMATICA
nmax=50; CoefficientList[Series[Product[(1+x^k)^(3*k+2), {k, 1, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Mar 07 2015
STATUS
approved