OFFSET
1,5
FORMULA
G.f.: Sum_{k>0} k*(k+1)/2 * x^(3*k-1) / (1 - x^(3*k-1)).
a(n) = Sum_{d|n, d==2 mod 3} binomial((d+1)/3+1,2).
MATHEMATICA
a[n_] := DivisorSum[n, Binomial[(#+1)/3+1, 2] &, Mod[#, 3] == 2 &]; Array[a, 100] (* Amiram Eldar, Jun 30 2023 *)
PROG
(PARI) a(n) = sumdiv(n, d, (d%3==2)*binomial((d+1)/3+1, 2));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jun 30 2023
STATUS
approved