杨氏矩阵

引入

杨氏矩阵(Young tableau),又名杨表,是一种常用于表示论和舒伯特演算中的组合对象.

杨表是一种特殊的矩阵.它便于对称群和一般线性群的群表示和性质研究.杨表由剑桥大学数学家阿尔弗雷德·杨(Alfred Young)于 1900 年首次提出,于 1903 年被德国数学家弗罗贝尼乌斯(Ferdinand Georg Frobenius)应用于对称群的研究.

注释

表示论(Representation theory)是数学的一个分支.它通过将元素表示为向量空间的线性变换来研究抽象代数结构.舒伯特演算(Schubert calculus)是代数几何的一个分支,于 19 世纪由赫尔曼·舒伯特为了解决射影几何的计数问题而引入.

定义

杨图

杨图(Young diagram,使用点表示时又称 Ferrers 图,在 分拆数 一节中有相关介绍)是一个有限的框或单元格集合,左对齐排列,行长按非递增顺序排列.如果把杨图每行的方格数列出,我们得到了一个非负整数 nn(总方格数)的 整数分拆(integer partition)λ\lambda.因此,我们可以将杨图的形状看作 λ\lambda,因为它携带与其整数拆分相同的信息.

杨图之间的包含关系定义了整数分拆上的一个 偏序 关系,此关系拥有 格 的结构,被称为 杨格(Young's lattice).如果把杨图各列的方格数列出,则会得到整数分拆 λ\lambda 的「共轭分拆」,或「转置分拆」,它所对应到的杨图可由原本的杨图沿主对角线作镜射对称而得.

杨图每个方格的位置由分别代表 行数 与 列数 的两个座标点决定.列的顺序由左向右,行的顺序则按方格数的由多向少的方向.此处需要注意,根据习惯不同存在着两种不同的杨图画法:第一个将方格数较少的行排在方格数较多的行的下方,第二种画法将各行由大到小一层一层往上叠.由于前一种画法主要由英语国家使用,而后者通常被法语国家使用,习惯上我们分别称它们为英式画法和法式画法.

以下表格中分别为整数分拆 (5,4,1)(5,4,1) 对应的杨图不同画法:

  • 英式画法:
    配图:young diagram 1
  • 法式画法:
    配图:young diagram 2

杨表

定义

杨表(Young tableau)是通过用取自某个字母表的符号填充杨氏图的框来获得的,这通常需要是一个全序集和.填入的元素写作 x1x_{1},x2x_{2},x3x_{3},…\ldots.但为了方便起见,都直接填入正整数.

杨表最初应用于对称群的表示理论时,允许在杨图的 nn 的方格中任意填入 11 到 nn 中相异的正整数.但现在的研究大多使用「标准」的杨表,即上述条件中各行与各列中方格的数字皆为严格递增的.由 nn 个方格的相异杨表数个数形成 对和数:

注释

对和数(involution number/telephone number)是在数学中是一个整数序列,用来计算 nn 条电话线中每条线路最多可以连接到另一条线路时可以相互连接的方法个数.它还可以用来描述完全图 nn 个顶点上的匹配数,nn 个对合元素的排列数,Hermite 多项式系数的绝对值之和,含有 nn 个格子的标准杨表的个数,以及不可约对称群的度数之和.

1,1,2,4,10,26,76,232,764,2620,9496,…1, 1, 2, 4, 10, 26, 76, 232, 764, 2620, 9496, \ldots(OEIS 中的数列 A000085)

在其他应用中,杨图也可以被填入相同的数字.若填法的同列数字严格递增,且同行数字单调递增,则该杨表被称为是 半标准的(Semistandard Young Tableaux, 有时称为列严格).杨表中个数字出现的次数记录下来得到的序列被视为杨表的 权重.因此,标准杨表的权重必然是 (1,1,…,1)(1,1,\ldots,1).因为在标准杨表中,11 到 nn 的每个正整数恰好各出现一次.

标准杨表的插入算法

排列的性质可以由杨表直观地表现出来.RSK 插入算法 就提供了一个将杨表和排列联系起来的途径.它由 Robinson, Schensted 和 Knuth 提出.

令 SS 是一个杨表,定义 S←xS \leftarrow x 表示将 xx 从第一行插入杨表中,具体如下:

  1. 在当前行中找到最小的比 xx 大的数 yy.
  2. 如果找到了,用 xx 去替换 yy,移到下一行,令 x←yx \leftarrow y 重复操作 1.
  3. 如果找不到,就把 xx 放在该行末尾并退出.记 xx 在第 ss 行第 tt 列,(s,t)(s, t) 必定是一个边角.一个格子 (s,t)(s, t) 是边角当且仅当 (s+1,t)(s + 1, t) 和 (s,t+1)(s, t + 1) 都不存在格子.

例如,将 33 插入杨表 (2,5,9)(6,7)(8)(2, 5, 9)(6, 7)(8) 的步骤为:

配图:young tableau insert

变体

非完全严格标准的杨表有许多变体(Variations).例如行严格杨表要求同行数字严格递增,且同列数字单调递增,即列严格杨表的共轭.此外,在平面分拆(plane partitions)理论中,习惯上会将上述的定义中的递增改为递减.其他变体例如带状杨表,会先将一些方块打包成群,然后要求各群的方块必须填入相同数字.

斜杨表

给定两个杨图 λ=(λ1,λ2…)\lambda = (\lambda{1}, \lambda{2} \ldots)、μ=(μ1,μ2,…)\mu = (\mu{1}, \mu{2},\ldots),满足 λ\lambda 包含 μ\mu,即 μi≤μi\mu{i} \leq \mu{i} 对所有 ii.定义「斜杨图」λ/μ\lambda/\mu 为 λ\lambda 中所有方格减去 μ\mu 中的所有方格 即 λ\lambda 差集 μ\mu.在斜杨图的各方格中填入元素就形成了 斜杨表(Skew tableaux).

例如,下图为整数分拆 (5,4,1)(5,4,1) 对应的一个标准斜杨表:

配图:skew tableau

同理,若满足同一列中的数字严格递增,且同一行中的数字单调递增,则该斜杨表被称作 半标准斜杨表;若半标准斜杨表满足各方格不重复的填入数字 11 到 nn(方格总数),则该斜杨表被称作 标准斜杨表.注意,由不同的 λ\lambda 和 μ\mu 可得到相同的 λ/μ\lambda/\mu.虽然大部分斜杨表的性质都只依赖于取完差集的方格,但是仍然部分运算依赖于 λ\lambda 和 μ\mu 的选取.因此,λ/μ\lambda/\mu 必须被视为包含两个元素信息:λ\lambda 和 μ\mu.当 μ\mu 是空分拆(00 的唯一一种分拆)时,斜杨表 λ/μ\lambda/\mu 就变成杨表 λ\lambda.

应用

杨表常用于在组合学、表示理论和代数几何中,用各种不同计算杨表个数的方法得到舒尔函数的定义及相关的恒等式.在信息学竞赛中,常有考察杨表勾长公式的题目.

勾长

给定一个共有 nn 个方格的杨表 πλ\pi_{\lambda},把 11 到 nn 的 nn 个数字填入杨表中,使得每行从左到右,每列从下到上都是递增的.用 dimπλdim_{\pi_{\lambda}} 表示可以这样填的方法个数.

对于杨表中的一个方格 vv,定义其 勾长 hook(v)\mathrm{hook}(v) 等于同行右边的方格数加上同列上面的方格数,再加 1(即方格本身).

勾长公式

如果用 dimλdim_{\lambda} 表示这样的方法个数,勾长公式 就是方法个数等于 n!n! 除以所有方格的勾长的乘积.

dim⁡πλ=n!∏x∈Y(λ)hook(x).\dim \pi _{\lambda}={\frac {n!}{\prod_{{x\in Y(\lambda)}}{\mathrm {hook}}(x)}}.
配图:young tableau 2

所以对于整数分拆 10=5+4+110 = 5 + 4 + 1 的杨表,如上图所示.有

dim⁡πλ=10!7⋅5⋅4⋅3⋅1⋅5⋅3⋅2⋅1⋅1=288.\dim \pi _{\lambda }={\frac {10!}{7\cdot 5\cdot 4\cdot 3\cdot 1\cdot 5\cdot 3\cdot 2\cdot 1\cdot 1}}=288.

种方法.

例题

子序列问题

对于杨表 PP, 定义对于一个从 11 到 nn 的排列 X=x1,…,xnX = x_{1}, \ldots , x_{n}.

  1. PXP_{X} 中第一行的长度即为排列 XX 的 最长上升子序列(LIS) 长度.注意,PP 的第一行并不一定是 LIS 本身,所以不能直接利用杨表性质解决「LIS 划分」之类的问题.

  2. 对于一个排列 XX 和它产生的杨表 PXP_{X},若 XRX^R 是 XX 的翻转,那么 XRX^R 产生的杨表 PXRP_{X^R} 即为 PXP_{X} 交换行列得到.

    例如,对于排列 X=1,5,7,2,8,6,3,4X = 1, 5, 7, 2, 8, 6, 3, 4 和 XR=4,3,6,8,2,7,5,1X^R = 4, 3, 6, 8, 2, 7, 5, 1, 我们可得到如下杨表 PXP_{X}:

    配图:young tableau LIS
  3. 杨表 PXP_{X} 中的第一列长度即为排列 XX 的 最长下降子序列(LDS) 长度.

定义长度不超过 kk 的 LIS/LDSLIS/LDS 长度为 k−LISk-LIS 和 k−LDSk-LDS, 此类问题我们同样可以用杨表来解决.对于 1−LIS1-LIS,显而易见最长的 1−LIS1-LIS 子序列就是该序列的 LDSLDS,这也正是杨表的第一列;同样可得,杨表前 kk 列的长度就是最长的 k−LISk-LIS 子序列的长度.证明如下:

对于一个排列 XX 和它的 mm 行杨表 PP,令排列 X∗X^∗ 为 (Pm,1…,Pm,λm,Pm−1,1…,P1,1…P1,λ1)(P_{m,1}\ldots,P_{m,λ_{m}},P_{m−1,1}\ldots,P_{1,1}\ldots P_{1,λ_{1}})(即将杨表从下往上每行依次写在后面).那么 XX 一定可以通过交换操作转化成 X∗X^*.

所以,最长 k−LISk-LIS 子序列长度可以表示成 F(k)=∑i=1mmin⁡(k,λi)F(k)=\sum_{i=1}^{m} \min(k,\lambda_{i}),即前 kk 列的长度和.

CTSC2017 最长上升子序列

有一个长为 nn 的数列 bb.对于序列 Bm=(b1,b2,…,bm)B_{m} = (b_{1}, b_{2},\ldots, b_{m}),设 CC 是 BmB_{m} 的子序列,且 CC 的最长上升子序列的长度不超过 kk,询问 CC 的长度最大值.

解题思路

多个询问考虑使用扫描线的方法.这样我们就需要维护每个前缀的杨表.如果使用以上结论,可以发现问题变成了如何快速维护杨表前 kk 列的长度之和.如果直接维护,复杂度是 O(n2log⁡n)O(n^2 \log n) 的不能接受.考虑维护前 n\sqrt{n} 列和前 n\sqrt{n} 行.

可以发现,杨表一定不会完全覆盖这个 W×HW \times H 的矩形.如果 K≤WK \leq W,那么可以直接得答案;如果 K>WK > W,那么大于 WW 的部分一定在 HH 行内.所以可以考虑如何同时维护前 n\sqrt{n} 列和前 n\sqrt{n} 行.将这个排列翻转一下就可以得到杨表的翻转,所以只需要再同时维护 −Ai-A_{i} 即可,复杂度为 O(nnlog⁡n)O(n \sqrt{n} \log n).

BJWC2018 最长上升子序列

现在有一个长度为 nn 的随机排列,求它的最长上升子序列长度的期望.

CF1268B 杨氏多米诺骨牌

给定一个具有 nn 列长度 a1,a2,…,an (a1≥a2≥…≥an≥1)a_{1} ,a_{2},\ldots,a_{n}\,(a_{1} \geq a_{2} \geq \ldots \geq a_{n} \geq 1) 的直方图.a=[3,2,2,2,1]a=[3,2,2,2,1] 的杨图.找到可以在此直方图中绘制的最大数量的非重叠多米诺骨牌(1×21 \times 2 或 2×12 \times 1 矩形).

参考资料与拓展阅读

  1. Young Tableau - from Wolfram MathWorld
  2. Young tableau - Wikipedia
  3. Hook length formula - Wikipedia
  4. 袁方舟,《浅谈杨氏矩阵在信息学竞赛中的应用》IOI2019, 中国国家候选队论文集,202-229

本文内容来自 OI Wiki《杨氏矩阵》;作者与修订历史见原文。AIOJ 对排版、站内链接和资源路径做了格式适配。

内容依 CC BY-SA 4.0 与 SATA 许可共享;附加条款以原文声明为准。感谢 OI Wiki 作者与贡献者;请按 SATA 约定前往 OI Wiki 仓库 Star 支持。