OFFSET
0,7
COMMENTS
In general, if m >= 1 and g.f. = Product_{k>=1} 1/(1-x^(m*k-m+1))^(m*k-m+1), then a(n, m) ~ exp(c*m + 3 * 2^(-2/3) * m^(-1/3) * Zeta(3)^(1/3) * n^(2/3)) * 2^(19*m/36 + 1/(6*m) - 1) * m^(17*m/36 + 5/(6*m) - 3/2) * Pi^(m/2 - 1) * Zeta(3)^(1/(6*m) + m/36) / (sqrt(3) * Gamma(1/m)^(m-1) * n^(1/2 + 1/(6*m) + m/36)), where c = Integral_{x=0..infinity} exp((2*m-1)*x) / (x*(exp(m*x) - 1)^2) + (1/12 - (m-1)^2/(2*m^2))/(x*exp(x)) - 1/(m^2*x^3) - (m-1)/(m^2*x^2) dx. - Vaclav Kotesovec, Apr 17 2017
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..2000
FORMULA
a(n) ~ 2^(301/180) * 5^(37/36) * Pi^(3/2) * Zeta(3)^(31/180) * exp(5*c + 3 * 2^(-2/3) * 5^(-1/3) * Zeta(3)^(1/3) * n^(2/3)) / (sqrt(3) * Gamma(1/5)^4 * n^(121/180)), where c = Integral_{x=0..inf} ((-71/(exp(x)*300) + 1/(exp(x)*(1 - exp(-5*x))^2) - 1/(25*x^2) - 4/(25*x))/x) dx = 0.186382690624752630391368364629918483384424086341764409146923686... - Vaclav Kotesovec, Apr 16 2017
MATHEMATICA
nmax = 50; CoefficientList[Series[Product[1/(1-x^(5*k-4))^(5*k-4), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Apr 16 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Apr 15 2017
STATUS
approved