OFFSET
0,8
FORMULA
G.f. A_k(x) of column k satisfies A_k(x) = ( 1 + x * A_k(x)^(2/k) * (1 + x * 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 A001002.
B(x)^k = B(x)^(k-1) + x * B(x)^(k+1) + x^2 * B(x)^(k+2). So T(n,k) = T(n,k-1) + T(n-1,k+1) + T(n-2,k+2) for n > 1.
EXAMPLE
Square array begins:
1, 1, 1, 1, 1, 1, 1, ...
0, 1, 2, 3, 4, 5, 6, ...
0, 3, 7, 12, 18, 25, 33, ...
0, 10, 26, 49, 80, 120, 170, ...
0, 38, 105, 210, 363, 575, 858, ...
0, 154, 444, 927, 1672, 2761, 4290, ...
0, 654, 1944, 4191, 7810, 13325, 21385, ...
PROG
(PARI) T(n, k, t=2, u=1) = if(k==0, 0^n, k*sum(r=0, n, binomial(t*r+u*(n-r)+k, r)*binomial(r, n-r)/(t*r+u*(n-r)+k)));
matrix(7, 7, n, k, T(n-1, k-1))
CROSSREFS
KEYWORD
AUTHOR
Seiichi Manyama, Nov 21 2024
STATUS
approved