OFFSET
0,5
FORMULA
G.f. A_k(x) of column k satisfies A_k(x) = ( 1 + x * A_k(x)^(3/k) * (1 + A_k(x)^(1/k)) )^k for k > 0.
G.f. of column k: B(x)^k where B(x) is the g.f. of A144097.
B(x)^k = B(x)^(k-1) + x * B(x)^(k+2) + x * B(x)^(k+3). So T(n,k) = T(n,k-1) + T(n-1,k+2) + T(n-1,k+3) for n > 0.
EXAMPLE
Square array begins:
1, 1, 1, 1, 1, 1, 1, ...
0, 2, 4, 6, 8, 10, 12, ...
0, 14, 32, 54, 80, 110, 144, ...
0, 134, 324, 578, 904, 1310, 1804, ...
0, 1482, 3696, 6810, 11008, 16490, 23472, ...
0, 17818, 45316, 85278, 140936, 216002, 314700, ...
0, 226214, 583152, 1113854, 1870352, 2914790, 4320608, ...
PROG
(PARI) T(n, k, t=3, u=1) = if(k==0, 0^n, k*sum(r=0, n, binomial(n, r)*binomial(t*n+u*r+k, n)/(t*n+u*r+k)));
matrix(7, 7, n, k, T(n-1, k-1))
CROSSREFS
KEYWORD
AUTHOR
Seiichi Manyama, Nov 20 2024
STATUS
approved