A210753 Triangle of coefficients of polynomials u(n,x) jointly generated with A210754; see the Formula section.
1, 2, 2, 3, 7, 4, 4, 16, 20, 8, 5, 30, 61, 52, 16, 6, 50, 146, 198, 128, 32, 7, 77, 301, 575, 584, 304, 64, 8, 112, 560, 1408, 1992, 1616, 704, 128, 9, 156, 966, 3060, 5641, 6328, 4272, 1600, 256, 10, 210, 1572, 6084, 14002, 20330, 18880, 10912, 3584
Offset: 1
Examples
First five rows: 1 2...2 3...7....4 4...16...20...8 5...30...61...52...16 First three polynomials u(n,x): 1, 2 + 2x, 3 + 7x + 4x^2.
Programs
-
Mathematica
u[1, x_] := 1; v[1, x_] := 1; z = 16; u[n_, x_] := (x + 1)*u[n - 1, x] + x*v[n - 1, x] + 1; v[n_, x_] := (x + 1)*u[n - 1, x] + (x + 1)*v[n - 1, x] + 1; Table[Expand[u[n, x]], {n, 1, z/2}] Table[Expand[v[n, x]], {n, 1, z/2}] cu = Table[CoefficientList[u[n, x], x], {n, 1, z}]; TableForm[cu] Flatten[%] (* A210753 *) Table[Expand[v[n, x]], {n, 1, z}] cv = Table[CoefficientList[v[n, x], x], {n, 1, z}]; TableForm[cv] Flatten[%] (* A210754 *) Table[u[n, x] /. x -> 1, {n, 1, z}] (* A007070 *) Table[v[n, x] /. x -> 1, {n, 1, z}] (* A035344 *)
Formula
u(n,x)=(x+1)*u(n-1,x)+x*v(n-1,x)+1,
v(n,x)=(x+1)*u(n-1,x)+(x+1)*v(n-1,x)+1,
where u(1,x)=1, v(1,x)=1.
Comments