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

1, 1, 12, 279, 11536, 746525, 69768036, 8902181575, 1487939919936, 315597946293657, 82839437215344100, 26366747854082944451, 10006618140321691249296, 4464690010732922712332149, 2313871692128866349730705924, 1378552938661073773617331110975

Offset

0,3

Links

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

Formula

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

From Vaclav Kotesovec, Feb 19 2022: (Start)

a(n) ~ exp(1) * n! * n^n.

a(n) ~ sqrt(2*Pi) * n^(2*n + 1/2) / exp(n-1). (End)

Mathematica

Join[{1}, Table[n!*Sum[n^(n - k)*(n - k)^k/k!, {k, 0, n}], {n, 1, 20}]] (* Vaclav Kotesovec, Feb 19 2022 *)

Prog

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

Crossrefs

Main diagonal of A351761.

Cf. A332408, A351768, A351779.

Sequence in context: A296383 A296640 A239780 * A295037 A295528 A166337

Adjacent sequences: A351762 A351763 A351764 * A351766 A351767 A351768

Keyword

nonn

Author

Seiichi Manyama, Feb 18 2022

Status

approved