OFFSET
1,10
COMMENTS
In general, for d > 0, if g.f. = Sum_{k>=1} k * x^(k*(d*k - d + 2)/2) / (1 - x^k), then Sum_{k=1..n} a(k) ~ 2^(3/2) * n^(3/2) / (3*sqrt(d)).
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..10000
FORMULA
Sum_{k=1..n} a(k) ~ n^(3/2) / 3.
MATHEMATICA
Table[Sum[If[n > 4*k*(k-1), k, 0], {k, Divisors[n]}], {n, 1, 100}]
nmax = 100; Rest[CoefficientList[Series[Sum[k*x^(k*(4*k - 3))/(1 - x^k), {k, 1, Sqrt[nmax/4] + 1}], {x, 0, nmax}], x]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Oct 23 2024
STATUS
approved