OFFSET
1,2
COMMENTS
Inverse Möbius transform of A002417.
LINKS
N. J. A. Sloane, Transforms
FORMULA
G.f.: Sum_{k>=1} x^k*(1 + 3*x^k)/(1 - x^k)^5.
G.f.: Sum_{k>=1} k*A000292(k)*x^k/(1 - x^k).
L.g.f.: -log(Product_{k>=1} (1 - x^k)^A000292(k)) = Sum_{n>=1} a(n)*x^n/n.
Dirichlet g.f.: (zeta(s-4) + 3*zeta(s-3) + 2*zeta(s-2))*zeta(s)/6.
a(n) = (2*sigma_2(n) + 3*sigma_3(n) + sigma_4(n))/6.
a(n) = Sum_{d|n} A002417(d).
Sum_{k=1..n} a(k) ~ Zeta(5) * n^5 / 30. - Vaclav Kotesovec, Feb 02 2019
MATHEMATICA
Table[Sum[d Binomial[d + 2, 3], {d, Divisors[n]}], {n, 40}]
nmax = 40; Rest[CoefficientList[Series[Sum[x^k (1 + 3 x^k)/(1 - x^k)^5, {k, 1, nmax}], {x, 0, nmax}], x]]
Table[(2 DivisorSigma[2, n] + 3 DivisorSigma[3, n] + DivisorSigma[4, n])/6, {n, 40}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 14 2018
STATUS
approved