%I #31 Feb 15 2022 08:18:49
%S 1,1,1,2,2,1,4,4,3,1,8,8,7,4,1,16,16,15,11,5,1,32,32,31,26,16,6,1,64,
%T 64,63,57,42,22,7,1,128,128,127,120,99,64,29,8,1,256,256,255,247,219,
%U 163,93,37,9,1
%N Binomial transform of A128174.
%C Row sums = A045623: (1, 2, 5, 12, 28, 64, 144, ...).
%C A128176 = A128174 * A007318.
%C Riordan array ((1-x)/(1-2x),x/(1-x)). - _Paul Barry_, Oct 02 2010
%C Fusion of polynomial sequences p(n,x) = (x+1)^n and q(n,x) = x^n + x^(n-1) + ... + x + 1; see A193722 for the definition of fusion. - _Clark Kimberling_, Aug 04 2011
%F A007318 * A128174 as infinite lower triangular matrices.
%F Antidiagonals of an array in which the first row = (1, 1, 2, 4, 8, 16, ...); and (n+1)-th row = partial sums of n-th row.
%F exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(4 + 4*x + 3*x^2/2! + x^3/3!) = 4 + 8*x + 15*x^2/2! + 26*x^3/3! + 42*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - _Peter Bala_, Dec 22 2014
%e First few rows of the triangle:
%e 1;
%e 1, 1;
%e 2, 2, 1;
%e 4, 4, 3, 1;
%e 8, 8, 7, 4, 1;
%e 16, 16, 15, 11, 5, 1;
%e 32, 32, 31, 26, 16, 6, 1;
%e 64, 64, 63, 57, 42, 22, 7, 1;
%e ...
%e From _Paul Barry_, Oct 02 2010: (Start)
%e Production matrix is
%e 1, 1;
%e 1, 1, 1;
%e 0, 0, 1, 1;
%e 0, 0, 0, 1, 1;
%e 0, 0, 0, 0, 1, 1;
%e 0, 0, 0, 0, 0, 1, 1;
%e 0, 0, 0, 0, 0, 0, 1, 1;
%e 0, 0, 0, 0, 0, 0, 0, 1, 1;
%e 0, 0, 0, 0, 0, 0, 0, 0, 1, 1;
%e ...
%e Matrix logarithm is
%e 0;
%e 1, 0;
%e 1, 2, 0;
%e 1, 1, 3, 0;
%e 1, 1, 1, 4, 0;
%e 1, 1, 1, 1, 5, 0;
%e 1, 1, 1, 1, 1, 6, 0;
%e 1, 1, 1, 1, 1, 1, 7, 0;
%e 1, 1, 1, 1, 1, 1, 1, 8, 0;
%e 1, 1, 1, 1, 1, 1, 1, 1, 9, 0;
%e 1, 1, 1, 1, 1, 1, 1, 1, 1, 10, 0;
%e ... (End)
%e .
%e First few rows of the array:
%e 1, 1, 2, 4, 8, 16, ...
%e 1, 2, 4, 8, 16, 32, ...
%e 1, 3, 7, 15, 31, 63, ...
%e 1, 4, 11, 26, 57, 120, ...
%e 1, 5, 16, 42, 99, 219, ...
%e ... (End)
%p A193820 := (n,k) -> `if`(k=0 or n=0, 1, A193820(n-1,k-1)+A193820(n-1,k));
%p A128175 := (n,k) -> A193820(n-1,n-k);
%p seq(print(seq(A128175(n,k),k=0..n)),n=0..10); # _Peter Luschny_, Jan 22 2012
%t z = 10; a = 1; b = 1;
%t p[n_, x_] := (a*x + b)^n
%t q[0, x_] := 1
%t q[n_, x_] := x*q[n - 1, x] + 1; q[n_, 0] := q[n, x] /. x -> 0;
%t t[n_, k_] := Coefficient[p[n, x], x^k]; t[n_, 0] := p[n, x] /. x -> 0;
%t w[n_, x_] := Sum[t[n, k]*q[n + 1 - k, x], {k, 0, n}]; w[-1, x_] := 1
%t g[n_] := CoefficientList[w[n, x], {x}]
%t TableForm[Table[Reverse[g[n]], {n, -1, z}]]
%t Flatten[Table[Reverse[g[n]], {n, -1, z}]] (* A193820 *)
%t TableForm[Table[g[n], {n, -1, z}]]
%t Flatten[Table[g[n], {n, -1, z}]] (* A128175 *)
%t (* _Clark Kimberling_, Aug 06 2011 *)
%t (* function dotTriangle[] is defined in A128176 *)
%t a128175[r_] := dotTriangle[Binomial, If[EvenQ[#1 + #2], 1, 0]&, r]
%t TableForm[a128174[7]] (* triangle *)
%t Flatten[a128174[9]] (* data *) (* _Hartmut F. W. Hoft_, Mar 15 2017 *)
%Y Cf. A045623, A128176, A007318.
%K nonn,tabl
%O 1,4
%A _Gary W. Adamson_, Feb 17 2007