%I #47 Oct 03 2023 13:14:46
%S 1,4,1,36,12,1,504,192,24,1,9576,3960,600,40,1,229824,100656,17160,
%T 1440,60,1,6664896,3048192,563976,54600,2940,84,1,226606464,107255232,
%U 21095424,2256576,142800,5376,112,1,8837652096,4302305280,887785920,102332160,7254576,325584,9072,144,1
%N Triangle read by rows, the inverse Bell transform of n!*binomial(4,n) (without column 0).
%C Previous name was: Triangle of numbers related to triangle A049223; generalization of Stirling numbers of second kind A008277, Bessel triangle A001497.
%C T(n, m) = S2p(-4; n,m), a member of a sequence of triangles including S2p(-1; n,m) = A001497(n-1,m-1) (Bessel triangle) and ((-1)^(n-m))*S2p(1; n,m) = A008277(n,m) (Stirling 2nd kind). T(n, 1) = A008546(n-1).
%C For the definition of the Bell transform see A264428 and the link. - _Peter Luschny_, Jan 16 2016
%H G. C. Greubel, <a href="/A011801/b011801.txt">Rows n = 1..50 of the triangle, flattened</a>
%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>, arXiv:quant-ph/0402027, 2004.
%H Wolfdieter Lang, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL3/LANG/lang.html">On generalizations of Stirling number triangles</a>, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
%H Peter Luschny, <a href="https://oeis.org/wiki/User:Peter_Luschny/BellTransform">The Bell transform</a>
%H <a href="/index/Be#Bessel">Index entries for sequences related to Bessel functions or polynomials</a>
%F T(n, m) = n!*A049223(n, m)/(m!*5^(n-m)).
%F T(n+1, m) = (5*n-m)*T(n, m) + T(n, m-1), for n >= m >= 1, with T(n, m) = 0, for n < m, and T(n, 0) = 0, T(1, 1) = 1.
%F E.g.f. of n-th column: (1/n!)*( 1 - (1-5*x)^(1/5) )^n.
%F Sum_{k=1..n} T(n, k) = A028575(n).
%e Triangle starts:
%e 1;
%e 4, 1;
%e 36, 12, 1;
%e 504, 192, 24, 1;
%e 9576, 3960, 600, 40, 1;
%e 229824, 100656, 17160, 1440, 60, 1;
%e 6664896, 3048192, 563976, 54600, 2940, 84, 1;
%e 226606464, 107255232, 21095424, 2256576, 142800, 5376, 112, 1;
%t (* First program *)
%t T[n_, m_] /; n>=m>=1:= T[n, m]= (5*(n-1)-m)*T[n-1, m] + T[n-1, m-1]; T[n_, m_] /; n<m=0; T[_, 0]=0; T[1,1]=1;
%t Table[T[n,m], {n,10}, {m,n}]//Flatten (* _Jean-François Alcover_, Jun 20 2018 *)
%t (* Second program *)
%t rows = 10;
%t b[n_, m_]:= BellY[n, m, Table[k! Binomial[4, k], {k, 0, rows}]];
%t T= Table[b[n, m], {n,rows}, {m,rows}]//Inverse//Abs;
%t A011801= Table[T[[n, m]], {n,rows}, {m,n}]//Flatten (* _Jean-François Alcover_, Jun 22 2018 *)
%o (Sage) # uses[inverse_bell_matrix from A264428]
%o # Adds 1,0,0,0, ... as column 0 at the left side of the triangle.
%o inverse_bell_matrix(lambda n: factorial(n)*binomial(4, n), 8) # _Peter Luschny_, Jan 16 2016
%o (Magma)
%o function T(n,k) // T = A011801
%o if k eq 0 then return 0;
%o elif k eq n then return 1;
%o else return (5*(n-1)-k)*T(n-1,k) + T(n-1,k-1);
%o end if;
%o end function;
%o [T(n,k): k in [1..n], n in [1..12]]; // _G. C. Greubel_, Oct 03 2023
%Y Cf. A008277, A008546, A049223, A264428.
%Y Cf. A028575 (row sums).
%Y Triangles with the recurrence T(n,k) = (m*(n-1)-k)*T(n-1,k) + T(n-1,k-1): A010054 (m=1), A001497 (m=2), A004747 (m=3), A000369 (m=4), this sequence (m=5), A013988 (m=6).
%K easy,nonn,tabl
%O 1,2
%A _Wolfdieter Lang_
%E New name from _Peter Luschny_, Jan 16 2016