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