OFFSET
1,2
COMMENTS
If n is prime or semiprime, then a(n) = 2n+1.
MAPLE
with(numtheory): a:=n->sigma(n)+phi(n)-mobius(n): seq(a(n), n=1..100);
MATHEMATICA
Table[DivisorSigma[1, n] + EulerPhi[n] - MoebiusMu[n], {n, 100}]
PROG
(PARI) a(n) = {my(f = factor(n)); sigma(f) + eulerphi(f) - moebius(f); } \\ Amiram Eldar, Nov 14 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wesley Ivan Hurt, Dec 22 2016
STATUS
approved
