排列组合

引入

排列组合是组合数学中的基础.排列就是指从给定个数的元素中取出指定个数的元素进行排序;组合则是指从给定个数的元素中仅仅取出指定个数的元素,不考虑排序.排列组合的中心问题是研究给定要求的排列和组合可能出现的情况总数.排列组合与古典概率论关系密切.

在高中初等数学中,排列组合多是利用列表、枚举等方法解题.

加法 & 乘法原理

加法原理

完成一个工程可以有 nn 类办法,ai(1≤i≤n)a_i(1 \le i \le n) 代表第 ii 类方法的数目.那么完成这件事共有 S=a1+a2+⋯+anS=a_1+a_2+\cdots +a_n 种不同的方法.

乘法原理

完成一个工程需要分 nn 个步骤,ai(1≤i≤n)a_i(1 \le i \le n) 代表第 ii 个步骤的不同方法数目.那么完成这件事共有 S=a1×a2×⋯×anS = a_1 \times a_2 \times \cdots \times a_n 种不同的方法.

排列与组合基础

排列数

从 nn 个不同元素中,任取 mm(m≤nm\leq n,mm 与 nn 均为自然数,下同)个元素按照一定的顺序排成一列,叫做从 nn 个不同元素中取出 mm 个元素的一个排列;从 nn 个不同元素中取出 mm(m≤nm\leq n) 个元素的所有排列的个数,叫做从 nn 个不同元素中取出 mm 个元素的排列数,用符号 Anm\mathrm A_n^m(或者是 Pnm\mathrm P_n^m)表示.

排列的计算公式如下:

Anm=n(n−1)(n−2)⋯(n−m+1)=n!(n−m)!\mathrm A_n^m = n(n-1)(n-2) \cdots (n-m+1) = \frac{n!}{(n - m)!}

n!n! 代表 nn 的阶乘,即 6!=1×2×3×4×5×66! = 1 \times 2 \times 3 \times 4 \times 5 \times 6.

公式可以这样理解:nn 个人选 mm 个来排队 (m≤nm \le n).第一个位置可以选 nn 个,第二位置可以选 n−1n-1 个,以此类推,第 mm 个(最后一个)可以选 n−m+1n-m+1 个,得:

Anm=n(n−1)(n−2)⋯(n−m+1)=n!(n−m)!\mathrm A_n^m = n(n-1)(n-2) \cdots (n-m+1) = \frac{n!}{(n - m)!}

全排列:nn 个人全部来排队,队长为 nn.第一个位置可以选 nn 个,第二位置可以选 n−1n-1 个,以此类推得:

Ann=n(n−1)(n−2)⋯3×2×1=n!\mathrm A_n^n = n(n-1)(n-2) \cdots 3 \times 2 \times 1 = n!

全排列是排列数的一个特殊情况.

组合数

从 nn 个不同元素中,任取 m≤nm \leq n 个元素组成一个集合,叫做从 nn 个不同元素中取出 mm 个元素的一个组合;从 nn 个不同元素中取出 m≤nm \leq n 个元素的所有组合的个数,叫做从 nn 个不同元素中取出 mm 个元素的组合数,用符号 (nm)\dbinom{n}{m} 来表示,读作「nn 选 mm」.

组合数计算公式

(nm)=Anmm!=n!m!(n−m)!\dbinom{n}{m} = \frac{\mathrm A_n^m}{m!} = \frac{n!}{m!(n - m)!}

如何理解上述公式?我们考虑 nn 个人选 mm 个出来(m≤nm \le n),不排队,不在乎顺序.如果在乎顺序那么就是 Anm\mathrm A_n^m,如果不在乎那么就要除掉重复,那么重复了多少?同样选出来的 mm 个人,他们还要「全排」得 m!m!,所以得:

(nm)×m!=Anm(nm)=Anmm!=n!m!(n−m)!\begin{aligned} \dbinom{n}{m} \times m! &= \mathrm A_n^m\\ \dbinom{n}{m} &= \frac{\mathrm A_n^m}{m!} = \frac{n!}{m!(n-m)!} \end{aligned}

组合数也常用 Cnm\mathrm C_n^m 表示,即 Cnm=(nm)\displaystyle \mathrm C_n^m=\binom{n}{m}.现在数学界普遍采用 (nm)\dbinom{n}{m} 的记号而非 Cnm\mathrm C_n^m.

组合数也被称为「二项式系数」,下文二项式定理将会阐述其中的联系.

