cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A361032 Square array read by ascending antidiagonals: T(n,k) = F(n) * (4*k)!/(k!*(k + n + 1)!^3), where F(n) = (1/8)*(4*n + 4)!/(n + 1)!; n, k >= 0.

This page as a plain text file.
%I A361032 #16 Jan 05 2024 14:27:22
%S A361032 3,315,9,46200,280,280,7882875,17325,3675,17325,1466593128,1513512,
%T A361032 116424,116424,1513512,288592936632,162954792,5885880,2134440,5885880,
%U A361032 162954792,59064793444800,20193091776,399072960,67953600,67953600,399072960,20193091776,12445136556298875
%N A361032 Square array read by ascending antidiagonals: T(n,k) = F(n) * (4*k)!/(k!*(k + n + 1)!^3), where F(n) = (1/8)*(4*n + 4)!/(n + 1)!; n, k >= 0.
%C A361032 The Catalan numbers A000108 are given by the formula Catalan(k) = (2*k)!/(k!*(k + 1)!). Gessel (1992) considered generalized Catalan numbers defined by Catalan(n,k) = J(n) * (2*k)!/(k!*(k + n + 1)!), where J(n) = (2*n + 2)!/(2*(n + 1)!) = (2^n)*Product_{j = 0..n} (2*j + 1) is chosen so that these numbers are always integers. Gessel's generalized Catalan numbers are particular cases of super ballot numbers. See A135573 for a table of these generalized Catalan numbers.
%C A361032 For this table we carry out an analogous construction using the numbers B(k) = (4*k)!/k!^4 = A008977(k) in place of the central binomial numbers (2*k)!/k!^2. We define sequences {B(n,k) : k >= 0}, n = 0, 1, 2, ..., by B(n,k) = F(n) * (4*k)!/(k!*(k + n + 1)!^3), where choosing F(n) = (4*n + 4)!/(8*(n + 1)!) = (1/2)*(4^n)*Product_{j = 0..n} (4*j + 1)*(4*j + 2)*(4*j + 3) appears to produce integer values for these quantities. The rows of the square array below are the sequences  {B(0,k)}, {B(1,k)}, {B(2,k)}, ....
%C A361032 An alternative expression is B(n,k) =  G(n,k) * B(n+k+1), where G(n,k) = (1/8)*Product_{j = 0..n} ( (4*j + 1)*(4*j + 2)*(4*j + 3)/((4*k + 4*j + 1)*(4*k + 4*j + 2)*(4*k + 4*j + 3)) ).
%H A361032 Andrew Howroyd, <a href="/A361032/b361032.txt">Table of n, a(n) for n = 0..1325</a> (first 51 antidiagonals)
%H A361032 Ira M. Gessel, <a href="http://people.brandeis.edu/~gessel/homepage/papers/superballot.pdf">Super ballot numbers</a>, J. Symbolic Comp., 14 (1992), 179-194.
%F A361032 T(n,k) = (-1)^k*(1/8)*256^(n+k+1)*binomial(n+1/4, n+k+1)*binomial(n+1/2, n+k+1)* binomial(n+3/4, n+k+1).
%F A361032 P-recursive: (n + k + 1)^3*T(n,k) = 4*(4*k - 1)*(4*k - 2)*(4*k - 3)*T(n,k-1) with T(n,0) = (1/8)*(4*n + 4)!/(n + 1)!^4 = (1/8)*A008977(n+1).
%F A361032 (n + k + 1)^3*T(n,k) = 4*(4*n + 1)*(4*n + 2)*(4*n + 3)*T(n-1,k) with T(0,k) = 3*(4*k)!/(k!*(k+1)!^3) = A361033(k).
%F A361032 T(n,k) = (1/2) * (1/(2*Pi))^3 * 256^(n+k+1) * Integral_{x = 0..1} (1 - x)^(n+1/4)*x^(k-1/4) dx * Integral_{x = 0..1} (1 - x)^(n+1/2)*x^(k-1/2) dx * Integral_{x = 0..1} (1 - x)^(n+3/4)*x^(k-3/4) dx.
%e A361032 The square array with rows n >= 0 and columns k >= 0 begins:
%e A361032   n\k|          0          1          2           3           4 ...
%e A361032   ----------------------------------------------------------------------
%e A361032    0 |          3          9        280       17325     1513512 ...
%e A361032    1 |        315        280       3675      116424     5885880 ...
%e A361032    2 |      46200      17325     116424     2134440    67953600 ...
%e A361032    3 |    7882875    1513512    5885880    67953600  1449322875 ...
%e A361032    4 | 1466593128  162954792  399072960  3086579925 46235189000 ...
%e A361032    5 |  ...
%e A361032   ...
%e A361032 As a triangle:
%e A361032  Row
%e A361032   0 |             3
%e A361032   1 |           315          9
%e A361032   2 |         46200        280      280
%e A361032   3 |       7882875      17325     3675    17325
%e A361032   4 |    1466593128    1513512   116424   116424  1513512
%e A361032   5 |  288592936632  162954792  5885880  2134440  5885880  162954792
%p A361032 # as a square array
%p A361032 T := proc (n, k) (-1)^k*(1/8)*256^(n+1+k)*binomial(n+1/4, n+1+k)*binomial(n+2/4, n+1+k)* binomial(n+3/4, n+1+k); end proc:
%p A361032 for n from 0 to 10 do seq(T(n, k), k = 0..10); end do;
%p A361032 # as a triangle
%p A361032 T := proc (n, k) (-1)^k*(1/8)*256^(n+1+k)*binomial(n+1/4, n+1+k)*binomial(n+2/4, n+1+k)* binomial(n+3/4, n+1+k); end proc:
%p A361032 for n from 0 to 10 do seq(T(n-k, k), k = 0..n); end do;
%o A361032 (PARI) T(n,k) = (-1)^k*(1/8)*256^(n+k+1)*binomial(n+1/4, n+k+1)*binomial(n+1/2, n+k+1)* binomial(n+3/4, n+k+1) \\ _Andrew Howroyd_, Jan 05 2024
%Y A361032 A361033 (row 0), A361034 (row 2), A361035 (row 3).
%Y A361032 Cf. A008977, A135573, A361027.
%K A361032 nonn,tabl,easy
%O A361032 0,1
%A A361032 _Peter Bala_, Mar 04 2023