OFFSET
0,5
LINKS
Seiichi Manyama, Antidiagonals n = 0..139, flattened
FORMULA
G.f. A_k(x) of column k satisfies A_k(x) = 1/(1 - x) + k*x*A_k(x)^2.
A_k(x) = 2/( 1 - x + sqrt((1 - x) * (1 - (4*k+1)*x)) ).
T(n,k) = 1 + k * Sum_{j=0..n-1} T(j,k) * T(n-1-j,k).
(n+1) * T(n,k) = 2 * ((2*k+1) * n - k) * T(n-1,k) - (4*k+1) * (n-1) * T(n-2,k) for n > 1.
E.g.f. of column k: exp((2*k+1)*x) * (BesselI(0,2*k*x) - BesselI(1,2*k*x)). - Ilya Gutkovskiy, Feb 01 2021
EXAMPLE
Square array begins:
1, 1, 1, 1, 1, 1, ...
1, 2, 3, 4, 5, 6, ...
1, 5, 13, 25, 41, 61, ...
1, 15, 71, 199, 429, 791, ...
1, 51, 441, 1795, 5073, 11571, ...
1, 188, 2955, 17422, 64469, 181776, ...
MATHEMATICA
T[n_, k_] := Sum[If[j == k == 0, 1, k^j] * Binomial[n, j] * CatalanNumber[j], {j, 0, n}]; Table[T[k, n - k], {n, 0, 9}, {k, 0, n}] // Flatten (* Amiram Eldar, Feb 01 2021 *)
PROG
(PARI) {T(n, k) = sum(j=0, n, k^j*binomial(n, j)*(2*j)!/(j!*(j+1)!))}
(PARI) {T(n, k) = 1+k*sum(j=0, n-1, T(j, k)*T(n-1-j, k))}
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Seiichi Manyama, Jan 31 2021
STATUS
approved