A368725 a(n) = n! * Sum_{k=0..n} (-1)^(n-k) * k^n / k!.

1, 1, 2, 9, 100, 1065, 10626, 224161, 4598504, 46288017, 2509940710, 84061763841, -1602021820596, 164372205860473, 5216105126641514, -883395389739028095, 79008645559978113616, -1023235751229436800735, -651030746777115881959602, 113943411938145511923004513

Offset

0,3

Links

Table of n, a(n) for n=0..19.

Eric Weisstein's World of Mathematics, Bell Polynomial.

Wikipedia, Touchard polynomials.

Formula

E.g.f.: Sum_{k>=0} (k * x)^k / (k! * (1 + k * x)).

a(n) = n! * [x^n] B_n(x) * exp(x) / (1+x), where B_n(x) = Bell polynomials.

a(n) ~ A000587(n) * (-1)^n * exp(-1) * n!. - Vaclav Kotesovec, Jul 18 2025

Mathematica

Join[{1}, Table[n!*Sum[(-1)^(n-k)*k^n/k!, {k, 0, n}], {n, 1, 20}]] (* Vaclav Kotesovec, Jul 18 2025 *)

Prog

(PARI) a(n) = n!*sum(k=0, n, (-1)^(n-k)*k^n/k!);

Crossrefs

Main diagonal of A368724.

Cf. A000587, A256016.

Sequence in context: A357825 A187647 A322645 * A277180 A013520 A378045

Adjacent sequences: A368722 A368723 A368724 * A368726 A368727 A368728

Keyword

sign

Author

Seiichi Manyama, Jan 04 2024

Status

approved