OFFSET
1,2
COMMENTS
Not multiplicative: a(3)*a(7) <> a(21), for example.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
FORMULA
G.f.: Sum_{n>0} n*x^n/(1-x^(4*n)). - Vladeta Jovovic, Nov 14 2002
G.f.: Sum_{k>0} x^(4*k-3) / (1 - x^(4*k-3))^2. - Seiichi Manyama, Jun 29 2023
from Amiram Eldar, Nov 05 2023: (Start)
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = A222183. (End)
MAPLE
A050460 := proc(n)
a := 0 ;
for d in numtheory[divisors](n) do
if (n/d) mod 4 = 1 then
a := a+d ;
end if;
end do:
a;
end proc:
seq(A050460(n), n=1..40) ; # R. J. Mathar, Dec 20 2011
MATHEMATICA
a[n_] := DivisorSum[n, Boole[Mod[n/#, 4] == 1]*#&]; Array[a, 70] (* Jean-François Alcover, Dec 01 2015 *)
PROG
(PARI) a(n)=sumdiv(n, d, if(n/d%4==1, d)) \\ Charles R Greathouse IV, Dec 04 2013
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Dec 23 1999
STATUS
approved