%I #16 Mar 28 2020 06:38:09
%S 2,3,5,4,11,9,5,19,26,14,6,29,55,50,20,7,41,99,125,85,27,8,55,161,259,
%T 245,133,35,9,71,244,476,574,434,196,44,10,89,351,804,1176,1134,714,
%U 276,54,11,109,485,1275,2190,2562,2058,1110,375,65,12,131,649,1925
%N Triangular array: the fission of (p(n,x)) by (q(n,x)), where p(n,x)=x*p(n-1,x)+n+1 with p(0,x)=1, and q(n,x)=(x+1)^n.
%C See A193842 for the definition of fission of two sequences of polynomials or triangular arrays.
%e First six rows:
%e 2
%e 3...5
%e 4...11....9
%e 5...19...26...14
%e 6...29...55...50...20
%e 7...41...99...125..85...27
%p # The function 'fission' is defined in A193842.
%p p := (n,x) -> `if`(n=0,1,x*p(n-1,x)+n+1);
%p q := (n,x) -> (x+1)^n;
%p A193971_row := n -> fission(p, q, n);
%p for n from 0 to 5 do A193971_row(n) od; # _Peter Luschny_, Jul 23 2014
%t z = 11;
%t p[0, x_] := 1; p[n_, x_] := x*p[n - 1, x] + n + 1;
%t q[n_, x_] := (x + 1)^n
%t p1[n_, k_] := Coefficient[p[n, x], x^k];
%t p1[n_, 0] := p[n, x] /. x -> 0;
%t d[n_, x_] := Sum[p1[n, k]*q[n - 1 - k, x], {k, 0, n - 1}]
%t h[n_] := CoefficientList[d[n, x], {x}]
%t TableForm[Table[Reverse[h[n]], {n, 0, z}]]
%t Flatten[Table[Reverse[h[n]], {n, -1, z}]] (* A193971 *)
%t TableForm[Table[h[n], {n, 0, z}]]
%t Flatten[Table[h[n], {n, -1, z}]] (* A193972 *)
%o (Sage) # uses[fission from A193842]
%o p = lambda n,x: x*p(n-1,x)+n+1 if n > 0 else 1
%o q = lambda n,x: (x+1)^n
%o A193971_row = lambda n: fission(p, q, n);
%o for n in range(7): A193971_row(n) # _Peter Luschny_, Jul 23 2014
%Y Cf. A193842, A193972.
%K nonn,tabl
%O 0,1
%A _Clark Kimberling_, Aug 10 2011
|