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