|
|
A362281
|
|
a(n) = n! * Sum_{k=0..floor(n/2)} n^k * binomial(n-k,k)/(n-k)!.
|
|
2
|
|
|
1, 1, 5, 19, 241, 1601, 32581, 308995, 8655809, 106673761, 3805452901, 57704760851, 2500580809585, 45018720191329, 2295683481085541, 47848514992963651, 2806491306922172161, 66464103165835330625, 4407449313521981148229, 116893033842508769526931
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
a(n) = n! * [x^n] exp(x + n*x^2).
E.g.f.: exp( sqrt( -LambertW(-2*x^2)/2 ) ) / (1 + LambertW(-2*x^2)).
a(n) ~ (1 + (-1)^n/exp(sqrt(2))) * 2^((n-1)/2) * n^n / exp(n/2 - 1/sqrt(2)). - Vaclav Kotesovec, Apr 15 2023
|
|
PROG
|
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(exp(sqrt(-lambertw(-2*x^2)/2))/(1+lambertw(-2*x^2))))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|