OFFSET
0,3
COMMENTS
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