A355376 a(n) = Sum_{k=0..n} k! * (-k)^(n-k) * Stirling2(n,k).

1, 1, 1, -5, -29, 271, 3091, -61025, -744029, 34875871, 211095331, -37415273345, 300267009571, 61080483836191, -2133136977892829, -119576844586022465, 11752559492568148771, 94348367247493654111, -68793303669649907424989, 2764486881197709482575615

Offset

0,4

Links

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

Formula

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

Mathematica

a[n_] := Sum[k! * (-k)^(n - k) * StirlingS2[n, k], {k, 0, n}]; a[0] = 1; Array[a, 20, 0] (* Amiram Eldar, Jun 30 2022 *)

Prog

(PARI) a(n) = sum(k=0, n, k!*(-k)^(n-k)*stirling(n, k, 2));
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(sum(k=0, N, (1-exp(-k*x))^k/k^k)))

Crossrefs

Cf. A229234, A355373, A355375.

Sequence in context: A181356 A177440 A292567 * A332517 A332469 A112799

Adjacent sequences: A355373 A355374 A355375 * A355377 A355378 A355379

Keyword

sign

Author

Seiichi Manyama, Jun 30 2022

Status

approved