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