A367493 a(n) = Sum_{d|n} (n-d)^n.

0, 1, 8, 97, 1024, 20450, 279936, 7509953, 144295424, 4570291850, 100000000000, 4491754172274, 106993205379072, 5221973073321002, 171975117132398592, 8931527427394008833, 295147905179352825856, 20290116242888952838355, 708235345355337676357632, 51879761166564630630389778

Offset

1,3

Links

Table of n, a(n) for n=1..20.

Formula

a(n) = Sum_{k=0..n} (-1)^k*binomial(n,k)*n^(n-k)*sigma_k(n).

Mathematica

a[n_]:=Sum[(n-d)^n, {d, Divisors[n]}]; Array[a, 20] (* Stefano Spezia, Nov 20 2023 *)

Prog

(Python)
from sympy import divisors
def A367493(n): return sum((n-d)**n for d in divisors(n, generator=True))

Crossrefs

Cf. A094471, A367326, A367327, A367368.

Sequence in context: A320398 A218451 A262777 * A365845 A292747 A083182

Adjacent sequences: A367490 A367491 A367492 * A367494 A367495 A367496

Keyword

nonn

Author

Chai Wah Wu, Nov 20 2023

Status

approved