login
Triangle of coefficients of polynomials v(n,x) jointly generated with A210803; see the Formula section.
3

%I #18 Jan 27 2020 01:32:14

%S 1,2,2,5,8,3,14,27,18,5,41,88,79,40,8,122,284,310,215,80,13,365,912,

%T 1152,980,510,156,21,1094,2917,4144,4091,2660,1150,294,34,3281,9296,

%U 14578,16176,12393,6752,2461,544,55,9842,29526,50436,61638,53730

%N Triangle of coefficients of polynomials v(n,x) jointly generated with A210803; see the Formula section.

%C Row n ends with F(n), where F=A000045 (Fibonacci numbers).

%C Column 1: A007051.

%C Row sums: A000302 (powers of 4).

%C Alternating row sums: 1,0,0,0,0,0,0,0,0,...

%C For a discussion and guide to related arrays, see A208510.

%C Essentially the same triangle as given by (2, 1/2, 3/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (2, -1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - _Philippe Deléham_ Jul 11 2012

%F u(n,x) = u(n-1,x) + x*v(n-1,x) + 1, v(n,x) = (x-1)*u(n-1,x) + (x+3)*v(n-1,x), where u(1,x)=1, v(1,x)=1.

%F T(n,k) = 4*T(n-1,k) + T(n-1,k-1) - 3*T(n-2,k) - 2*T(n-2,k-1) + T(n-2,k-2), T(1,0) = 1, T(2,0) = T(2,1) = 2, T(3,0) = 5, T(3,1) = 8, T(3,2) = 3, T(n,k) = 0 if k < 0 or if k >= n. - _Philippe Deléham_, Jul 11 2012

%F G.f.: (-1+2*x-x*y)*x*y/(-1+4*x+x*y-3*x^2-2*x^2*y+x^2*y^2). - _R. J. Mathar_, Aug 12 2015

%e First five rows:

%e 1;

%e 2, 2;

%e 5, 8, 3;

%e 14, 27, 18, 5;

%e 41, 88, 79, 40, 8;

%e First three polynomials v(n,x):

%e 1

%e 2 + 2x

%e 5 + 8x + 3x^2

%t u[1, x_] := 1; v[1, x_] := 1; z = 16;

%t u[n_, x_] := u[n - 1, x] + (x + j)*v[n - 1, x] + c;

%t d[x_] := h + x; e[x_] := p + x;

%t v[n_, x_] := d[x]*u[n - 1, x] + e[x]*v[n - 1, x] + f;

%t j = 0; c = 0; h = -1; p = 3; f = 0;

%t Table[Expand[u[n, x]], {n, 1, z/2}]

%t Table[Expand[v[n, x]], {n, 1, z/2}]

%t cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];

%t TableForm[cu]

%t Flatten[%] (* A210803 *)

%t cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];

%t TableForm[cv]

%t Flatten[%] (* A210804 *)

%t Table[u[n, x] /. x -> 1, {n, 1, z}] (* A047849 *)

%t Table[v[n, x] /. x -> 1, {n, 1, z}] (* A000302 *)

%t Table[u[n, x] /. x -> -1, {n, 1, z}] (* A000007 *)

%t Table[v[n, x] /. x -> -1, {n, 1, z}] (* A000007 *)

%Y Cf. A210803, A208510.

%K nonn,tabl

%O 1,2

%A _Clark Kimberling_, Mar 27 2012