符号化方法
解析组合试图从一个较为机械化的方式帮助我们将组合计数问题从模型直接转为生成函数。
—— $\text{EntropyIncreaser}$
解析组合为我们提供了一套能处理一系列组合结构计数和渐进估计的方法，分为解析和组合两个部分。
解析的部分主要讨论了如何近似生成函数的系数，而组合的部分则着眼于操作组合结构以便计数。符号化方法就是组合部分中操作组合结构的一种机械方法，通过描述组合结构来声明对应的生成函数。通过这种方式，我们就可以跳过转移方程，直接写出生成函数，也能够得到关于生成函数和组合意义的联系的更深刻理解。
阅读下文可能需要有基础的生成函数知识。也许是前置知识。
本文目前的内文字符数约为 30.8k，可能需要一点时间加载 \(\LaTeX\)。
注意本文内容可能有疏漏或错误，如果读者发现敬请在评论区指出不足。
\(1.\) 无标号体系
\(1.1\) 记号与定义
组合对象是满足某一性质的可数对象，例如二叉树、烯烃或字符串。组合类是一系列组合对象组成的可数集合。
下面将采用美术花体（\mathcal）大写字母表示一个组合类。
每个组合类 \(\mathcal A\) 都定义了一个大小函数 \(f: \mathcal A\to \mathbb N\)，满足对任何 \(k\in \mathbb N\) 只有有限个组合对象 \(a\in \mathcal A\) 满足 \(f(a) = k\)。\(f(a)\) 表示了组合对象 \(a\) 的固有属性，可能是树高、节点数或串长等。常记 \(f(a)\) 为 \(\lvert a\rvert\)（在某些需要指明组合类的场合下记作 \(\lvert a \rvert_{\mathcal A}\)）。
我们将 \(a\in \mathcal A \text{ s.t. } \lvert a \rvert = n\) 的全体记作 \(\mathcal A_n\)，并以此构造一个计数序列 \(A\)，满足 \(A[n] = \lvert \mathcal A_n\rvert\)。
两个组合类 \(\mathcal A, \mathcal B\) 在组合意义上同构记作 \(\mathcal A = \mathcal B\) 或 \(\mathcal A \cong \mathcal B\)，但仅在 \(\mathcal A\) 和 \(\mathcal B\) 不同构时才用后一种记号。
我们在描述组合对象时只会关注我们需要的单一性质，这种抽象有助于我们理解下文中将不同对象组合为新的对象的操作。总的来说，我们忽视组合对象除了“大小”外的所有信息，这样很多组合操作都可以被简单地映射在“并列”操作上。
一个组合类 \(\mathcal A\) 的 \(\text{OGF}\) 即其对应的计数序列 \(A\) 的 \(\text{OGF}\)，为一个形式幂级数
它也可以被等价地写作
我们称 \(z^{\lvert a \rvert}\) 为 \(a\) 的幂表示，这里的 \(z\) 仍然只是占位元。
由以上内容不难验证组合类的笛卡尔积和 \(\text{OGF}\) 的乘法 \(\times\) 同构，集合不交并和 \(\text{OGF}\) 的加法 \(+\) 同构（处理了重复元素的加法会在下方重新定义）。
下面对集合的笛卡尔积也记作 \(\times\)，集合并也记作 \(+\)。
组合类都是由一些更本质的组合类构造出来的；一个构造是从一组组合类映射到一个组合类的函数。
当我们说一个构造可以被翻译为生成函数运算时，我们想要表达的就是这种组合类间的映射能和生成函数间的运算对应。一个经典的例子是 \(\exp\) 的组合意义；下文中还会有更多类似的例子。
下文中需要对某些表性相似的元素加以合并，因此在这里引入等价关系和等价类的定义。
对于组合类 \(\mathcal A\)，设其中任意对象都是 \(\mathcal B\) 中某些对象的组合。我们定义 \(a\in \mathcal A\) 对 \(\mathcal B\) 的拆分为 \(a = (b_1, b_2, \dots, b_m), \forall b_i \in \mathcal B\)，这拆分一般应当是唯一的。记 \(a\) 在 \(\mathcal B\) 上的容为 \(\text{cap}_{\mathcal B}(a) = m\)，即拆分中的元素个数。
我们定义一个等价关系是一个置换群列，用加粗大写字母（\textbf）表示。假设 \(\textbf G\) 是一个置换群列，则可以记 \(\textbf G = \{G_0, G_1, G_2,\dots\}\)，其中 \(G_k\) 是一个置换群，满足任意 \(g\in G_k\) 都是一个大小为 \(k\) 的元素上的置换。若一个等价关系 \(\textbf G\) 是针对组合类 \(\mathcal A\) 中元素对 \(\mathcal B\) 的拆分定义的，则称 \(\textbf G\) 以 \(\mathcal B\) 为拆分集，记作 \(\textbf G_{\mathcal B}\)。
我们定义 \(a_1, a_2\) 对 \(\mathcal B\) 的拆分在 \(\textbf G\) 上等价（本质相同）当且仅当
- \(\text{cap}_{\mathcal B}(a_1) = \text{cap}_{\mathcal B}(a_2)\)；
- \(\exists g\in G_{\text{cap}_{\mathcal B}(a_1)}, g(a_1) = a_2\)。
注意这里 \(g\) 的作用是将 \(a_1\) 的拆分重排。
下面常将两个拆分 \((a_1, a_2, \dots a_n)\) 和 \((b_1, b_2, \dots, b_n)\) 在 \(\textbf G\) 上等价记作 \((a_1, a_2, \dots a_n)\textbf G(b_1, b_2, \dots, b_n)\)。下文中也会以下标变换的方式说明对应的 \(g\in G_n\) 的形式。
对于组合类 \(\mathcal A\) 和等价关系 \(\textbf G\)，取一个 \(\mathcal A\) 的子集 \(\mathcal S\)，若 \(\mathcal S\) 中任意两个对象等价，且任意 \(\mathcal S\) 中对象和 \(\complement_{\mathcal A}\mathcal S\) 中对象不等价，则我们称 \(\mathcal S\) 是 \(\textbf G\) 下的等价类。常记 \(\lvert \mathcal S\rvert\) 为等价类中一个元素的大小。需要注意的是，一般等价类中任意组合对象大小相等。
对于组合类 \(\mathcal A\) 和等价关系 \(\textbf G_{\mathcal B}\)，定义 \(\mathcal A / \textbf G_{\mathcal B}\) 得到了一个新的组合类，满足任意原等价类在新的组合类中有且仅有一个元素作为代表；这里常取在某一性质（例如字典序）下最小元素。有时也定义新的组合类中的对象是原组合类中的等价类，即新组合类是“组合类的组合类”。
在下面对组合类的构造过程中，新组合类的对象常对原组合类进行拆分，因此在语义明显时不声明原组合类为等价关系的拆分集。
\(1.2\) 基础无标号类
下面引入两种特殊的组合对象和对应的组合类。
我们记 \(\epsilon\) 为中性对象（neutral object），对应的组合类记作 \(\mathcal E = \{\epsilon \}\)，称作中性类（neutral class）。
恒有 \(\lvert \epsilon \rvert = 0\)，因此有中性类的 \(\text{OGF}\) \(E(z) = 1\)。
需要注意的是，两个中性对象 \(\epsilon_1, \epsilon_2\) 可能不同。记一个下标为 \(k\) 的中性对象 \(\epsilon_k\) 对应的类为 \(\mathcal E_k\)。
我们有
对于任意组合类 \(\mathcal A\)，规定 \(\mathcal A^0 = \mathcal E\)。但是由于组合意义，\(\mathcal E\) 不能视作笛卡尔积的单位元。
定义了中性对象后，我们就可以定义不需要考虑交集的集合并了。由于集合论中需要满足集合内无重合元素，我们可以对两个类分别乘上不同的中性类，这样两个类中可能相同的对象在并集中也彼此不同。
我们定义两个类 \(\mathcal A, \mathcal B\) 的集合并 \(\mathcal A + \mathcal B\) 为
其中 \(\cup\) 为集合的不交并。
我们记 \(\bullet\) 为原子对象（atom object），对应的组合类记作 \(\mathcal Z = \{\bullet \}\)，称作原子类（atom class）。
恒有 \(\lvert \bullet \rvert = 1\)，因此有原子类的 \(\text{OGF}\) \(Z(z) = z\)。
原子对象常用于合并数个组合对象，经典例子是作为树的根出现。
\(1.3\) 例题
\(\textbf{例 1}\)：01 串
我们如何求得长度为 \(n\) 的 01 串的个数？
设 01 串的组合类为 \(\mathcal S\)，我们要求的就是 \(\lvert\mathcal S_n\rvert\)。可以写出
对于某个 01 串，其要么为空串，要么是 \(\mathtt{0},\mathtt{1}\) 接在另一个 01 串前面。可以写出
翻译成生成函数就是 \(S(z) = 1 + (z + z) S(z)\)，也就能得到 \(S(z) = \dfrac{1}{1 - 2z}\)。答案即为 \([z^n]S(z) = 2^n\)。
\(\textbf{例 2}\)：卡特兰数
我们如何构造长度为 \(n\) 的生成函数？
卡特兰数有一个很经典的组合意义：\(Cat(n)\) 是 \(n\) 个点的二叉树个数。
对于一棵二叉树，它要么为空，要么形如 \(左子树-根-右子树\)。根可以视作 \(\bullet\)，两侧则又是二叉树。
假设二叉树的组合类为 \(\mathcal C\)，可以写出
翻译成生成函数就是 \(C(z) = 1 + zC(z)^2\)，这和我们先前得到的形式相同。
\(2.\) 经典的无标号构造
\(2.1\) 集合的 \(\text{Sequence}\) 构造
\(\text{Sequence}\) 构造生成了所有可能的有序组合。
我们定义
且要求 \(\mathcal A_0 = \varnothing\)，也就是 \(\mathcal A\) 中没有大小为 \(0\) 的对象。我们还可以将这构造生成的组合类写作
这映射翻译为生成函数即为
其中 \(Q\) 为 Pólya 准逆（quasi-inversion）。
\(\textbf{例 1}\)：写出 \(\{\text{a}, \text{b}\}\) 的 \(\text{Sequence}\) 构造：
\(\textbf{例 2}\)：有序有根树计数
我们可以使用 \(\text{Sequence}\) 构造来定义有序有根树（不同儿子之间的顺序有意义的有根树）的计数。
设对应的组合类是 \(\mathcal T\)，我们可以用一个原子对象和自身的 \(\text{Sequence}\) 描述自身，也就是
翻译为生成函数就是
从这里也可以看出，\(T(z) = \dfrac{1 - \sqrt{1 - 4z}}{2}\)，因此大小为 \(n\) 的有序有根树的大小为移位的卡特兰数。
\(2.2\) 集合的 \(\text{Multiset}\) 构造
\(\text{Multiset}\) 构造生成了所有可能的组合，但是不区分对象内部的元素的顺序。
我们仍然可以使用 \(\text{Sequence}\) 构造描述 \(\text{Multiset}\) 构造，但是由于顺序原因，我们可以递推解决，每次只拿出一个元素作 \(\text{Sequence}\) 构造。假设 \(\mathcal A = \left\{\alpha_0, \alpha_1, \dots, \alpha_k \right\}\)，则可以递推地作 \(\text{MSET}\) 构造如下：
即
且要求 \(\mathcal A_0 = \varnothing\)。可以得到等价的
其中 \((a_1, a_2, \dots, a_n)\textbf{R}(b_1, b_2, \dots, b_n)\) 当且仅当存在一个置换 \(\sigma\)，对于任意 \(j\) 满足 \(b_j = a_{\sigma(j)}\)。
翻译为生成函数就是
作 \(\ln - \exp\) 可以得到
其中 \(\text{Exp}\) 被称作为 Pólya 指数，又称 Euler 变换。
\(\textbf{例 1}\)：整数的分拆
对每个 \(1\le i \le 10^5\) 求 \(f(i)\)，其中 \(f(k)\) 是将 \(k\) 进行分拆的方案数。
设全体正整数类为 \(\mathcal I\)，可以发现
\(\text{SEQ}_{\ge 1}\) 是有限制的构造，在下方讲到。当然求上面的东西对答案没啥帮助，我们可以直接构造 \(I(z) = \sum_{i\ge 1} z^i\)。
所求的就是 \(\text{MSET}(\mathcal I)\)。
\(\textbf{例 2}\)：无标号无根树计数
我们能否对无标号无根树的计数构造生成函数？
假设无标号有根树的组合类是 \(\mathcal T\)，则可以写出
也就是
对于无根的情况，论文 The Number of Trees, Richard Otter 说明了无根树的 \(\text{OGF}\) 是
这点也可以通过对奇偶次项系数进行不同的讨论得到。
这个数列即为 A55。
\(\textbf{例 3}\)：Pólya 指数的逆变换
我们能否根据 \(\text{Exp}\) 的结构得到它的逆变换（不考虑常数项）？
记给定幂级数为 \(F(z)\)。我们需要求的就是一个序列 \(A\)，满足
两边取对数。得到
记 \(\ln F(z)\) 为 \(f\)，可以发现
作莫比乌斯反演得到
这就是逆运算的构造。
\(2.3\) 集合的 \(\text{Powerset}\) 构造
\(\text{Powerset}\) 构造生成了所有可能的子集。
可以发现，对于每个元素我们都可以选或不选，这会导出两类不同的子集。因此我们可以递归地定义 \(\text{PSET}\) 构造：
可以写作
要求 \(\mathcal A_0 = \varnothing\)。
容易发现 \(\text{PSET}(\mathcal A)\subset \text{MSET}(\mathcal A)\)，因为 \(\mathcal E + \{\alpha\} = \text{SEQ}_{k\le 1}(\{\alpha\})\)。
翻译为生成函数就是
\(2.4\) 集合的 \(\text{Amplification}\) 构造
\(\text{Amplification}\) 构造生成了所有可能的并列。
我们可以如此定义 \(\text{AMP}\) 构造：
要求 \(\mathcal A_0 = \varnothing\)。
也就是说，任意多个 \(\mathcal A\) 中的元素 \(\alpha\) 自身的并列组成了 \(\text{AMP}(\mathcal A)\)。
翻译成生成函数就是
\(2.5\) 集合的 \(\text{Cycle}\) 构造
\(\text{Cycle}\) 构造生成了所有可能的组合，但是不区分仅轮换不同的组合。
我们定义这个构造为
其中 \(\textbf S\) 为等价关系，\((a_1, a_2, \dots a_n) \textbf S(b_1, b_2, \dots, b_n)\) 当且仅当存在循环移位 \(\tau\) 使得对于任意 \(j\) 有 \(b_j = a_{\tau(j)}\)。
翻译为生成函数是
\(\text{Log}\) 被称作为 Pólya 对数。公式的证明可以参见 The Cycle Construction, P. Flajolet and M. Soria，这里不再展开。采用群论相关知识我们也可以获得相同的结果。
\(\textbf{例 1}\)：列举 \(\text{CYC}(\{\text{a}, \text{b}\})\) 中长度为 \(4\) 的对象。
它们分别为
由于 \(\text{aaab}\textbf S\text{aaba}\textbf S\text{abaa}\textbf S\text{baaa}, \text{aabb}\textbf S\text{abba}\textbf S\text{bbaa}\textbf S\text{baab}, \text{abbb}\textbf S\text{bbba}\textbf S\text{bbab}\textbf S\text{babb}, \text{abab}\textbf S\text{baba}\)，这几个等价类都只有一个元素出现在 \(\text{CYC}(\{\text{a}, \text{b}\})\) 中。
这里采用字典序最小的字符串作为代表。
\(2.6\) 有限制的构造
对于上述构造，我们并没有限制每个对象组成部分的个数。
这里以 \(\text{SEQ}\) 构造为例。我们若在 \(\text{SEQ}\) 的下标加上作用于整数的谓词用于约束其组成部分的个数，如 \(\text{SEQ}_{= k}, \text{SEQ}_{\ge k}, \text{SEQ}_{\in [1, k]}\)，则表示构造出的等价类中每个对象的拆分中相同元素需要被下标上对应的谓词限制。常简写 \(\text{SEQ}_{= k}\) 为 \(\text{SEQ}_{k}\)。
令一个构造 \(\text{CONS}\) 为上述构造中的一种，并设 \(\mathcal A = \text{CONS}_k(\mathcal B)\)，则我们需要对 \(\forall \alpha \in \mathcal A\) 有 \(\alpha = \left\{(\beta_1, \beta_2, \dots, \beta_k) \mid \forall \beta_i \in \mathcal B \right\}\)。
这种构造我们在先前已经充分接触过了，其翻译到生成函数上就是用新占位元来标识组成部分的个数。
组合意义上，我们定义 \(\chi\) 函数作用于一个元素上标识其组成部分的个数，也就是它需要被对应的逻辑谓词限制。用在上面的例子上就是 \(\chi(\alpha) = k\)。
延续先前的记号，令
翻译为生成函数即为
我们只需要提取 \(t^k\) 项系数即可得到对应的表达式。对应 \(\text{SEQ}_k\) 就能得到
也就是 \(A(z) = B(z)^k\)。同样能得到 \(\text{SEQ}_{\ge k}(\mathcal B)\) 能引出 \(A(z) = \dfrac{B(z)^k}{1 - B(z)}\)。
对于 \(\text{MSET}_k\) 有
这个得提取系数，得到
对于有限制的构造，没有要求 \(\mathcal B_0 = \varnothing\)。
常用的有限制构造
我们或许也可以通过 Pólya 定理导出结论。讲解可以看 HandWiki、Ency 或 x义x.
\(\textbf{例 1}\)：烷基计数
计数 \(n\) 个节点的根节点度数不超过 \(3\)，其余节点度数不超过 \(4\) 的无序有根树。
考虑我们把更小的满足条件的树的根连接到新根节点上时，会新建一条边，这也不会让原来的根的度数超过 \(4\)。因此假设满足条件的树的组合类为 \(\mathcal T\)，我们可以写出
或者直接按能接上一棵空树来计数。假设 \(\hat{\mathcal T} = \mathcal T + \mathcal E\)，我们可以写出
\(\textbf{例 2}\)：有序树的进一步探讨
承接 \(2.1\) 例 \(2\)，我们限制每个节点子树的数量。
令 \(\Omega\subseteq \mathbb N\) 为一个包含 \(0\) 的集合。我们定义组合类 \(\mathcal T^{\Omega}\)，其中的对象被称作 \(\Omega\) 树。一棵 \(\Omega\) 树需要满足任意节点的儿子数量 \(\omega \in \Omega\)。在接下来的讨论中，常要用到 \(\Omega\) 的特征函数
举几个例子。\(\Omega = \{0, 2\}\) 对应的组合对象是满二叉树，其中每个节点要么没有儿子，要么有两个儿子。它的特征函数是 \(\phi(t) = 1 + t^2\)。\(\Omega = \{0, 1, 2\}\) 对应的组合对象是一般的二叉树，特征函数是 \(\phi(t) = 1 + t + t^2\)；\(\Omega = \{0, 3\}\) 对应的组合对象是三叉树，特征函数是 \(\phi(t) = 1 + t^3\)。对于一般的树有 \(\Omega = \mathbb N\)，特征函数是 \(\phi(t) = (1 - t)^{-1}\)。
我们可以通过 \(\Omega\) 的特征函数 \(\phi\) 递归定义 \(\mathcal T^{\Omega}\) 的 \(\text{OGF}\) \(T^{\Omega}(z)\)：
这形式的证明需要从 \(\text{SEQ}\) 构造着手。一个节点限制儿子的数量后子树的选择方案可以施 \(\text{SEQ}_{\in \Omega}\) 构造，对组合类我们能写出 \(\mathcal T^{\Omega} = \mathcal Z\times \text{SEQ}_{\in \Omega} (\mathcal T^{\Omega})\)。之前的讨论可以推广到这里，这种构造翻译成生成函数就是复合 \(\phi\)，也就是 \(T^{\Omega}(z) = z\times \phi\left(T^{\Omega}(z)\right)\)。
上面的形式启发我们通过拉格朗日反演提取系数。我们知道 \(\dfrac{T^{\Omega}(z)}{\phi\left(T^{\Omega}(z)\right)} = z\)，这就有 \(\left(T^{\Omega}(z)\right)^{\langle -1\rangle} = \dfrac{z}{\phi(z)}\)。可以得到
如上的思路在 \(\Omega\) 为多重集时自然成立。例如，\(\Omega = \{0, 1, 1, 3\}\) 给出了一种一叉-三叉树，这类树可以有两种一叉节点，分别染了不同的颜色；这集合的特征函数就是 \(\phi(z) = 1 + 2z + z^3\)。
由上面的讨论不难看出，如果 \(\Omega\) 包含了 \(r\) 个不同的元素，则对 \(\Omega\) 树的计数可以用一个 \(r - 1\) 重的二项式系数求和表示。
我们也可以用这种方法计数 \(k\) 阶 \(\Omega\) 树森林，其组合类定义为 \(\mathcal F = \text{SEQ}_k(\mathcal T^{\Omega})\)。关注生成函数能得到 \(F(z) = T^{\Omega}(z)^k\)。
同样应用拉格朗日反演可以得到
当 \(\Omega = \mathbb N\) 时，我们有 \(k\) 阶一般树森林的计数为
系数对应的数列又被称为 ballot numbers。
当 \(\Omega = \{0, 1, 2\}\) 时，令 \(\mathcal M\) 为对应的一般二叉树的组合类，不难写出 \(\mathcal M = \mathcal Z\times \text{SEQ}_{\le 2}(\mathcal M)\)，也就能得到
可以用拉格朗日反演导出
系数对应的数列（移位）又被称为 Motzkin numbers。
可能的例题：P7592 数树
\(\textbf{例 3}\)：被括号包含（bracketed）
这个问题被论文 Vier combinatorische Probleme, E. Schröder 系统地阐述后广为人知，它也是这篇论文中的第二个问题。
我们需要计数的是 \(n\) 个相同的字母 \(x\) 有多少种被括号包含的合法方案。合法方案是递归地定义的：
- \(x\) 本身是合法方案。
- 若 \(\sigma_1, \sigma_2, \dots, \sigma_k\) 分别是合法方案，且 \(k \ge 2\)，则 \((\sigma_1 \sigma_2 \cdots \sigma_k)\) 是合法的，也就是并列后在最外层加一对括号。
例如，\((((xx)x(xxx))((xx)(xx)x))\) 是合法的，\(((((x)(x))))\) 就不是。
我们令 \(\mathcal S\) 表示合法方案对应的组合类，大小函数表示一个方案中 \(x\) 的个数。上面的递归定义可以被形式化地表示为
翻译得到
考虑一个合法序列和一棵树的关系。可以发现，每个 \(x\) 代表的就是这棵树的叶子节点，而一对括号代表了树的内部节点，我们需要保证内部节点的度数 \(\ge 2\)。叶子节点在这里对大小函数有贡献，内部节点则没有。
注意这里树的形式和例 \(2\) 中描述的形式不同，大小函数的定义也不同，无法直接应用上面的构造。但是拉格朗日反演仍然可以应用，这里有
提取 \(z^n\) 项系数可以得到
这个其实微分有限，化成组合数形式也只是为了得到新的组合意义。
系数对应的数列是 A1003。
\(2.7\) 隐式构造
在很多情况下，我们会见到隐式构造的组合类 \(\mathcal X\)，它作为一个映射的原像出现，而像已经被定义了。我们希望刻画这种组合类，并将他显式地放在像的位置。这就促使我们在组合类运算中定义逆。
我们举几个例子，其中 \(\mathcal A, \mathcal B\) 已知，\(\mathcal X\) 是新定义的：
这里用生成函数的语言更容易定义逆。分别地，我们有
\(\text{MSET}\) 的逆的形式已经在 \(2.2\) 例 \(3\) 讨论了。
通过组合类间关系隐式构造的组合类被称为隐式结构（implicit structures）。
\(\textbf{例 1}\)：排列的连续段分解
一个排列 \(\sigma = \sigma_1 \sigma_2 \cdots \sigma_n\) 是可分解的，当且仅当其下标可以被分割为大于一个不交的区间 \([l_i, r_i] = \{k \mid k\in \mathbb N \land l_i\le k \le r_i \}\)，满足
- 这些区间的并集是 \([1,n] = \{k \mid k\in \mathbb N \land 1\le k \le n \}\)；
- 对任何一个区间 \([l,r]\)，\(\forall i \in [l, r], \exist j\in [l, r], \sigma_i = j\)；
- 对任意两个区间 \([l_i, r_i], [l_j, r_j]\) 满足 \(i < j\)，有 \(\forall x \in [l_i, r_i], \forall y \in [l_j, r_j], \sigma_x < \sigma_y\)；
- 区间数量极大。
举例来说，\(\sigma = \{2,5,4,1,3,6,8,10,7,9\}\) 的分解就是 \([1, 5], [6, 6], [7, 10]\)，用图像表示就是
我们称无法执行上述操作的排列是不可分解的排列，例如 \(\{2,4,1,3\}, \{5,4,3,2,1\}\)。现在请构造不可分解的排列的组合类 \(\mathcal I\)。
我们取出排列 \(\sigma\) 的分解中的一段 \([l,r]\) 来观察。不难发现，排列 \(\{\sigma_l - l + 1, \sigma_{l+1} - l + 1, \dots, \sigma_r - l + 1 \}\) 是一个新的 \(1\sim r - l + 1\) 的排列，且这个排列是不可再分解的。我们如果只关注排列的大小，则一个可分解的排列是一系列不可分解的排列的并列。
假设任意排列的组合类是 \(\mathcal P\)，则我们能得到
这隐式定义了 \(\mathcal I\)，可以写出
系数对应的数列是 A3319。
\(\textbf{例 2}\)：系数在有限域内的素多项式
我们固定一个质数 \(p\)，下面所有系数间的运算都在模 \(p\) 数域 \(\mathbb F_p\) 内运算。记多项式环 \(\mathbb F_p[x]\) 是系数在 \(\mathbb F_p\) 内的多项式的全体组成的环。
我们在这里讨论的是首一多项式，即最高次项系数为 \(1\) 的多项式。记 \(\mathbb F_p[x]\) 集合内的首一多项式组成了组合类 \(\mathcal P\)，大小函数定义为多项式最高次项。记 \(\mathcal A = \mathbb F_p\)。由于一个多项式可以由它的系数唯一确定，我们可以将它视作由 \(\mathcal A\) 内对象组成的序列，因此有
这也能得到有 \(p^n\) 个最高次项为 \(x^n\) 的首一多项式。
众所周知，多项式域 \(\mathbb F_p[x]\) 上可以作欧几里得法求最大公因子，因此多项式域是 \(\text{UFD}\)。我们称没有除常数外的因子的多项式是素多项式，是不可分解的。举例来说，在 \(\mathbb F_3[x]\) 上有
我们令 \(\mathcal I\) 为素多项式对应的组合类。唯一分解性质说明，\(\mathcal P\) 在组合意义上和 \(\mathcal I\) 的 \(\text{Multiset}\) 构造同构，也就是
这样我们就能写出
我们也自然能通过上面的方法以及提取系数得到
特别的，\(\mathcal I_n = O(p^n / n)\)。这个结论已经被高斯知晓。因此 \(n\) 度的多项式中任选一个多项式是素多项式的概率为 \(1/n\)。
\(3.\) 有标号体系
\(3.1\) 概述
有标号体系是无标号体系的自然拓展，大部分组合对象的性质都可以自然地引入有标号体系中。
有标号体系和无标号体系最大的区别在于，有标号体系所关注的对象都有着两两不同的标号，彼此可以区分。不失一般性地，我们将标号集合视作正整数集合。举例来说，一个排列可以被视作一系列两两不同正整数依次排开，而循环分解将这排列视作循环有向图的组合，图的每个顶点都有着正整数标号。
有标号体系内的操作基于一类特殊的乘法：有标号乘（labelled product）。它是对无标号体系中笛卡尔积的自然模拟。
同样的，我们也需要得到对无标号体系中 \(\text{OGF}\) 的有标号模拟，这就是指数生成函数，简称 \(\text{EGF}\)。在下面的叙述中我们可以看到，\(\text{EGF}\) 可以很好地描述计数序列。
新的构造使得我们能关注更多的结构，特别是在有顺序的性质上。
\(3.2\) 记号与定义
我们所处理的对象和无标号体系一样，是有限的组合类与组合对象。不同的是，我们现在处理的组合对象是有标签的。每个组合对象都[有一种独特的颜色/被分配到一个正整数标号]，以保证组合对象是两两不同的。
一个大小为 \(n\) 的弱标号对象是一种有 \(n\) 个部分的离散结构，其 \(n\) 个部分分别有一个标号，这些标号组成了一个大小为 \(n\) 的整数集。我们称满足这条件的对象是弱标号的。等价地，可以说每个部分都有标号，这种叙述隐含了标号是 \(\in \mathbb Z\) 的彼此不同的数。一个大小为 \(n\) 的有标号（又称强标号，一般在需要区分弱标号对象的语境下使用）对象需要是弱标号的，且其标号组成的集合是整数区间 \([1, n] = \{x\mid x\in \mathbb N, 1\le x\le n\}\)。一个例子是一张有 \(n\) 个节点的图，它的一个部分是一个带标号节点。
一个有标号类是一个由有标号的组合对象组成的组合类，这里把 \(1.\) 中介绍过的无标号的组合类称作无标号类，以示区分。对有标号类的记号仍然沿用美术花体（\mathcal）大写字母，其中的一个有标号对象常记作对应的小写字母/希腊字母，对应的生成函数用对应的一般字体大写字母记。有标号类也是一个组合类，大小函数、计数序列和同构等定义与 \(1.\) 中介绍的无标号类相同。
一个有标号类 \(\mathcal A\) 的 \(\text{EGF}\) 可以用它的计数序列 \(A[n]\) 或其中元素表示，写作
这里的 \(z\) 是用于标记组合对象大小的占位元。
这里作一下关于 \(\text{EGF}\) 的提取系数的声明。提取系数需要满足 \(A[n] = n!\times [z^n]A(z)\)，也就是 \([z^n]A(z) = A[n]/n!\)。
\(3.3\) 基础有标号类与操作
在无标号体系中，我们已经见到了中性类和原子类发挥的作用，不妨在有标号体系中也定义这样的结构。
我们记 \(\epsilon\) 为中性对象（neutral object），对应的组合类记作 \(\mathcal E = \{\epsilon \}\)，称作中性类（neutral class）。
恒有 \(\lvert \epsilon \rvert = 0\)，因此有中性类的 \(\text{EGF}\) \(E(z) = 1\)。
我们记 \(①\) 为原子对象（atom object），对应的组合类记作 \(\mathcal Z = \{① \}\)，称作原子类（atom class）。我们令原子对象中唯一的元素标号恒为 \(1\)，这使得原子对象是强标号的。
恒有 \(\lvert① \rvert = 1\)，因此有原子类的 \(\text{EGF}\) \(Z(z) = z\)。
有标号类的集合并可以类似无标号类的集合并一样定义，这里也记作 \(+\)。\(\mathcal A + \mathcal B\) 翻译到生成函数上就是 \(A(z) + B(z)\)。
而有标号类的乘法不好得到良定义，这是由于笛卡尔积的结构可能导出非良标号对象，也就是两个有标号对象可能存在标号重叠的情况。我们不采用笛卡尔积，而是定义新的运算：有标号乘法，这运算可以自然地翻译为 \(\text{EGF}\) 的乘法。记 \(\mathcal A\) 与 \(\mathcal B\) 的有标号乘法为 \(\mathcal A \times \mathcal B\)，它翻译到生成函数上是 \(A(z)\times B(z)\)。注意这里没有定义 \(\mathcal A,\mathcal B\) 的笛卡尔积，我们也不使用 \(\star\) 代替 \(\times\)。
首先介绍二项式卷积，也就是 \(\text{EGF}\) 在做乘法运算时系数的映射结构。假设 \(F(z), G(z), H(z)\) 分别是序列 \(f, g, h\) 的 \(\text{EGF}\)，且 \(F(z) = G(z)\times H(z)\)，能发现它们的系数满足
同样的，如果我们有 \(F(z) = G_1(z)\times G_2(z)\times\cdots\times G_r(z)\)，则
这里的系数是多项式系数，有
这种结构对我们下面的讨论很有帮助，是 \(\text{EGF}\) 应用的核心所在。
我们接下来要讨论的是重标号操作。由于我们需要让两个强标号对象合并后仍然是强标号的，我们需要进行重标号操作，这种操作要求新的标号的大小关系中保留原标号之间的大小关系。下面列出的两种重标号方式是重要的：
- 减缩。
对于一个大小为 \(n\) 的弱标号结构，减缩操作将每个点的标号映射到整数区间 \([1, n]\) 内，保留了原标号的大小关系，相当于作了离散化。我们记对结构 \(\alpha\) 做减缩操作得到 \(\rho(\alpha)\)。
举例来说，序列 \(\langle 7, 3, 9, 2 \rangle\) 减缩得到 \(\langle 3, 2, 4, 1 \rangle\)。 - 增扩。
定义该操作需要一个增扩函数 \(e : [1, n]\to \mathbb Z\)，满足 \(\forall i\in [1, n], e(i) > i\)。对于一个强标号对象 \(\alpha\)，定义它增扩后得到一个弱标号对象 \(\hat\alpha\)，\(\alpha\) 中标号为 \(k\) 的部分在 \(\hat\alpha\) 中标号为 \(e(k)\)。我们记对结构 \(\alpha\) 做增扩操作得到 \(e(\alpha)\)。
举例来说，序列 \(\langle 3, 2, 4, 1 \rangle\) 增扩可以得到 \(\langle 33, 22, 44, 11 \rangle\) 或 \(\langle 7, 3, 9, 2 \rangle\)。
我们可以通过这些操作设计出一种适于有标号类的乘法，其最初由 D. Foata 提出，取名为 partitional product。有标号乘法的思路就是对合并后的部分进行重标号，以避免出现同一个标号出现多次。
我们取两个有标号对象 \(\beta \in \mathcal B\) 和 \(\gamma \in \mathcal C\)，我们定义这两个对象的有标号乘法（简称乘法）记作 \(\beta \times \gamma\)，定义为一个集合，其包含所有由 \((\beta, \gamma)\) 生成的强标号的有序对 \((\beta', \gamma ')\)。具体地，我们定义
这个定义同样可以通过增扩导出，不再赘述。
注意到有标号乘积导出的构造中每个对象都是强标号的。假设 \(\beta, \gamma\) 都是强标号的，大小分别为 \(n_1, n_2\)，则分配标号的方案总共有 \(\dbinom{n_1 + n_2}{n_1}\) 种，也就是说
随后我们就能对有标号类自然地定义有标号乘法了。
设 \(\mathcal B, \mathcal C\) 是两个有标号类，则 \(\mathcal B\) 和 \(\mathcal C\) 的有标号乘法记作 \(\mathcal B \times \mathcal C\)，定义为
随后我们可以自然地描述计数序列的关系了。假设 \(\mathcal A = \mathcal B \times \mathcal C\)，则我们能知道
可以发现这转移和 \(\text{EGF}\) 的乘法同构。因此我们构造出了一种有意义的有标号乘法，使得 \(\text{EGF}\) 可以用来描述组合类的计数序列。
通过这个构造，我们就可以构造有标号集合、有标号序列和有标号循环了，方法和无标号的构造法类似。在这过程中我们也会注意可容许性。
\(4.\) 经典的有标号构造
\(4.1\) 有标号 \(\text{Sequence}\) 构造
有标号 \(\text{Sequence}\) 构造生成了所有可能的有标号序列。
一个有标号类 \(\mathcal B\) 的 \(k\) 次（有标号）幂被定义为 \(\mathcal A = \mathcal B\times \mathcal B\times \cdots \times \mathcal B\)，这里有 \(k\) 个 \(\mathcal B\) 作乘法。翻译为生成函数为 \(A(z) = B(z)^k\)。
\(\mathcal B\) 的 \(k\) 次（有标号）幂记作 \(\text{SEQ}_k(\mathcal B)\)。
我们定义
记 \(\mathcal B = \text{SEQ}(\mathcal A)\)。翻译为生成函数即为
这里仍然需要 \(\mathcal A_0 \neq \varnothing\)。
\(4.2\) 有标号 \(\text{Set}\) 构造
有标号 \(\text{Set}\) 构造生成了所有可能的有标号序列，但是不区分组成对象的元素的顺序。
由有标号类 \(\mathcal B\) 中对象作为元素的大小为 \(k\) 的集合的全体就是有标号类 \(\mathcal B\) 的 \(\text{Set}_k\) 构造。它是无标号类中 \(\text{Multiset}, \text{Powerset}\) 构造的模拟，定义为 \(\text{SET}_k(\mathcal B) = \text{SEQ}_k(\mathcal B) / \textbf R\)，其中 \((a_1, a_2, \dots, a_n)\textbf{R}(b_1, b_2, \dots, b_n)\) 当且仅当存在一个置换 \(\sigma\)，对于任意 \(j\) 满足 \(b_j = a_{\sigma(j)}\)。
由有标号类 \(\mathcal B\) 中对象作为元素的集合的全体就是有标号类 \(\mathcal B\) 的 \(\text{Set}\) 构造。我们定义
对于一个有标号类，有标号 \(\text{Set}_k\) 构造中的一个元素对应着有标号 \(\text{Sequence}\) 构造中的 \(k!\) 个元素，这可以通过计数置换得到。因此可以将这构造翻译为生成函数：
记 \(\mathcal B = \text{SET}_k(\mathcal A)\)。\(\text{Set}_k\) 构造翻译为生成函数即为
记 \(\mathcal B = \text{SET}(\mathcal A)\)。\(\text{Set}\) 构造翻译为生成函数即为
这里仍然需要 \(\mathcal A_0 \neq \varnothing\)。
注意到无标号 \(\text{Powerset}\) 构造和无标号 \(\text{Multiset}\) 构造在带标号后同构，\(\text{Powerset}\) 构造可以并入，统称为有标号 \(\text{Set}\) 构造。
\(4.3\) 有标号 \(\text{Cycle}\) 构造
\(\text{Cycle}\) 构造生成了所有可能的有标号序列，但是不区分仅轮换不同的序列。
有标号 \(\text{Cycle}\) 构造是无标号类中 \(\text{Cycle}\) 构造的模拟，定义为 \(\text{CYC}_k(\mathcal B) = \text{SEQ}_k(\mathcal B) / \textbf S\)，其中 \((a_1, a_2, \dots a_n) \textbf S(b_1, b_2, \dots, b_n)\) 当且仅当存在循环移位 \(\tau\) 使得对于任意 \(j\) 有 \(b_j = a_{\tau(j)}\)。
同上地定义 \(\text{CYC}(\mathcal B) = \sum_{k\ge 0}\text{CYC}_k(\mathcal B)\)。
可以发现，\(\text{CYC}_k(\mathcal B)\) 中任意一个对象在 \(\text{SEQ}_k(\mathcal B)\) 中出现了恰好 \(k\) 次，因为每个对象都可以旋转 \(k\) 次。因此可以将这构造翻译为生成函数：
记 \(\mathcal B = \text{CYC}_k(\mathcal A)\)。\(\text{Cycle}_k\) 构造翻译为生成函数即为
记 \(\mathcal B = \text{CYC}(\mathcal A)\)。\(\text{Cycle}\) 构造翻译为生成函数即为
这里仍然需要 \(\mathcal A_0 \neq \varnothing\)。
\(4.4\) 有标号 \(\text{Pointing}\) 构造
\(\text{Pointing}\) 构造生成了所有可能的“有根”结构。
对于一个有标号类 \(\mathcal A\)，\(\text{Pointing}\) 构造生成了一个有标号类 \(\mathcal B\)，使得对于任意 \(\mathcal A_n\)，我们都用 \(n\) 个彼此不同的中性对象分别单独标明 \(\forall \alpha \in \mathcal A_n\) 中一个单位大小的元素，即对于 \(\alpha\) 对应的 \(\{\beta\} \subseteq \mathcal B\)，有 \(\{\beta\} = \{\alpha \}\times \{\epsilon_1, \epsilon_2, \dots, \epsilon_n\}\)。容易发现 \(\forall \alpha \in \mathcal A_n\) 都生成了 \(n\) 个不同的 \(\beta \in \mathcal B\)；当 \(\mathcal A\) 是有标号无根树类时，\(\mathcal B\) 就是有标号有根树类。
\(\text{PNT}(\mathcal A)\) 又记作 \(\Theta \mathcal A\)。
容易翻译为生成函数：若 \(\mathcal A = \text{PNT}(\mathcal B)\) 则
这就是 \(A = \vartheta B\)。
\(4.5\) 有标号 \(\text{Substitution}\) 构造
\(\text{Substitution}\) 构造生成了所有可能的替换生成结构。
\(\text{Substitution}\) 构造是二元运算，对有序二元组 \((\mathcal A, \mathcal B)\) 的 \(\text{Substitution}\) 构造常记作 \(\mathcal A \circ \mathcal B\)。\(\text{Substitution}\) 构造将前一个有标号类中的元素的每个单位元替换为后一个组合类中的元素。形式化地，定义 \(\mathcal S \boxtimes \mathcal T = \{(s, t) \mid s \in \mathcal S, t \in \mathcal T , |(s, t)| = |t|\}\)，则有
这翻译成生成函数就是复合：
我们也可以通过这构造来导出其他一些构造。
举例来说，设排列的有标号类是 \(\mathcal P\)，任意有标号类 \(\mathcal A\) 的 \(\text{Sequence}\) 构造就是 \(\mathcal P\circ \mathcal A\)。由于 \(P(x) = \dfrac{1}{1 - x}\)，我们可以得到 \(\text{Sequence}\) 构造的经典表述 \(\dfrac{1}{1 - A(x)}\)。
在多项式复合的复杂度被优化至 \(O(n\log^2 n)\) 后，该构造的适用性逐渐提升。
参考资料：
《Analytic Combinatorics》, Philippe Flajolet and Robert Sedgewick ;
oi-wiki 符号化方法, hly1204 et al. ;
多项式计数杂谈, command_block ;
组合结构符号化学习笔记, x义x ;
抄袭 x义x 的 Symbolic Method 讲义, alpha1022 ;
待补：增加有标号例题
以下是博客签名，与正文无关。
请按如下方式引用此页：
本文作者 joke3579，原文链接：https://www.cnblogs.com/joke3579/p/symbolic_method.html。
遵循 CC BY-NC-SA 4.0 协议。
请读者尽量不要在评论区发布与博客内文完全无关的评论，视情况可能删除。
