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