|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
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
|
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|