%I #34 Jun 03 2024 22:50:23
%S 1,1,1,2,1,1,5,7,1,1,15,87,34,1,1,52,1657,2971,209,1,1,203,43833,
%T 513559,163121,1546,1,1,877,1515903,149670844,326922081,12962661,
%U 13327,1,1,4140,65766991,66653198353,1346634725665,363303011071,1395857215,130922,1,1
%N Triangle of certain generalized Bell numbers.
%C Let B_{n}(x) = sum_{j>=0}(exp(j!/(j-n)!*x-1)/j!) and
%C S(n,k) = k! [x^k] taylor(B_{n}(x)), where [x^k] denotes the
%C coefficient of x^k in the Taylor series for B_{n}(x).
%C Then S(n,k) (n>0, k>=0) is the square array representation of the triangle.
%C To illustrate the cross-references of T(n,k) when written as a square array.
%C n\k A000012,A000012,A002720,A069948,A182925,A182924,...
%C 0: A000012: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
%C 1: A000110: 1, 1, 2, 5, 15, 52, 203, 877, 4140, ...
%C 2: A020556: 1, 1, 7, 87, 1657, 43833, 1515903, ...
%C 3: A069223: 1, 1, 34, 2971, 513559, 149670844, ...
%C 4: A071379: 1, 1, 209, 163121, 326922081, ...
%C 5: A090209: 1, 1, 1546, 12962661, 363303011071,...
%C 6: ... 1, 1, 13327, 1395857215, 637056434385865,...
%C Note that the sequence T(0,k) is not included in the data.
%C - _Peter Luschny_, Mar 27 2011
%D P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, Phys. Lett. A 309 (2003) 198-205.
%D M. Schork, On the combinatorics of normal ordering bosonic operators and deforming it, J. Phys. A 36 (2003) 4651-4665.
%H P. Blasiak, K. A. Penson and A. I. Solomon, <a href="http://www.arXiv.org/abs/quant-ph/0402027">The general boson normal ordering problem.</a>
%H W. Lang, <a href="/A090210/a090210.txt">First 8 rows</a>.
%H K. A. Penson, P. Blasiak, A. Horzela, A. I. Solomon and G. H. E. Duchamp,<a href="http://arxiv.org/abs/0904.0369">Laguerre-type derivatives: Dobinski relations and combinatorial identities</a>, J. Math. Phys. 50, 083512 (2009).
%F a(n, m) = Bell(m;n-(m-1)), n>= m-1 >=0, with Bell(m;k) := Sum_{p=m..m*k} S2(m;k, p), where S2(m;k, p) := (((-1)^p)/p!) * Sum_{r=m..p} ((-1)^r)*binomial(p, r)*fallfac(r, m)^k; with fallfac(n, m) := A008279(n, m) (falling factorials) and m<=p<=k*m, k>=1, m=1, 2, ..., else 0. From eqs.(6) with r=s->m and eq.(19) with S_{r, r}(n, k)-> S2(r;n, k) of the Blasiak et al. reference. [Corrected by _Sean A. Irvine_, Jun 03 2024]
%F a(n, m) = (Sum_{k>=m} fallfac(k, m)^(n-(m-1)))/exp(1), n>=m-1>=0, else 0. From eq.(26) with r->m of the Schork reference which is rewritten eq.(11) of the original Blasiak et al. reference.
%F E.g.f. m-th column (no leading zeros): (Sum_{k>=m} exp(fallfac(k, m)*x)/k!) + A000522(m)/m!)/exp(1). Rewritten from the top of p. 4656 of the Schork reference.
%e Triangle begins:
%e 1;
%e 1, 1;
%e 2, 1, 1;
%e 5, 7, 1, 1;
%e 15, 87, 34, 1, 1;
%e 52, 1657, 2971, 209, 1, 1;
%e 203, 43833, 513559, 163121, 1546, 1, 1;
%p A090210_AsSquareArray := proc(n,k) local r,s,i;
%p if k=0 then 1 else r := [seq(n+1,i=1..k-1)]; s := [seq(1,i=1..k-1)];
%p exp(-x)*n!^(k-1)*hypergeom(r,s,x); round(evalf(subs(x=1,%),99)) fi end:
%p seq(lprint(seq(A090210_AsSquareArray(n,k),k=0..6)),n=0..6);
%p # _Peter Luschny_, Mar 30 2011
%t t[n_, k_] := t[n, k] = Sum[(n+j)!^(k-1)/(j!^k*E), {j, 0, Infinity}]; t[_, 0] = 1;
%t Flatten[ Table[ t[n-k+1, k], {n, 0, 8}, {k, n, 0, -1}]][[1 ;; 43]] (* _Jean-François Alcover_, Jun 17 2011 *)
%Y Cf. A000110, A020556, A069223, A071379, A090209, A002720, A069948, A182924, A182925, A182933.
%Y T(n,n) gives A070227.
%K nonn,easy,tabl
%O 1,4
%A _Wolfdieter Lang_, Dec 01 2003