|
%I #18 Dec 31 2023 10:34:38
%S 1,1,-4,11,21,-674,6551,-33479,-174114,7478121,-117699599,1090997976,
%T 865365421,-302755297739,7922094623596,-127940743443649,
%U 974028543402401,21377262410290446,-1179125036786673989,31760741865879345821,-552216474702144564074,2814873629049018241701
%N Expansion of e.g.f. exp((1 - exp(-5*x))/5).
%H Robert Israel, <a href="/A318180/b318180.txt">Table of n, a(n) for n = 0..459</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/BellPolynomial.html">Bell Polynomial</a>
%F a(n) = Sum_{k=0..n} (-5)^(n-k)*Stirling2(n,k).
%F a(0) = 1; a(n) = Sum_{k=1..n} (-5)^(k-1)*binomial(n-1,k-1)*a(n-k).
%F a(n) = (-5)^n*BellPolynomial_n(-1/5). - _Peter Luschny_, Aug 20 2018
%p seq((-5)^n*BellB(n,-1/5),n=0..30); # _Robert Israel_, Nov 11 2020
%t nmax = 21; CoefficientList[Series[Exp[(1 - Exp[-5 x])/5], {x, 0, nmax}], x] Range[0, nmax]!
%t Table[Sum[(-5)^(n - k) StirlingS2[n, k], {k, 0, n}], {n, 0, 21}]
%t a[n_] := a[n] = Sum[(-5)^(k - 1) Binomial[n - 1, k - 1] a[n - k], {k, 1, n}]; a[0] = 1; Table[a[n], {n, 0, 21}]
%t Table[(-5)^n BellB[n, -1/5], {n, 0, 21}] (* _Peter Luschny_, Aug 20 2018 *)
%o (PARI) my(x = 'x + O('x^25)); Vec(serlaplace(exp((1 - exp(-5*x))/5))) \\ _Michel Marcus_, Nov 11 2020
%Y Cf. A005011, A008548, A009235, A014182, A317996, A318179, A318181.
%K sign
%O 0,3
%A _Ilya Gutkovskiy_, Aug 20 2018
|
