login
Rectangular array: (row n) = b**c, where b(h) = h, c(h) = 2*(n-1+h)-1, n>=1, h>=1, and ** = convolution.
5

%I #4 Jun 21 2012 22:28:36

%S 1,5,3,14,11,5,30,26,17,7,55,50,38,23,9,91,85,70,50,29,11,140,133,115,

%T 90,62,35,13,204,196,175,145,110,74,41,15,285,276,252,217,175,130,86,

%U 47,17,385,375,348,308,259,205,150,98,53,19,506,495,465,420

%N Rectangular array: (row n) = b**c, where b(h) = h, c(h) = 2*(n-1+h)-1, n>=1, h>=1, and ** = convolution.

%C Principal diagonal: A007585

%C Antidiagonal sums: A002417

%C row 1, (1,2,3,4,5,...)**(1,3,5,7,9,...): A000330

%C row 2, (1,2,3,4,5,...)**(3,5,7,9,...): A051925

%C row 3, (1,2,3,4,5,...)**(5,7,9,11,...): (2*k^3 + 15*k^2 + 13*k)/6

%C row 4, (1,2,3,4,5,...)**(7,9,11,13,...): (2*k^3 + 21*k^2 + 19*k)/6

%C For a guide to related arrays, see A213500.

%F T(n,k) = 4*T(n,k-1)-6*T(n,k-2)+4*T(n,k-3)-T(n,k-4).

%F G.f. for row n: f(x)/g(x), where f(x) = (2*n - 1) - (2*n - 3)*x and g(x) = (1 - x )^4.

%e Northwest corner (the array is read by falling antidiagonals):

%e 1....5....14...30....55....91

%e 3....11...26...50....85....133

%e 5....17...38...70....115...175

%e 7....23...50...90....145...217

%e 9....29...62...110...175...259

%e 11...35...74...130...205...301

%t b[n_] := n; c[n_] := 2 n - 1;

%t t[n_, k_] := Sum[b[k - i] c[n + i], {i, 0, k - 1}]

%t TableForm[Table[t[n, k], {n, 1, 10}, {k, 1, 10}]]

%t Flatten[Table[t[n - k + 1, k], {n, 12}, {k, n, 1, -1}]]

%t r[n_] := Table[t[n, k], {k, 1, 60}] (* A213750 *)

%t d = Table[t[n, n], {n, 1, 40}] (* A007585 *)

%t s1 = Table[s[n], {n, 1, 50}] (* A002417 *)

%t FindLinearRecurrence[s1]

%t FindGeneratingFunction[s1, x]

%Y Cf. A213500.

%K nonn,tabl,easy

%O 1,2

%A _Clark Kimberling_, Jun 20 2012