OFFSET
1,2
COMMENTS
a(A004613(n))=0.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..10000
FORMULA
G.f.: Sum_{k>=1} x^(2*k)*(2 + 3*x^k + 4*x^(2*k) + 2*x^(4*k) + x^(5*k))/(1 - x^(4*k))^2. - Ilya Gutkovskiy, Sep 12 2019
MATHEMATICA
sd[n_]:= Total[Select[Divisors[n], !IntegerQ[(# - 1) / 4]&]]; Array[sd, 100] (* Vincenzo Librandi, May 17 2013 *)
Table[DivisorSum[n, #&, (!IntegerQ[(#-1)/4]&)], {n, 80}] (* Harvey P. Dale, Nov 30 2019 *)
PROG
(PARI) for(n=1, 100, print1(sumdiv(n, d, if(d%4!=1, d))", "))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Ralf Stephan, Apr 02 2003
STATUS
approved