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

1, 0, 1, 26, 1089, 70124, 6495985, 821315214, 136115947009, 28651724077976, 7470040450004001, 2363470644596843330, 892244303052345224641, 396227360441775922668036, 204487588996059177697597969, 121370399839482643287189048374

Offset

0,4

Links

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

Formula

a(n) = n! * [x^n] exp(-x/(1 - n*x)) / (1 - n*x).

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

a(n) ~ sqrt(2*Pi) * BesselJ(0,2) * n^(2*n + 1/2) / exp(n). - Vaclav Kotesovec, Dec 18 2019

Mathematica

Join[{1}, Table[n! Sum[(-1)^k Binomial[n, k] n^(n - k)/k!, {k, 0, n}], {n, 1, 15}]]
Join[{1}, Table[n^n n! LaguerreL[n, 1/n], {n, 1, 15}]]
Table[n! SeriesCoefficient[Exp[-x/(1 - n x)]/(1 - n x), {x, 0, n}], {n, 0, 15}]

Prog

(Magma) [Factorial(n)*&+[(-1)^k*Binomial(n, k)*n^(n-k)/Factorial(k):k in [0..n]]:n in [0..15]]; // Marius A. Burtea, Dec 18 2019

Crossrefs

Cf. A009940, A025166, A061711, A277423, A307885, A318224, A319392, A330260.

Sequence in context: A160261 A357376 A241874 * A037138 A091429 A200721

Adjacent sequences: A330494 A330495 A330496 * A330498 A330499 A330500

Keyword

nonn

Author

Ilya Gutkovskiy, Dec 18 2019

Status

approved