OFFSET
0,5
COMMENTS
In general, if m > 1 and g.f. = Product_{k>=1} 1/(1-x^(m*k-1))^(m*k-1), then a(n, m) ~ exp(c*m + 3 * 2^(-2/3) * m^(-1/3) * Zeta(3)^(1/3) * n^(2/3)) * (2*Zeta(3))^(1/(6*m) + m/36) / (sqrt(3) * Gamma(1 - 1/m) * m^(1/2 - 5/(6*m) + m/36) * n^(1/2 + 1/(6*m) + m/36)), where c = Integral_{x=0..infinity} exp((m+1)*x) / (x*(exp(m*x)-1)^2) + (1/12 - 1/(2*m^2))/(x*exp(x)) - 1/(m^2*x^3) - 1/(m^2*x^2) dx. - Vaclav Kotesovec, Apr 17 2017
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..10000
FORMULA
a(n) ~ exp(5*c + 3*2^(-2/3)*5^(-1/3)*Zeta(3)^(1/3)*n^(2/3)) * (2*Zeta(3))^(31/180) / (sqrt(3) * 5^(17/36) * Gamma(4/5) * n^(121/180)), where c = Integral_{x=0..inf} ((19/(exp(x)*300) + 1/(exp(4*x)*(1-exp(-5*x))^2) - 1/(25*x^2) - 1/(25*x))/x) dx = -0.12699586713882325294527057473113580561183418857868946729897216431919... - Vaclav Kotesovec, Apr 15 2017
MATHEMATICA
nmax = 100; CoefficientList[Series[Product[1/(1-x^(5*k-1))^(5*k-1), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Apr 15 2017 *)
PROG
(PARI) x='x+O('x^100); Vec(prod(k=0, 100, 1/(1 - x^(5*k + 4))^(5*k + 4))) \\ Indranil Ghosh, Apr 15 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Apr 15 2017
STATUS
approved