|
|
A371038
|
|
E.g.f. satisfies A(x) = exp(x^2*A(x)) / (1-x).
|
|
1
|
|
|
1, 1, 4, 18, 132, 1140, 12720, 164640, 2514960, 43500240, 850076640, 18418609440, 439831909440, 11457415569600, 323707663319040, 9854548934630400, 321709145793235200, 11209975693710393600, 415330670608805952000, 16303720885477254028800
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
E.g.f.: LambertW( -x^2/(1-x) ) / (-x^2).
a(n) = n! * Sum_{k=0..floor(n/2)} (k+1)^(k-1) * binomial(n-k,n-2*k)/k!.
a(n) ~ exp(2) * sqrt(1 + 4*exp(1) - sqrt(1 + 4*exp(1))) * 2^(n + 3/2) * n^(n-1) / ((1 + 2*exp(1) - sqrt(1 + 4*exp(1)))*(-1 + sqrt(1 + 4*exp(1)))^(n+1)). - Vaclav Kotesovec, Mar 12 2024
|
|
PROG
|
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(lambertw(-x^2/(1-x))/(-x^2)))
(PARI) a(n) = n!*sum(k=0, n\2, (k+1)^(k-1)*binomial(n-k, n-2*k)/k!);
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|