OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..388
FORMULA
a(prime(n)^m) = 1 - prime(n)^prime(n) = -A088730(n) for m > 0.
G.f.: Sum_{k>=1} mu(k)*k^k*x^k/(1 - x^k). - Ilya Gutkovskiy, Apr 20 2019
EXAMPLE
a(6) = 1 - 2^2 - 3^3 + 6^6 = 46626.
MATHEMATICA
Array[DivisorSum[#, MoebiusMu[#]*#^# &] &, 21] (* Michael De Vlieger, Apr 21 2019 *)
PROG
(PARI) {a(n) = sumdiv(n, d, moebius(d)*d^d)}
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, Apr 20 2019
STATUS
approved