OFFSET
1,2
FORMULA
a(n) = Sum_{d|n} d * ((d') mod 2).
MAPLE
A353956 := proc(n)
add(d*modp(A003415(d), 2), d=numtheory[divisors](n)) ;
end proc:
seq(A353956(n), n=1..80) ; # R. J. Mathar, Jul 05 2022
MATHEMATICA
d[1] = 0; d[n_] := n * Plus @@ ((Last[#]/First[#]) & /@ FactorInteger[n]); a[n_] := DivisorSum[n, # &, OddQ[d[#]] &]; Array[a, 100] (* Amiram Eldar, May 12 2022 *)
PROG
(PARI) ad(n) = vecsum([n/f[1]*f[2]|f<-factor(n+!n)~]); \\ A003415
a(n) = sumdiv(n, d, if (ad(d) %2, d)); \\ Michel Marcus, May 12 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, May 12 2022
STATUS
approved