A271362
Number T(n,k) of series-reduced free trees with n nodes of which exactly k>=3 are leaves, k+1 <= n <= 2k-2.
Original entry on oeis.org
1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 4, 3, 1, 4, 6, 3, 1, 2, 10, 9, 4, 1, 8, 17, 12, 4, 1, 4, 22, 30, 16, 5, 1, 15, 47, 44, 20, 5, 1, 6, 53, 91, 67, 25, 6, 1, 32, 127, 158, 91, 30, 6, 1, 11, 121, 282, 258, 126
Offset: 4
Irregular triangle begins
n \ k 3 4 5 6 7 8
4 1;
5 1;
6 1, 1;
7 1, 1;
8 1, 2, 1;
9 2, 2, 1;
10 2, 4, 3, 1;
11 4, 6, 3, 1;
12 2, 10, 9, 4, 1;
13 8, 17, 12, 4, 1;
14 4, 22, 30, 16, 5, 1;
15 15, 47, 44, 20, 5, 1;
...
-
\\ using files hitree4.txt etc from McKay.
nL(n, Tr) = { my(E = strsplit(Tr, " "), u_v, Deg = vectorsmall(n));
for(j = 1, n-1, u_v = strsplit(E[j], " "); u_v = eval(u_v);
Deg[ u_v[1]+1 ]++; Deg[ u_v[2]+1 ]++); sum(v = 1, n, Deg[v] == 1)
};
Rows(r1, r2) = {my(F, C, nF); for(n = r1, r2,
F = readstr(Str("hitree", n, ".txt")); C = vectorsmall(n-1);
for(i = 1, #F, nF = nL(n, F[i]); C[nF]++ );
print1(n" "); for(i=1, #C, if(C[i] > 0, print1(C[i]", "))); print() )
}; \\ Washington Bomfim, Jul 09 2021
A345970
Irregular triangle T(n,k) read by rows in which n-th row lists in colex order all series-reduced tree degree sequences D of n nodes encoded as t = Product_{d in D} prime(d); n >= 4, 1 <= k <= A002865(n-2).
Original entry on oeis.org
40, 112, 352, 400, 832, 1120, 2176, 3520, 3136, 4000, 4864, 8320, 9856, 11200, 11776, 21760, 23296, 30976, 35200, 31360, 40000, 29696, 48640, 60928, 73216, 83200, 98560, 87808, 112000, 63488, 117760, 136192, 191488, 173056, 217600, 232960, 309760, 275968, 352000, 313600, 400000
Offset: 4
Triangle begins:
n \ k| 1 2 ... n \ k| 1 2 ...
-----+------------- -----+-----------------------------------
4 | 40; 4 | [3,1,1,1];
5 | 112; 5 | [4,1,1,1,1];
6 | 352, 400; <=> 6 | [5,1,1,1,1,1], [3,3,1,1,1,1];
7 | 832, 1120; 7 | [6,1,1,1,1,1,1], [4,3,1,1,1,1,1];
... ...
Row n = 7 follows from table
.
+---------------------+------------------+---------------------------+
| Partitions of n-2 = | | |
| 5 without parts 1 | Degree sequences | Encoding |
+---------------------+------------------+---------------------------+
| [5] | 6,1,1,1,1,1,1 | prime(6) * 2^6 |
| [2, 3] | 4,3,1,1,1,1,1 | prime(4) * prime(3) * 2^5 |
+---------------------+------------------+---------------------------+
-
Row(n) = {my(j=0, V = vector(numbpart(n-2) - numbpart(n-3)));
forpart(P=n-2, V[j++] = prod(k=1,#P, prime(P[k]+1)) << (n-#P),[2, n-2]); V};
-
Decode(t) = {my(V = [], i = 1, p); while(t > 1, p = prime(i); while(t % p == 0, t /= p; V = concat(V, Vec(i)) ); i++); vecsort(V, (x,y)->y-x) };
Showing 1-2 of 2 results.
Comments