OFFSET
1,10
LINKS
FORMULA
a(n) = Sum_{d|n, d<n} [1 == d mod 4]*d.
G.f.: Sum_{k>=1} (4*k-3) * x^(8*k-6) / (1 - x^(4*k-3)). - Ilya Gutkovskiy, Apr 14 2021
Sum_{k=1..n} a(k) = c * n^2 + O(n*log(n)), where c = Pi^2/48 - 1/8 = 0.0806167... . - Amiram Eldar, Nov 27 2023
MATHEMATICA
a[n_] := DivisorSum[n, # &, # < n && Mod[#, 4] == 1 &]; Array[a, 100] (* Amiram Eldar, Nov 27 2023 *)
PROG
(PARI) A293901(n) = sumdiv(n, d, (d<n)*(1==(d%4))*d);
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Antti Karttunen, Oct 19 2017
STATUS
approved