OFFSET
0,3
LINKS
Robert Israel, Table of n, a(n) for n = 0..232
FORMULA
a(n) ~ BesselI(1, 2) * sqrt(2*Pi) * n^(2*n-1/2) / exp(n). - Vaclav Kotesovec, Oct 01 2017
a(n) = n! * Sum_{k=1..n} n^(n-k) * binomial(n-1,k-1)/k! for n > 0. - Seiichi Manyama, Feb 03 2021
MAPLE
S:=series(exp(x/(1-n*x)), x, 31):
seq(coeff(S, x, n)*n!, n=0..30); # Robert Israel, Oct 01 2017
MATHEMATICA
Table[n! SeriesCoefficient[Exp[x/(1 - n x)], {x, 0, n}], {n, 0, 16}]
Join[{1}, Table[n! SeriesCoefficient[Product[Exp[n^k x^(k + 1)], {k, 0, n}], {x, 0, n}], {n, 1, 16}]]
Join[{1}, Table[Sum[n^(n - k) n!/k! Binomial[n - 1, k - 1], {k, n}], {n, 1, 16}]]
Join[{1}, Table[n^n (n - 1)! Hypergeometric1F1[1 - n, 2, -1/n], {n, 1, 16}]]
PROG
(PARI) {a(n) = if(n==0, 1, n!*sum(k=1, n, n^(n-k)*binomial(n-1, k-1)/k!))} \\ Seiichi Manyama, Feb 03 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Oct 01 2017
STATUS
approved