OFFSET
0,5
EXAMPLE
Triangle begins:
k=0 1 2 3 4 5 6 7
n=0 [1]
n=1 [0, 1]
n=2 [0, 2, 2]
n=3 [0, 5, 6, 6]
n=4 [0, 14, 22, 36, 24]
n=5 [0, 42, 90, 150, 240, 120]
n=6 [0, 132, 378, 648, 1560, 1800, 720]
n=7 [0, 429, 1638, 3318, 8400, 16800, 15120, 5040]
...
T(3,2) = 6 counts the three leaf permutations of each of the following trees:
__o__ __o__
/ | \ / | \
(1) (1) (2) (1) (2) (2)
PROG
(PARI)
subsets(S) = {my(s=List()); for(i=0, 2^(#S) -1, my(x=List()); for(j=1, #S, if(bitand(i, 1<<(j-1)), listput(x, S[j]))); listput(s, Vec(x))); Vec(s)}
C_aB(B) = {my(S = subsets(B)); sum(i=1, #S, (1/(1-x*z*#S[i]))*(-1)^(#B-#S[i]))}
D(k, N, B) = {if(k>N, 1, substpol(C_aB(B), z, 1 + D(k+1, N-#B+1, B)))}
Dx(N, B) = {Vec(1+D(1, N, B)+ O('x^(N+1)))}
T(max_row) = {my( N = max_row+1, v = vector(N, i, if(i==1, 1, 0))~); for(k=1, N, v=matconcat([v, Dx(N+1, vector(k, i, i))~])); vector(N, n, vector(n, k, v[n, k]))}
T(8)
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
John Tyler Rascoe, Jun 18 2025
STATUS
approved