特别地,规定当 m>nm>n 时,Anm=(nm)=0\mathrm A_n^m=\dbinom{n}{m}=0.

插板法

插板法(Stars and bars)是用于求一类给相同元素分组的方案数的一种技巧,也可以用于求一类线性不定方程的解的组数.

正整数和的数目

问题一:现有 nn 个 完全相同 的元素,要求将其分为 kk 组,保证每组至少有一个元素,一共有多少种分法?

考虑拿 k−1k - 1 块板子插入到 nn 个元素两两形成的 n−1n - 1 个空里面.

因为元素是完全相同的,所以答案就是 (n−1k−1)\dbinom{n - 1}{k - 1}.

本质是求 x1+x2+⋯+xk=nx_1+x_2+\cdots+x_k=n 的正整数解的组数.

非负整数和的数目

问题二:如果问题变化一下,每组允许为空呢?

显然此时没法直接插板了,因为有可能出现很多块板子插到一个空里面的情况,非常不好计算.

我们考虑创造条件转化成有限制的问题一,先借 kk 个元素过来,在这 n+kn + k 个元素形成的 n+k−1n + k - 1 个空里面插板,答案为

(n+k−1k−1)=(n+k−1n)\binom{n + k - 1}{k - 1} = \binom{n + k - 1}{n}

虽然不是直接求的原问题,但这个式子就是原问题的答案,可以这么理解:

开头我们借来了 kk 个元素,用于保证每组至少有一个元素,插完板之后再把这 kk 个借来的元素从 kk 组里面拿走.因为元素是相同的,所以转化过的情况和转化前的情况可以一一对应,答案也就是相等的.

由此可以推导出插板法的公式:(n+k−1n)\dbinom{n + k - 1}{n}.

本质是求 x1+x2+⋯+xk=nx_1+x_2+\cdots+x_k=n 的非负整数解的组数(即要求 xi≥0x_i \ge 0).

不同下界整数和的数目

问题三:如果再扩展一步,要求对于第 ii 组,至少要分到 ai,∑ai≤na_i,\sum a_i \le n 个元素呢?

本质是求 x1+x2+⋯+xk=nx_1+x_2+\cdots+x_k=n 的解的数目,其中 xi≥aix_i \ge a_i.

类比无限制的情况,我们借 ∑ai\sum a_i 个元素过来,保证第 ii 组至少能分到 aia_i 个.也就是令

xi′=xi−aix_i^{\prime}=x_i-a_i

得到新方程:

(x1′+a1)+(x2′+a2)+⋯+(xk′+ak)=nx1′+x2′+⋯+xk′=n−a1−a2−⋯−akx1′+x2′+⋯+xk′=n−∑ai\begin{aligned} (x_1^{\prime}+a_1)+(x_2^{\prime}+a_2)+\cdots+(x_k^{\prime}+a_k)&=n\\ x_1^{\prime}+x_2^{\prime}+\cdots+x_k^{\prime}&=n-a_1-a_2-\cdots-a_k\\ x_1^{\prime}+x_2^{\prime}+\cdots+x_k^{\prime}&=n-\sum a_i \end{aligned}

其中

xi′≥0x_i^{\prime}\ge 0

然后问题三就转化成了问题二,直接用插板法公式得到答案为

(n−∑ai+k−1n−∑ai)\binom{n - \sum a_i + k - 1}{n - \sum a_i}

不相邻的排列

1∼n1 \sim n 这 nn 个自然数中选 kk 个,这 kk 个数中任何两个数都不相邻的组合有 (n−k+1k)\dbinom {n-k+1}{k} 种.

二项式定理

在进入排列组合进阶篇之前,我们先介绍一个与组合数密切相关的定理——二项式定理.

二项式定理阐明了一个展开式的系数:

(a+b)n=∑i=0n(ni)an−ibi(a+b)^n=\sum_{i=0}^n\binom{n}{i}a^{n-i}b^i

证明可以采用数学归纳法,利用 (nk)+(nk−1)=(n+1k)\dbinom{n}{k}+\dbinom{n}{k-1}=\dbinom{n+1}{k} 做归纳.

二项式定理也可以很容易扩展为多项式的形式:

设 nn 为正整数,xix_i 为实数,

(x1+x2+⋯+xt)n=∑满足n1+⋯+nt=n的非负整数解(nn1,n2,⋯ ,nt)x1n1x2n2⋯xtnt(x_1 + x_2 + \cdots + x_t)^n = \sum_{满足 n_1 + \cdots + n_t=n 的非负整数解} \binom{n}{n_1,n_2,\cdots,n_t} x_1^{n_1}x_2^{n_2}\cdots x_t^{n_t}

