A090210 - OEIS

%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