A344430 a(n) = Sum_{k=1..n} mu(k) * k^k.
1, -3, -30, -30, -3155, 43501, -780042, -780042, -780042, 9999219958, -275312450653, -275312450653, -303150419042906, 10808856406515110, 448702746787374485, 448702746787374485, -826791559139549389692, -826791559139549389692
Offset: 1
Keywords
Programs
-
Mathematica
a[n_] := Sum[MoebiusMu[k] * k^k, {k,1,n}]; Array[a, 20] (* Amiram Eldar, May 19 2021 *) Accumulate[Table[MoebiusMu[n]n^n,{n,20}]] (* Harvey P. Dale, Jan 25 2022 *)
-
PARI
a(n) = sum(k=1, n, moebius(k)*k^k);
-
Python
from sympy import mobius def A344430(n): return sum(mobius(k)*k**k for k in range(1,n+1)) # Chai Wah Wu, Apr 05 2023