OFFSET
1,5
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
MAPLE
A050451 := proc(n) local a, d; a := 0 ; for d in numtheory[divisors](n) do if d mod 4 = 1 then a := a+d^3 ; end if; end do; a ; end proc:
seq(A050451(n), n=1..100) ; # R. J. Mathar, Jan 07 2011
MATHEMATICA
a[n_] := DivisorSum[n, #^3 &, Mod[#, 4]==1 &]; Array[a, 30] (* Amiram Eldar, Dec 06 2018 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Dec 23 1999
STATUS
approved