OFFSET
1,3
MATHEMATICA
d[1] = 0; d[n_] := n * Plus @@ ((Last[#]/First[#]) & /@ FactorInteger[n]); a[n_] := DivisorSum[n, # * (-1)^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, d*(-1)^ad(d)); \\ Michel Marcus, May 12 2022
CROSSREFS
KEYWORD
sign
AUTHOR
Wesley Ivan Hurt, May 12 2022
STATUS
approved