A376576 E.g.f. A(x) satisfies A(x) = exp(x*A(x)/(1 - x^3)).
1, 1, 3, 16, 149, 1656, 22567, 372184, 7141689, 156630448, 3871782251, 106504501104, 3227742350197, 106879926110296, 3839600650843791, 148746681984864856, 6181806007303273073, 274355581868776940256, 12951023558423725477459, 647956009961120527442272
Offset: 0
Keywords
Links
- Eric Weisstein's World of Mathematics, Lambert W-Function.
Programs
-
PARI
a(n) = n!*sum(k=0, n\3, (n-3*k+1)^(n-3*k-1)*binomial(n-2*k-1, k)/(n-3*k)!);
-
PARI
my(N=20, x='x+O('x^N)); Vec(serlaplace(exp(-lambertw(-x/(1-x^3)))))
Formula
a(n) = n! * Sum_{k=0..floor(n/3)} (n-3*k+1)^(n-3*k-1) * binomial(n-2*k-1,k)/(n-3*k)!.
E.g.f.: exp( -LambertW(-x/(1-x^3)) ).
From Vaclav Kotesovec, Oct 10 2024: (Start)
E.g.f.: -LambertW(-x/(1-x^3))*(1-x^3)/x.
a(n) ~ sqrt(2^(2/3) * 3^(5/3) / ((2*(9 + sqrt(81 + 12*exp(3))))^(1/3) - 2*exp(1)*(3/(9 + sqrt(81 + 12*exp(3))))^(1/3)) - 2*exp(1)) * 2^(2*n/3) * 3^(4*n/3) * ((9 + sqrt(81 + 12*exp(3)))^(1/3) / (2^(1/3) * (3*(9 + sqrt(81 + 12*exp(3))))^(2/3) - 6*exp(1)))^n * n^(n-1) / exp(n - 1/2). (End)