|
|
A256042
|
|
Triangle read by rows: number of idempotent basis elements of rank k in partition monoid P_n.
|
|
0
|
|
|
1, 0, 1, 0, 5, 1, 0, 43, 15, 1, 0, 529, 247, 30, 1, 0, 8451, 4795, 805, 50, 1, 0, 167397, 108871, 22710, 1985, 75, 1, 0, 3984807, 2855279, 697501, 76790, 4130, 105, 1, 0, 111319257, 85458479, 23520966, 3070501, 209930, 7658, 140, 1, 0, 3583777723, 2887069491, 871103269, 129732498, 10604811, 495054, 13062, 180, 1
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,5
|
|
COMMENTS
|
|
|
LINKS
|
|
|
EXAMPLE
|
Triangle begins:
1,
0, 1,
0, 5, 1,
0, 43, 15, 1,
0, 529, 247, 30, 1,
0, 8451, 4795, 805, 50, 1,
0, 167397, 108871, 22710, 1985, 75, 1,
0, 3984807, 2855279, 697501, 76790, 4130, 105, 1,
0, 111319257, 85458479, 23520966, 3070501, 209930, 7658, 140, 1,
...
|
|
MATHEMATICA
|
rows = 10;
f[n_, r_, s_] := f[n, r, s] = Module[{resu, m, a, b}, Which[n <= 0, 0, s == 1, StirlingS2[n, r], r == 1, StirlingS2[n, s], True, resu = s*f[n - 1, r - 1, s] + r*f[n - 1, r, s - 1] + r*s*f[n - 1, r, s]; Do[resu += Binomial[n - 2, m]*(b*(r - a) + a*(s - b))*f[m, a, b]*f[-m + n - 1, r - a, s - b], {m, n}, {a, r - 1}, {b, s - 1}]; resu]];
a33[n_] := Module[{b = 0}, Do[b += r*s*f[n, r, s], {r, n}, {s, n}]; b];
BellMatrix[f_, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len - 1}, {k, 0, len - 1}]];
M = BellMatrix[a33[# + 1]&, rows];
|
|
PROG
|
(Sage) # uses[bell_matrix from A264428]
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|