%I #8 May 05 2023 07:47:00
%S 1,1,0,1,2,0,1,6,24,0,1,12,120,720,0,1,20,360,5040,40320,0,1,30,840,
%T 20160,362880,3628800,0,1,42,1680,60480,1814400,39916800,479001600,0,
%U 1,56,3024,151200,6652800,239500800,6227020800,87178291200,0,1,72,5040,332640,19958400,1037836800,43589145600,1307674368000,20922789888000,0
%N Triangle read by rows, T(n, k) = (-1)^k * RisingFactorial(n, k) * FallingFactorial(k - n, k).
%F T(n, k) = Pochhammer(n, k) * Pochhammer(n - k, k).
%F T(n, k) = Gamma(n + k) / Gamma(n - k) if k != n.
%F T(n, k) = (-1)^k * binomial(k - n, k) * binomial(n + k - 1, k) * (k!)^2.
%F T(n, k) = binomial(n + k - 1, n - 1) * binomial(n - 1, n - 1 - k) *(k!)^2.
%e Triangle T(n, k) starts:
%e [0] 1;
%e [1] 1, 0;
%e [2] 1, 2, 0;
%e [3] 1, 6, 24, 0;
%e [4] 1, 12, 120, 720, 0;
%e [5] 1, 20, 360, 5040, 40320, 0;
%e [6] 1, 30, 840, 20160, 362880, 3628800, 0;
%e [7] 1, 42, 1680, 60480, 1814400, 39916800, 479001600, 0;
%e [8] 1, 56, 3024, 151200, 6652800, 239500800, 6227020800, 87178291200, 0;
%p T := (n, k) -> if n = k then 0^k else GAMMA(n + k) / GAMMA(n - k) fi:
%p seq(seq(T(n, k), k = 0..n), n = 0..9);
%t Table[Pochhammer[n, k]*Pochhammer[n - k, k], {n, 0, 9}, {k, 0, n}] // Flatten (* _Michael De Vlieger_, May 05 2023 *)
%Y Cf. A362846 (row sums), A010050 (main diagonal), A002378 (column 1).
%K nonn,tabl
%O 0,5
%A _Peter Luschny_, May 05 2023