OFFSET
1,2
COMMENTS
Inverse Moebius transform of centered triangular numbers (A005448).
FORMULA
G.f.: Sum_{k>=1} (3*k*(k - 1)/2 + 1) * x^k/(1 - x^k).
a(n) = 3 * (sigma_2(n) - sigma_1(n))/2 + d(n).
MATHEMATICA
nmax = 52; CoefficientList[Series[Sum[x^k (1 - x^(3 k))/(1 - x^k)^4, {k, 1, nmax}], {x, 0, nmax}], x] // Rest
Table[3 (DivisorSigma[2, n] - DivisorSigma[1, n])/2 + DivisorSigma[0, n], {n, 1, 52}]
PROG
(PARI) a(n)={sumdiv(n, d, 3*d*(d-1)/2 + 1)} \\ Andrew Howroyd, Aug 14 2019
(PARI) a(n)={3*(sigma(n, 2) - sigma(n))/2 + numdiv(n)} \\ Andrew Howroyd, Aug 14 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Aug 14 2019
STATUS
approved