OFFSET
0,3
COMMENTS
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
FORMULA
a(n) ~ exp(Pi*sqrt(n)) * (2*gamma + log(4*n/Pi^2)) / (8*Pi*sqrt(n)), where gamma is the Euler-Mascheroni constant A001620.
MATHEMATICA
nmax = 40; CoefficientList[Series[Sum[x^k/(1-x^k), {k, 1, nmax}] * Product[(1+x^k)/(1-x^k), {k, 1, nmax}], {x, 0, nmax}], x]
PROG
(PARI) my(N=44, q='q+O('q^N)); Vec( prod(k=1, N, (1+q^k)/(1-q^k)) * sum(k=1, N, 1*q^k/(1-q^k)) ) \\ Joerg Arndt, Jun 18 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, May 25 2018
STATUS
approved