login
A347251
a(n) = Sum_{d|n} mu(d)*mu(n/d)*d^n.
2
1, -5, -28, 16, -3126, 47450, -823544, 0, 19683, 10009766650, -285311670612, -2176786432, -302875106592254, 11112685048647250, 437893920912786408, 0, -827240261886336764178, -101560344088905, -1978419655660313589123980, -100000000000001048576
OFFSET
1,2
LINKS
FORMULA
If p is prime, a(p) = -1 - p^p.
MATHEMATICA
a[n_] := DivisorSum[n, MoebiusMu[#] * MoebiusMu[n/#] * #^n &]; Array[a, 20] (* Amiram Eldar, Aug 24 2021 *)
PROG
(PARI) a(n) = sumdiv(n, d, moebius(d)*moebius(n/d)*d^n);
CROSSREFS
Diagonal of A347227.
Sequence in context: A052470 A351749 A016087 * A321236 A351774 A344579
KEYWORD
sign
AUTHOR
Seiichi Manyama, Aug 24 2021
STATUS
approved