|
%I #18 Nov 25 2022 11:48:29
%S 1,4,1,5,4,1,24,11,4,1,35,52,18,4,1,192,123,84,26,4,1,315,660,285,120,
%T 35,4,1,1920,1545,1500,545,160,45,4,1,3465,9180,4680,2820,930,204,56,
%U 4,1,23040,22005,27180,11220,4740,1470,252,68,4,1,45045,142380
%N Triangular array: Row n shows the coefficients of polynomials p(n,x) defined in Comments.
%C The polynomial p(n,x) is the numerator of the rational function given by f(n,x) = x + (n + 1)/f(n-1,x), where f(0,x) = 1.
%C Row sums give A249059(n) for n >= 1.
%C First column is A249060 (n-th term = n!! for n >= 0).
%H Clark Kimberling, <a href="/A249057/b249057.txt">Table of n, a(n) for n = 0..5049</a>
%e f(0,x) = 1/1, so that p(0,x) = 1
%e f(1,x) = (4 + x)/1, so that p(1,x) = 4 + x;
%e f(2,x) = (5 + 4 x + x^2)/(1 + x), so that p(2,x) = 5 + 4 x + x^2.
%e First 6 rows of the triangle of coefficients:
%e 1
%e 4 1
%e 5 4 1
%e 24 11 4 1
%e 35 52 18 4 1
%e 192 123 84 26 4 1
%t z = 12; f[x_, n_] := x + (n+3)/f[x, n - 1];
%t f[x_, 0] = 1; t = Table[Factor[f[x, n]], {n, 0, z}];
%t u = Numerator[t]; TableForm[Rest[Table[CoefficientList[u[[n]], x], {n, 0, z}]]];
%t Flatten[CoefficientList[u, x]] (* A249057 sequence *)
%o (PARI) f(n) = if (n, x + (n + 3)/f(n-1), 1);
%o row(n) = Vecrev(numerator(f(n)), 0); \\ _Michel Marcus_, Nov 25 2022
%Y Cf. A249059, A249060.
%K nonn,tabl,easy
%O 0,2
%A _Clark Kimberling_, Oct 20 2014
|