其中的 (nn1,n2,⋯ ,nt)\dbinom{n}{n_1,n_2,\cdots,n_t} 是多项式系数,它的性质也很相似:

∑(nn1,n2,⋯ ,nt)=tn\sum{\binom{n}{n_1,n_2,\cdots,n_t}} = t^n

排列与组合进阶篇

接下来我们介绍一些排列组合的变种.

多重集的排列数 | 多重组合数

请大家一定要区分 多重组合数 与 多重集的组合数!两者是完全不同的概念!

多重集是指包含重复元素的广义集合.设 S={n1⋅a1,n2⋅a2,⋯ ,nk⋅ak}S=\{n_1\cdot a_1,n_2\cdot a_2,\cdots,n_k\cdot a_k\} 表示由 n1n_1 个 a1a_1,n2n_2 个 a2a_2,…,nkn_k 个 aka_k 组成的多重集,SS 的全排列个数为

n!∏i=1kni!=n!n1!n2!⋯nk!\frac{n!}{\prod_{i=1}^kn_i!}=\frac{n!}{n_1!n_2!\cdots n_k!}

相当于把相同元素的排列数除掉了.具体地,你可以认为你有 kk 种不一样的球,每种球的个数分别是 n1,n2,⋯ ,nkn_1,n_2,\cdots,n_k,且 n=n1+n2+…+nkn=n_1+n_2+\ldots+n_k.这 nn 个球的全排列数就是 多重集的排列数.多重集的排列数常被称作 多重组合数.我们可以用多重组合数的符号表示上式:

(nn1,n2,⋯ ,nk)=n!∏i=1kni!\binom{n}{n_1,n_2,\cdots,n_k}=\frac{n!}{\prod_{i=1}^kn_i!}

可以看出,(nm)\dbinom{n}{m} 等价于 (nm,n−m)\dbinom{n}{m,n-m},只不过后者较为繁琐,因而不采用.

多重集的组合数 1

设 S={n1⋅a1,n2⋅a2,⋯ ,nk⋅ak}S=\{n_1\cdot a_1,n_2\cdot a_2,\cdots,n_k\cdot a_k\} 表示由 n1n_1 个 a1a_1,n2n_2 个 a2a_2,…,nkn_k 个 aka_k 组成的多重集.那么对于整数 r(r<ni,∀i∈[1,k])r(r<n_i,\forall i\in[1,k]),从 SS 中选择 rr 个元素组成一个多重集的方案数就是 多重集的组合数.这个问题等价于 x1+x2+⋯+xk=rx_1+x_2+\cdots+x_k=r 的非负整数解的数目,可以用插板法解决,答案为

(r+k−1k−1)\binom{r+k-1}{k-1}

多重集的组合数 2

考虑这个问题:设 S={n1⋅a1,n2⋅a2,⋯ ,nk⋅ak,}S=\{n_1\cdot a_1,n_2\cdot a_2,\cdots,n_k\cdot a_k,\} 表示由 n1n_1 个 a1a_1,n2n_2 个 a2a_2,…,nkn_k 个 aka_k 组成的多重集.那么对于正整数 rr,从 SS 中选择 rr 个元素组成一个多重集的方案数.

这样就限制了每种元素的取的个数.同样的,我们可以把这个问题转化为带限制的线性方程求解:

∀i∈[1,k], xi≤ni, ∑i=1kxi=r\forall i\in [1,k],\ x_i\le n_i,\ \sum_{i=1}^kx_i=r

于是很自然地想到了容斥原理.容斥的模型如下:

  1. 全集:∑i=1kxi=r\displaystyle \sum_{i=1}^kx_i=r 的非负整数解.
  2. 属性:xi≤nix_i\le n_i.

于是设满足属性 ii 的集合是 SiS_i,Si‾\overline{S_i} 表示不满足属性 ii 的集合,即满足 xi≥ni+1x_i\ge n_i+1 的集合(转化为上面插板法的问题三).那么答案即为

∣⋂i=1kSi∣=∣U∣−∣⋃i=1kSi‾∣\left|\bigcap_{i=1}^kS_i\right|=|U|-\left|\bigcup_{i=1}^k\overline{S_i}\right|

根据容斥原理,有:

