|
%I #34 Dec 23 2021 06:02:57
%S 1,-1,2,-3,-20,370,-4074,34293,-138312,-2932533,106271090,-2192834490,
%T 32208497124,-206343936097,-7657279887698,412496622532785,
%U -12455477719752976,260294034150380430,-2256541295745391542,-122593550603339550843,8728842979656718306780
%N a(n) = n! * [x^n] exp(n*(1 - exp(x))).
%H Alois P. Heinz, <a href="/A292866/b292866.txt">Table of n, a(n) for n = 0..415</a>
%F a(n) = exp(n) * Sum_{k>=0} (-n)^k*k^n/k!. - _Ilya Gutkovskiy_, Jul 13 2019
%F a(n) = Sum_{k=0..n} (-n)^k * Stirling2(n,k). - _Seiichi Manyama_, Jul 28 2019
%F a(n) = BellPolynomial(n, -n). - _Peter Luschny_, Dec 23 2021
%p b:= proc(n, k) option remember; `if`(n=0, 1,
%p -(1+add(binomial(n-1, j-1)*b(n-j, k), j=1..n-1))*k)
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..30); # _Alois P. Heinz_, Sep 25 2017
%t Table[n!*SeriesCoefficient[E^(n*(1 - E^x)),{x,0,n}], {n,0,20}] (* _Vaclav Kotesovec_, Sep 25 2017 *)
%t a[n_] := BellB[n, -n]; Table[a[n], {n, 0, 20}] (* _Peter Luschny_, Dec 23 2021 *)
%o (Ruby)
%o def ncr(n, r)
%o return 1 if r == 0
%o (n - r + 1..n).inject(:*) / (1..r).inject(:*)
%o end
%o def A(k, n)
%o ary = [1]
%o (1..n).each{|i| ary << k * (0..i - 1).inject(0){|s, j| s + ncr(i - 1, j) * ary[j]}}
%o ary
%o end
%o def A292866(n)
%o (0..n).map{|i| A(-i, i)[-1]}
%o end
%o p A292866(20)
%o (PARI) {a(n) = sum(k=0, n, (-n)^k*stirling(n, k, 2))} \\ _Seiichi Manyama_, Jul 28 2019
%Y Main diagonal of A292861.
%Y Cf. A242817.
%K sign
%O 0,3
%A _Seiichi Manyama_, Sep 25 2017
|
