|
%I #26 Jun 24 2023 22:06:37
%S 1,1,2,8,41,252,1782,14121,123244,1169832,11960978,130742196,
%T 1518514076,18645970943,241030821566,3268214127548,46338504902485,
%U 685145875623056,10538790233183702,168282662416550040,2784205185437851772,47646587512911994120
%N a(n) = Sum_{k=1..n} C(n-1,k-1) * S2(n,k) for n>0, a(0)=1, where S2(n,k) = A048993(n,k) are Stirling numbers of the 2nd kind.
%H Alois P. Heinz, <a href="/A134055/b134055.txt">Table of n, a(n) for n = 0..518</a>
%F O.g.f.: Sum_{n>=0} (n*x)^n/(1-n*x)^n * exp(-n*x/(1-n*x)) / n!. - _Paul D. Hanna_, Nov 04 2012
%F From _Alois P. Heinz_, Jun 24 2023: (Start)
%F a(n) mod 2 = A037011(n) for n >= 1.
%F a(n) mod 2 = 1 <=> n in { A048297 } or n = 0. (End)
%e O.g.f.: A(x) = 1 + x + 2*x^2 + 8*x^3 + 41*x^4 + 252*x^5 + 1782*x^6 + 14121*x^7 +...
%e where
%e A(x) = 1 + x/(1-x)*exp(-x/(1-x)) + 2^2*x^2/(1-2*x)^2*exp(-2*x/(1-2*x))/2! + 3^3*x^3/(1-3*x)^3*exp(-3*x/(1-3*x))/3! + 4^4*x^4/(1-4*x)^4*exp(-4*x/(1-4*x))/4! +...
%e simplifies to a power series in x with integer coefficients.
%e Illustrate the definition of the terms by:
%e a(4) = 1*1 + 3*7 + 3*6 + 1*1 = 41;
%e a(5) = 1*1 + 4*15 + 6*25 + 4*10 + 1*1 = 252;
%e a(6) = 1*1 + 5*31 + 10*90 + 10*65 + 5*15 + 1*1 = 1782.
%p a:= proc(n) option remember; local b; b:=
%p proc(h, m) option remember; `if`(h=0,
%p binomial(n-1, m-1), m*b(h-1, m)+b(h-1, m+1) )
%p end; b(n, 0)
%p end:
%p seq(a(n), n=0..22); # _Alois P. Heinz_, Jun 24 2023
%t Flatten[{1,Table[Sum[Binomial[n-1,k-1] * StirlingS2[n,k],{k,1,n}],{n,1,20}]}] (* _Vaclav Kotesovec_, Aug 11 2014 *)
%o (PARI) a(n)=if(n==0,1,sum(k=1, n, binomial(n-1, k-1)*polcoeff(1/prod(i=0, k, 1-i*x +x*O(x^(n-k))), n-k)))
%o (PARI) a(n)=polcoeff(sum(k=0,n+1,(k*x)^k/(1-k*x)^k*exp(-k*x/(1-k*x+x*O(x^n)))/k!),n)
%o for(n=0,25,print1(a(n),", ")) \\ _Paul D. Hanna_, Nov 04 2012
%Y Cf. A037011, A048297, A048993, A122455, A218667, A218670, A245059, A245060.
%K nonn
%O 0,3
%A _Paul D. Hanna_, Oct 08 2007
%E An initial '1' was added and definition changed slightly by _Paul D. Hanna_, Nov 04 2012
|
