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