∣⋃i=1kSi‾∣=∑i∣Si‾∣−∑i,j∣Si‾∩Sj‾∣+∑i,j,k∣Si‾∩Sj‾∩Sk‾∣−⋯+(−1)k−1∣⋂i=1kSi‾∣=∑i(k+r−ni−2k−1)−∑i,j(k+r−ni−nj−3k−1)+∑i,j,k(k+r−ni−nj−nk−4k−1)−⋯+(−1)k−1(k+r−∑i=1kni−k−1k−1)\begin{aligned} \left|\bigcup_{i=1}^k\overline{S_i}\right| =&\sum_i\left|\overline{S_i}\right| -\sum_{i,j}\left|\overline{S_i}\cap\overline{S_j}\right| +\sum_{i,j,k}\left|\overline{S_i}\cap\overline{S_j}\cap\overline{S_k}\right| -\cdots\\ &+(-1)^{k-1}\left|\bigcap_{i=1}^k\overline{S_i}\right|\\ =&\sum_i\binom{k+r-n_i-2}{k-1} -\sum_{i,j}\binom{k+r-n_i-n_j-3}{k-1}+\sum_{i,j,k}\binom{k+r-n_i-n_j-n_k-4}{k-1} -\cdots\\ &+(-1)^{k-1}\binom{k+r-\sum_{i=1}^kn_i-k-1}{k-1} \end{aligned}

拿全集 ∣U∣=(k+r−1k−1)\displaystyle |U|=\binom{k+r-1}{k-1} 减去上式,得到多重集的组合数

Ans=∑p=0k(−1)p∑A(k+r−1−∑AnAi−pk−1)Ans=\sum_{p=0}^k(-1)^p\sum_{A}\binom{k+r-1-\sum_{A} n_{A_i}-p}{k-1}

其中 A 是充当枚举子集的作用,满足 ∣A∣=p, Ai<Ai+1|A|=p,\ A_i<A_{i+1}.

圆排列

nn 个人全部来围成一圈,所有的排列数记为 Qnn\mathrm Q_n^n.考虑其中已经排好的一圈,从不同位置断开,又变成不同的队列. 所以有

Qnn×n=Ann⟹Qn=Annn=(n−1)!\mathrm Q_n^n \times n = \mathrm A_n^n \Longrightarrow \mathrm Q_n = \frac{\mathrm A_n^n}{n} = (n-1)!

由此可知部分圆排列的公式:

Qnr=Anrr=n!r×(n−r)!\mathrm Q_n^r = \frac{\mathrm A_n^r}{r} = \frac{n!}{r \times (n-r)!}

组合数性质 | 二项式推论

由于组合数在 OI 中十分重要,因此在此介绍一些组合数的性质.

(nm)=(nn−m)(1)\binom{n}{m}=\binom{n}{n-m}\tag{1}

相当于将选出的集合对全集取补集,故数值不变.(对称性)

(nk)=nk(n−1k−1)(2)\binom{n}{k} = \frac{n}{k} \binom{n-1}{k-1}\tag{2}

由定义导出的递推式.

(nm)=(n−1m)+(n−1m−1)(3)\binom{n}{m}=\binom{n-1}{m}+\binom{n-1}{m-1}\tag{3}

组合数的递推式(杨辉三角的公式表达).我们可以利用这个式子,在 O(n2)O(n^2) 的复杂度下推导组合数.

(n0)+(n1)+⋯+(nn)=∑i=0n(ni)=2n(4)\binom{n}{0}+\binom{n}{1}+\cdots+\binom{n}{n}=\sum_{i=0}^n\binom{n}{i}=2^n\tag{4}

这是二项式定理的特殊情况.取 a=b=1a=b=1 就得到上式.

∑i=0n(−1)i(ni)=[n=0](5)\sum_{i=0}^n(-1)^i\binom{n}{i}=[n=0]\tag{5}

二项式定理的另一种特殊情况,可取 a=1,b=−1a=1, b=-1.式子的特殊情况是取 n=0n=0 时答案为 11.

∑i=0k(ni)(mk−i)=(m+nk)(6)\sum_{i=0}^k \binom{n}{i}\binom{m}{k-i} = \binom{m+n}{k}\tag{6}

拆组合数的式子,在处理某些数据结构题时会用到.被称为 范德蒙恒等式.

∑i=0n(ni)2=(2nn)(7)\sum_{i=0}^n\binom{n}{i}^2=\binom{2n}{n}\tag{7}

这是 (6)(6) 的特殊情况,取 n=k=mn=k=m 即可.

∑i=0ni(ni)=n2n−1(8)\sum_{i=0}^ni\binom{n}{i}=n2^{n-1}\tag{8}

