OFFSET
1,6
COMMENTS
From Robert G. Wilson v, Mar 26 2015: (Start)
a(n) = 0 for n = 1, 2, 4, 5, 8, 10, 13, 16, 17, 20, 25, ... (A072437).
a(n) = 1 for n = 3, 7, 11, 19, 23, 31, 43, 47, 59, 67, ... (A002145). (End)
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Robert G. Wilson v)
FORMULA
MATHEMATICA
a[n_] := Total[(n/Select[Divisors@ n, Mod[#, 4] == 3 &])^3]; Array[a, 64] (* Robert G. Wilson v, Mar 26 2015 *)
a[n_] := DivisorSum[n, #^3 &, Mod[n/#, 4] == 3 &]; Array[a, 50] (* Amiram Eldar, Nov 05 2023 *)
PROG
(PARI) a(n) = sumdiv(n, d, ((n/d % 4)== 3)* d^3); \\ Michel Marcus, Mar 26 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Dec 23 1999
EXTENSIONS
Offset changed from 0 to 1 by Robert G. Wilson v, Mar 27 2015
STATUS
approved