OFFSET
0,3
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
FORMULA
G.f.: Sum_{j>=1} (x^j/(1-x^j)) * Product_{k>=1} ((1+x^k)/(1-x^k))^k.
a(n) ~ A156616(n) * 3*n^(1/3) * (log(3*n^(1/3)/(2*r)) + gamma) / (2*r), where r = 3*(7*zeta(3))^(1/3) / 2^(4/3) and gamma is the Euler-Mascheroni constant A001620.
a(n) ~ exp(1/12 + 3*(7*zeta(3))^(1/3) * (n^(2/3) / 2^(4/3))) * (log((2*n/(7*zeta(3)))^(1/3)) + gamma) / (A * 2^(4/9) * sqrt(3*Pi) * (7*zeta(3))^(5/36) * n^(13/36)), where A is the Glaisher-Kinkelin constant (A074962).
MATHEMATICA
nmax = 40; CoefficientList[Series[Sum[x^j/(1-x^j), {j, 1, nmax}] * Product[((1+x^k)/(1-x^k))^k, {k, 1, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Jun 14 2026
STATUS
approved
