%I #21 Jun 28 2018 08:24:19
%S 4,12,16,24,144,64,24,816,1152,256,0,3360,12480,7680,1024,0,10080,
%T 100800,134400,46080,4096,0,20160,645120,1747200,1182720,258048,16384,
%U 0,20160,3306240,18305280,22364160,9117696,1376256,65536,0,0,13305600,159667200,341591040,235339776,63995904,7077888,262144,0,0,39916800
%N Exponential Riordan array (1,4*x+6*x^2+4*x^3+x^4).
%C The column of index 0 contains a 1 followed by zeros and is not incorporated here.
%C Also the Bell polynomials of the second kind B(n,k)(4,12,24,24).
%C If the argument vector is generalized to contain falling powers of a variable x, B(n,k)(4*x^3,12*x,24*x,24) =sum_{j=0..k} binomial(k,j) *sum_{i=j..n-k+j} 6^(i-j) *binomial(j,i-j) *binomial(k-j,n-3*k+3*j-i) *4^(4*k-n-2*j) *x^(4*k-n) *n!/k!.
%H M. Abbas, S. Bouroubi, <a href="http://dx.doi.org/10.1016/j.disc.2004.08.023">On new identities for Bell's polynomials</a>, Disc. Math 293. (1-3) (2005) 5-10
%H Vladimir Kruchinin, <a href="http://arxiv.org/abs/1104.5065">Derivation of Bell Polynomials of the Second Kind</a>, arXiv:1104.5065
%H John Riordan, <a href="http://dx.doi.org/10.1090/S0002-9904-1946-08621-8">Derivatives of composite functions</a>, Bull. Am. Math. Soc. 52 (1946) 664
%H Eric W. Weisstein, <a href="http://mathworld.wolfram.com/BellPolynomial.html">Bell Polynomial</a>, MathWorld
%F B(n,k) = (n!/k!)*sum_{j=0..k} binomial(k,j) *sum_{i=j..n-k+j} 6^(i-j) *binomial(j,i-j) *binomial(k-j,n-3*k+3*j-i)) *4^(4*k-n-2*j).
%e Triangle begins:
%e 4;
%e 12,16;
%e 24,144,64;
%e 24,816,1152,256;
%e 0,3360,12480,7680,1024;
%e 0,10080,100800,134400,46080,4096;
%e 0,20160,645120,1747200,1182720,258048,16384;
%e 0,20160,3306240,18305280,22364160,9117696,1376256,65536;
%p # The function BellMatrix is defined in A264428.
%p # Adds (1,0,0,0, ..) as column 0.
%p BellMatrix(n -> `if`(n<4,[4,12,24,24][n+1],0), 9); # _Peter Luschny_, Jan 29 2016
%t BellMatrix[f_, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len - 1}, {k, 0, len - 1}]];
%t B = BellMatrix[Function[n, If[n<4, {4, 12, 24, 24}[[n+1]], 0]], rows = 12];
%t Table[B[[n, k]], {n, 2, rows}, {k, 2, n}] // Flatten (* _Jean-François Alcover_, Jun 28 2018, after _Peter Luschny_ *)
%o (Maxima)
%o B(n,k):=n!/k!*sum(binomial(k,j)*sum(6^(i-j)*binomial(j,i-j)*binomial(k-j,n-3*k+3*j-i),i,j,n-k+j)*4^(4*k-n-2*j),j,0,k);
%Y Cf. A008279 (its row 4), A188138, A188062, A188066.
%K nonn,tabl
%O 1,1
%A _Vladimir Kruchinin_, Mar 04 2011
|