OFFSET
0,2
COMMENTS
In general, if g.f. = Product_{k>=1} ((1 + x^(2*k-1))/(1 - x^(2*k-1)))^(c2*k^2 + c1*k + c0) and c2 > 0, then a(n) ~ exp(Pi*sqrt(2) * c2^(1/4) * n^(3/4) / 3 + 7*(c1+c2) * Zeta(3) * sqrt(n) / (2*sqrt(c2) * Pi^2) + (Pi*(4*c0 + 2*c1 + c2) / (8*sqrt(2) * c2^(1/4)) - 49*(c1+c2)^2 * Zeta(3)^2 / (2^(3/2) * c2^(5/4) * Pi^5)) * n^(1/4) - (7*c0 + 21*c1/4 + c2 + 7*c0*c1/c2 + 7*c1^2/(2*c2)) * Zeta(3) / (4*Pi^2) + 22411*(c1+c2)^3 * Zeta(3)^3 / (196 * c2^2 * Pi^8) - (c1+c2)/24) * A^((c1+c2)/2) * (n^((c1+c2)/96 - 5/8) / (2^(c0/2 + (11*c1 + 5*c2)/48 + 9/4) * Pi^((c1+c2)/24) * c2^((c1+c2)/96 - 1/8))), where A is the Glaisher-Kinkelin constant A074962.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..5000
FORMULA
a(n) ~ exp(sqrt(2)*Pi * n^(3/4)/3 + 7*Zeta(3) * sqrt(n) / (2*Pi^2) + (Pi / (8*sqrt(2)) - 49*Zeta(3)^2 / (2^(3/2) * Pi^5)) * n^(1/4) + 22411*Zeta(3)^3 / (196*Pi^8) - Zeta(3)/(4*Pi^2) - 1/24) * sqrt(A) / (2^(113/48) * Pi^(1/24) * n^(59/96)), where A is the Glaisher-Kinkelin constant A074962.
MATHEMATICA
nmax = 50; CoefficientList[Series[Product[((1+x^(2*k-1))/(1-x^(2*k-1)))^(k^2), {k, 1, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Nov 08 2017
STATUS
approved