OFFSET
0,3
COMMENTS
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
FORMULA
G.f.: Sum_{j>=1} (j*x^j/(1-x^j))*Product_{k>=1} (1+x^k).
a(n) ~ 2*n*A000009(n) ~ exp(Pi*sqrt(n/3)) * n^(1/4) / (2*3^(1/4)).
MATHEMATICA
Table[Sum[DivisorSigma[1, k] * PartitionsQ[n-k], {k, 1, n}], {n, 0, 50}]
nmax = 50; CoefficientList[Series[Sum[j*x^j/(1-x^j), {j, 1, nmax}]*Product[1+x^k, {k, 1, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Sep 25 2016
STATUS
approved