带权和的一个式子,通过对 (4)(4) 对应的多项式函数求导可以得证.

∑i=0ni2(ni)=n(n+1)2n−2(9)\sum_{i=0}^ni^2\binom{n}{i}=n(n+1)2^{n-2}\tag{9}

与上式类似,可以通过对多项式函数求导证明.

∑l=0n(lk)=(n+1k+1)(10)\sum_{l=0}^n\binom{l}{k} = \binom{n+1}{k+1}\tag{10}

通过组合分析一一考虑 S={a1,a2,⋯ ,an+1}S=\{a_1, a_2, \cdots, a_{n+1}\} 的 k+1k+1 子集数可以得证,在恒等式证明中比较常用.被称为 朱世杰恒等式.

(nr)(rk)=(nk)(n−kr−k)(11)\binom{n}{r}\binom{r}{k} = \binom{n}{k}\binom{n-k}{r-k}\tag{11}

通过定义可以证明.

∑i=0n(n−ii)=Fn+1(12)\sum_{i=0}^n\binom{n-i}{i}=F_{n+1}\tag{12}

其中 FF 是斐波那契数列.

(n+kk)2=∑j=0k(kj)2(n+2k−j2k)(13)\binom{n+k}{k}^2=\sum_{j=0}^k\binom{k}{j}^2\binom{n+2k-j}{2k}\tag{13}

通过 (6)(6) 可以证明.被称为 李善兰恒等式.

二项式反演

记 fnf_n 表示恰好使用 nn 个不同元素形成特定结构的方案数,gng_n 表示从 nn 个不同元素中选出 i≥0i \geq 0 个元素形成特定结构的总方案数.

若已知 fnf_n 求 gng_n,那么显然有:

gn=∑i=0n(ni)fig_n = \sum_{i = 0}^{n} \binom{n}{i} f_i

若已知 gng_n 求 fnf_n,那么:

fn=∑i=0n(ni)(−1)n−igif_n = \sum_{i = 0}^{n} \binom{n}{i} (-1)^{n-i} g_i

上述已知 gng_n 求 fnf_n 的过程,就称为 二项式反演.

证明

将反演公式的 gig_i 展开得到:

fn=∑i=0n(ni)(−1)n−i[∑j=0i(ij)fj]=∑i=0n∑j=0i(ni)(ij)(−1)n−ifj\begin{aligned} f_n &= \sum_{i = 0}^{n} \binom{n}{i} (-1)^{n-i} \left[\sum_{j = 0}^{i} \binom{i}{j} f_j\right] \\ &= \sum_{i = 0}^{n}\sum_{j = 0}^{i}\binom{n}{i}\binom{i}{j} (-1)^{n-i}f_j \end{aligned}

先枚举 jj,再枚举 ii,得到:

fn=∑j=0n∑i=jn(ni)(ij)(−1)n−ifj=∑j=0nfj∑i=jn(ni)(ij)(−1)n−i\begin{aligned} f_n &= \sum_{j = 0}^{n}\sum_{i = j}^{n}\binom{n}{i}\binom{i}{j} (-1)^{n-i}f_j \\ &= \sum_{j = 0}^{n}f_j\sum_{i = j}^{n}\binom{n}{i}\binom{i}{j} (-1)^{n-i} \end{aligned}

使用 「组合数性质 | 二项式推论」 的公式 (11) 得到:

fn=∑j=0nfj∑i=jn(nj)(n−ji−j)(−1)n−i=∑j=0n(nj)fj∑i=jn(n−ji−j)(−1)n−i\begin{aligned} f_n &= \sum_{j = 0}^{n}f_j\sum_{i = j}^{n}\binom{n}{j}\binom{n - j}{i - j} (-1)^{n-i} \\ &= \sum_{j = 0}^{n}\binom{n}{j}f_j\sum_{i = j}^{n}\binom{n - j}{i - j} (-1)^{n-i} \end{aligned}

令 k=i−jk = i - j.则 i=k+ji = k + j,上式转换为:

fn=∑j=0n(nj)fj∑k=0n−j(n−jk)(−1)n−j−k1kf_n = \sum_{j = 0}^{n}\binom{n}{j}f_j\sum_{k = 0}^{n - j}\binom{n - j}{k} (-1)^{n-j-k}1^{k}

使用 「组合数性质 | 二项式推论」 的公式 (5) 得到:

fn=∑j=0n(nj)fj[n=j]=fnf_n = \sum_{j = 0}^{n}\binom{n}{j}f_j[n = j] = f_n

证毕.

练习

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

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