OFFSET
1,1
LINKS
FORMULA
MATHEMATICA
f[p_, e_] := NextPrime[p, -1]^e; s[1] = 1; s[n_] := If[EvenQ[n], n/2, Times @@ f @@@ FactorInteger[n]]; sinv[1] = 1; sinv[n_] := sinv[n] = -DivisorSum[n, sinv[#] * s[n/#] &, # < n &]; a[n_] := DivisorSum[n, # * MoebiusMu[#] * s[n/#] + # * sinv[n/#] &]; Array[a, 100] (* Amiram Eldar, Nov 18 2021 *)
CROSSREFS
KEYWORD
sign
AUTHOR
Antti Karttunen, Nov 18 2021
STATUS
approved