login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A371259
Triangle read by rows, (2, 3)-Lah numbers.
5
1, 36, 1, 1764, 100, 1, 112896, 9864, 200, 1, 9144576, 1099296, 34064, 344, 1, 914457600, 142159392, 6004512, 92200, 540, 1, 110649369600, 21385410048, 1156921920, 24075712, 213700, 796, 1
OFFSET
3,2
COMMENTS
The (2, 3)-Lah numbers T(n, k) count ordered 2-tuples (pi(1), pi(2)) of partitions of the set {1, ..., n} into k linearly ordered blocks (lists, for short) such that the numbers 1, 2, 3 are in distinct lists, and bl(pi(1)) = bl(pi(2)) where for i = {1, 2} and pi(i) = b(1)^i, b(2)^i, ..., b(k)^i, where b(1)^i, b(2)^i, ..., b(k)^i are the blocks of partition pi(i), bl(pi(i)) = {min(b(1))^i, min(b(2))^i, ..., min(b(k))^i} is the set of block leaders, i.e., of minima of the lists in partition pi(i).
The (2, 3)-Lah numbers T(n, k) are the (m, r)-Lah numbers for m=2 and r=3. More generally, the (m, r)-Lah numbers count ordered m-tuples (pi(1), pi(2), ..., pi(m)) of partitions of the set {1, 2, ..., n} into k linearly ordered blocks (lists, for short) such that the numbers 1, 2, ..., r are in distinct lists, and bl(pi(1)) = bl(pi(2)) = ... = bl(pi(m)) where for i = {1, 2, ..., m} and pi(i) = {b(1)^i, b(2)^i, ..., b(k)^i}, where b(1)^i, b(2)^i,..., b(k)^i are the blocks of partition pi(i), bl(pi(i)) = {min(b(1))^i, min(b(2))^i, ..., min (b(k))^i} is the set of block leaders, i.e., of minima of the lists in partition pi(i).
LINKS
A. Žigon Tankosič, The (l, r)-Lah Numbers, Journal of Integer Sequences, Article 23.2.6, vol. 26 (2023).
FORMULA
Recurrence relation: T(n, k) = T(n-1, k-1) + (n+k-1)^2*T(n-1, k).
Explicit formula: T(n, k) = Sum_{4 <= j(1) < j(2) < ... < j(n-k) <= n} (2j(1)-2)^2 * (2j(2)-3)^2 * ... * (2j(n-k)-(n-k+1))^2.
Special cases:
T(n, k) = 0 for n < k or k < 3, [corrected by Paolo Xausa, Jun 11 2024]
T(n, n) = 1,
T(n, 3) = (A143498(n, 3))^2 = ((n+2)!)^2/14400,
T(n, n-1) = 2^2 * Sum_{j=3..n-1} j^2.
EXAMPLE
Triangle begins:
1;
36, 1;
1764, 100, 1;
112896, 9864, 200, 1;
9144576, 1099296, 34064, 344, 1;
914457600, 142159392, 6004512, 92200, 540, 1;
110649369600, 21385410048, 1156921920, 24075712, 213700, 796, 1.
...
An example for T(4, 3). The corresponding partitions are
pi(1) = {(1),(2),(3,4)},
pi(2) = {(1),(2),(4,3)},
pi(3) = {(1),(3),(2,4)},
pi(4) = {(1),(3),(4,2)},
pi(5) = {(1,4),(2),(3)},
pi(6) = {(4,1),(2),(3)}, since A143498 for n=4, k=3 equals 6. Sets of their block leaders are bl(pi(1)) = bl(pi(2)) = bl(pi(3)) = bl(pi(4)) = bl(pi(5)) = bl(pi(6)) = {1,2,3}.
Compute the number of ordered 2-tuples (i.e., ordered pairs) of partitions pi(1), pi(2), ..., pi(6) such that partitions in the same pair share the same set of block leaders. As there are six partitions with the set of block leaders equal to {1,2,3}, T(4, 3) = 6^2 = 36.
MAPLE
T:= proc(n, k) option remember; `if`(k<3 or k>n, 0,
`if`(n=k, 1, T(n-1, k-1)+(n+k-1)^2*T(n-1, k)))
end:
seq(seq(T(n, k), k=3..n), n=3..10);
MATHEMATICA
A371259[n_, k_] := A371259[n, k] = Which[n < k || k < 3, 0, n == k, 1, True, A371259[n-1, k-1] + (n+k-1)^2*A371259[n-1, k]];
Table[A371259[n, k], {n, 3, 10}, {k, 3, n}] (* Paolo Xausa, Jun 11 2024 *)
PROG
(Python)
def T_Lah(n, k):
if k < 3 or k > n:
return 0
elif n == k == 3:
return 1
else:
return T_Lah(n-1, k-1) + ((n+k-1)**2) * T_Lah(n-1, k)
def print_triangle(rows):
for n in range(rows):
row_values = [T_Lah(n, k) for k in range(n+1)]
print(' '.join(map(str, row_values)).center(rows*10))
rows = 10
print_triangle(rows)
CROSSREFS
Sequence in context: A330084 A350385 A374498 * A037935 A159824 A285575
KEYWORD
nonn,tabl
AUTHOR
STATUS
approved