OFFSET
1,2
COMMENTS
If p is prime, a(p) = p * Sum_{d|p} d^mu(d) = p * (1^1 + p^(-1)) = p * (1 + 1/p) = p + 1.
EXAMPLE
a(6) = 6 * Sum_{d|6} d^mu(d) = 6 * (1^1 + 2^(-1) + 3^(-1) + 6^1) = 47.
MATHEMATICA
Table[n*Sum[i^MoebiusMu[i] (1 - Ceiling[n/i] + Floor[n/i]), {i, n}], {n, 80}]
PROG
(PARI) a(n) = n*sumdiv(n, d, d^moebius(d)); \\ Michel Marcus, May 25 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, May 24 2021
STATUS
approved