login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A355843 E.g.f. satisfies log(A(x)) = x * (exp(x) - 1) * A(x). 10
1, 0, 2, 3, 40, 185, 2556, 22057, 349616, 4519377, 83642860, 1439639201, 31015493928, 663158322697, 16468280168900, 418772642545545, 11847925722273376, 348085509493265825, 11091199095506163420, 368912674236287743633, 13099432280183074041560 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
LINKS
FORMULA
E.g.f.: exp( -LambertW(x * (1 - exp(x))) ).
E.g.f.: LambertW(x * (1 - exp(x))) / (x * (1 - exp(x))).
a(n) ~ sqrt(1 + exp(1+r)*r^2) * n^(n-1) / (exp(n-1) * r^n), where r = 0.528399250336668412340528181936966763473482889289226687323... is the root of the equation exp(1+r) - exp(1) = 1/r. - Vaclav Kotesovec, Jul 21 2022
a(n) = n! * Sum_{k=0..floor(n/2)} (k+1)^(k-1) * Stirling2(n-k,k)/(n-k)!. - Seiichi Manyama, Aug 28 2022
MATHEMATICA
nmax = 20; A[_] = 1;
Do[A[x_] = Exp[x*(Exp[x] - 1)*A[x]] + O[x]^(nmax+1) // Normal, {nmax}];
CoefficientList[A[x], x]*Range[0, nmax]! (* Jean-François Alcover, Mar 04 2024 *)
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(exp(-lambertw(x*(1-exp(x))))))
(PARI) a(n) = n!*sum(k=0, n\2, (k+1)^(k-1)*stirling(n-k, k, 2)/(n-k)!); \\ Seiichi Manyama, Aug 28 2022
CROSSREFS
Sequence in context: A323734 A347817 A088984 * A097170 A231797 A355235
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jul 18 2022
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 18:04 EDT 2024. Contains 371254 sequences. (Running on oeis4.)