OFFSET
0,4
FORMULA
a(n) ~ exp(Pi*sqrt(n)) * (log(4*n/Pi^2) + 2*gamma)^2 / (8*Pi^2), where gamma is the Euler-Mascheroni constant A001620.
MATHEMATICA
nmax = 40; CoefficientList[Series[Sum[x^j/(1-x^j), {j, 1, nmax}]^2 * Product[(1+x^k)/(1-x^k), {k, 1, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Jun 16 2026
STATUS
approved
