【组合数学：三】容斥原理

容斥原理及其应用

容斥原理

• 之前我们已经见过，对集合中对象个数的间接计数要比对这些对象直接计数容易，即利用减法原理
• 例子：$1\sim n$ 的全排列，且 $1$ 不在第一个位置的方案数
  解一：第一个位置选取除了 $1$ 的数，方案数 $n-1$ ，剩下的全排列 $(n-1)!$
  解二：全集：$n!$ ，第一个位置放 $1$ 的排列方案为 $(n-1)!$，所以方案为 $n!-(n-1)!$
• 例子：$1\sim 600$ 之间不能被 $6$ 整除的整数的个数
  解：$600-600/6=500$
• 上面是关于容斥定理的最简单的例子
  现在我们把它扩展到容斥定理
  设 $P_i$ 表示 $S$ 对象涉及到的性质
  并设 A i = { x : x 属 于 S 且 x 具 有 性 质 P i } A_i=\{x:x属于S且x具有性质P_i\} Ai​={x:x属于S且x具有性质Pi​}
  比如 $A_i\cap A_j$ 就是同时具有性质 $P_i,P_j$
  不具有任何性质的对象形成的子集则是 ${\cap }_i\overline{A_i}$
• 定理：集合 $S$ 中不具有性质 $P_1\sim P_m$ 的对象个数由下面的交错表达式给出：
  $$|{\cap }_i\overline{A_i}|=|S|-\sum |A_i|+\sum |A_i\cap A_j|-\sum |A_i\cap A_j\cap A_k|+\cdots +(-1{)}^m|{\cap }_i{A}_i|$$
  右边的项数为 $2^m$
  证明：若 $x$ 不具备任何一个性质 $P_i$ ，那么 $x$ 对式子的贡献是 $1-0+0\cdots =1$
  若 $x$ 具有 $y$ 个性质，$x$ 对于 $|S|$ 的贡献是 $1$
  它对 $\sum |A_i|$ 的贡献是 $C_n^1$
  它对 $\sum |A_i\cap A_j|$ 的贡献是 $C_n^2$
  于是，它的总的贡献就是 $\sum (-1{)}^i{C}_n^i=0$，于是得证
• 推论：集合 $S$ 中至少具有性质 $P_1\sim P_m$ 之一的对象个数由下式给出：
  $$|{\cup }_i{A}_i|=\sum |A_i|-\sum |A_i\cap A_j|+\cdots +(-1{)}^{m+1}|{\cap }_i{A}_i|$$
  证明：根据德摩根定理，即：
  $${\overline{\cup A}}_i=\cap \overline{A_i}$$
  所以等式左边 $|\cup A_i|=|S|-|\overline{\cup A_i}|=|S|-|\cap \overline{A_i}|$
  然后带入上面容斥原理式子即可得到
• 例子：$1\sim 1000$ 之间不能被 $5、6、8$ 整除的数的个数？
  解：设出 $A_i$ 表示 $S$ 中具有 $P_i$ 性质的集合
  容易得到我们最终的答案为：
  $$|\cap \overline{A_i}|$$
  容易得到：
  $$\begin{gathered} |A_1|=\lfloor \frac{1000}{5}\rfloor =200 \\ |A_2|=\lfloor \frac{1000}{6}\rfloor =166 \\ |A_3|=\lfloor \frac{1000}{8}\rfloor =125 \end{gathered}$$
  集合 $|A_i\cap A_j|$ 就是满足 $P_i,P_j$ 的性质的数字，也就是数字被 $lcm(数i,数j)$ 整除的数的个数
  $$\begin{gathered} |A_1\cap A_2|=\lfloor \frac{1000}{lcm(5,6)}\rfloor =33 \\ |A_1\cap A_3|=\lfloor \frac{1000}{lcm(5,8)}\rfloor =25 \\ |A_2\cap A_3|=\lfloor \frac{1000}{lcm(6,8)}\rfloor =41 \end{gathered}$$
  最后还有
  $$|A_1\cap A_2\cap A_3|=\lfloor \frac{1000}{lcm(5,6,8)}\rfloor =8$$
  根据容斥定理，得到：
  $$|\cap \overline{A_i}|=|S|-\sum |A_i|+\sum |A_i\cap A_j|-\sum |A_i\cap A_j\cap A_k|=600$$
• 后文记号：
  假设 ${\cap }_{i=1}^k{A}_i$ 的大小仅依赖于 $k$ 而不依赖于我们选了哪 $k$ 个集合，因此就存在常数 ${\alpha }_i$ 满足如下：
  $$\begin{array}{rl}{\alpha }_0& =|S|\\ {\alpha }_1& =|A_1|=|A_2|=\cdots =|A_m|\\ {\alpha }_2& =|A_1\cap A_2|=|A_2\cap A_3|=\cdots =|A_{m-1}\cap A_m|\\ ⋮& \\ {\alpha }_m& =|A_1\cap A_2\cap \cdots \cap A_m|\end{array}$$
  在这个情况下，容斥原理可以简化成：
  $$|\cap \overline{A_i}|=\sum _{i=0}^m(-1{)}^i{\alpha }_i{C}_m^i$$
• 例子：$0\sim 99999$ 之间有多少个数字含有 $2,5,8$ 至少一个的整数？
  解：我们可以把 $S$ 看做是 $0\sim 9$ 的一个多重集合的 $5$ 排列，每个数字的重数至少是 $5$
  设 $P_1,P_2,P_3$ 表示一个整数不包含 $2/5/8$ 的性质
  我们希望得到的是 $|\overline{A_1}\cup \overline{A_2}\cup \overline{A_3}|$
  根据上文的记号，我们有：
  $$\begin{gathered} {\alpha }_0=10^5 \\ {\alpha }_1=9^5 \\ {\alpha }_2=8^5 \\ {\alpha }_3=7^5 \end{gathered}$$
  于是我们带入容斥定理，就能得到答案：
  $$C_3^{0}10^5-C_3^1{9}^5+C_3^2{8}^5-C_3^3{7}^5$$

带重复的组合

• 之前，我们已经得到了，有 $k$ 种对象，每种对象都有无限重数的多重集合的 $r$ 组合的个数等于
  $$C_{k+r-1}^r$$
  在这里，我们利用上述公式与容斥原理的应用，给出寻找元素重数没有限制的多重集合的 $r$ 组合的方法
• 例子：确定多重集合 T = { 3 ⋅ a , 4 ⋅ b , 5 ⋅ c } T=\{3\cdot a,4\cdot b,5\cdot c\} T={3⋅a,4⋅b,5⋅c} 的 $10$ 组合的数目
  我们把容斥原理应用到多重集合 T ∗ = { ∞ ⋅ a , ∞ ⋅ b , ∞ ⋅ c } T^*=\{\infin \cdot a,\infin \cdot b,\infin \cdot c\} T∗={∞⋅a,∞⋅b,∞⋅c} 上
  设 $P_1$ 是 $T^{\ast }$ 的 $10$ 组合中 $a$ 出现多于 $3$ 次的性质
  设 $P_2$ 是 $T^{\ast }$ 的 $10$ 组合中 $b$ 出现多于 $4$ 次的性质
  设 $P_3$ 是 $T^{\ast }$ 的 $10$ 组合中 $c$ 出现多于 $5$ 次的性质
  容易得到，$T$ 的 $10$ 组合就是 $T^{\ast }$ 的 $10$ 组合中，不满足上述三个性质的数目。
  设 $A_i$ 由 $T^{\ast }$ 的 $10$ 组合中不具有性质 $P_i$ 组成的
  首先容易得到：
  $$|S|=C_{10+3-1}^{10}=66$$
  那么集合 $A_1$ 就是 $a$ 至少出现了 $4$ 次的 $10$ 组合，我们去掉这 $4$ 个 $a$ ，就变成了一个 $T^{\ast }$ 的 $6$ 组合的方案数：
  $$|A_1|=C_{3+6-1}^6=28$$
  类似的也能得到 $|A_2|,|A_3|$
  那么集合 $|A_1\cap A_2|$ 就是 $a$ 至少出现了 $4$ 次， $b$ 至少出现了 $5$ 次 的 $T^{\ast }$ 的 $10$ 组合
  我们可以去掉 $4$ 个 $a$ 和 $5$ 个 $b$ ，可以得到是一个 $T^{\ast }$ 的 $1$ 组合，就是：
  $$|A_1\cap A_2|=C_{1+3-1}^1=3$$
  其他的都是类似的做法
  然后使用容斥原理就可以得到我们要的答案：
  $$|\overline{A_1}\cap \overline{A_2}\cap \overline{A_3}|=|S|-\sum |A_i|+\sum |A_i\cap A_j|-\sum |A_i\cap A_j\cap A_k|$$
• 例子：满足
  $$\{\begin{array}{l}{x}_1+x_2+x_3+x_4=18\\ 1\le x_1\le 5\\ -2\le x_2\le 4\\ 0\le x_3\le 5\\ 3\le x_4\le 9\\ x_1,x_2,x_3,x_4\in Z\end{array}$$
  的解的个数？
  解：引入新变量 $y_1=x_1-1$ 等，把问题转化为：
  $$\{\begin{array}{l}{y}_1+y_2+y_3+y_4=16\\ 0\le y_1\le 4\\ 0\le y_2\le 6\\ 0\le y_3\le 5\\ 0\le y_4\le 6\\ y_1,y_2,y_3,y_4\in Z\end{array}$$
  的方案数
  然后方法类似于上面，直接容斥定理一波带走

错位排列

• 因为对象的实际性质与讨论不相干，我们可以取 $X$ 为集合 { 1 , 2 , 3 , ⋯   , n } \{1,2,3,\cdots,n\} {1,2,3,⋯,n} ，其中每个整数的位置都由它们在序列 $1,2,\cdots ,n$ 中的位置确定。
  { 1 , 2 , 3 , ⋯   , n } \{1,2,3,\cdots,n\} {1,2,3,⋯,n} 的一个错位排列是 { 1 , 2 , 3 , ⋯   , n } \{1,2,3,\cdots,n\} {1,2,3,⋯,n} 的一个排列 $i_1,i_2,\cdots ,i_n$ ，其中满足
  $$\forall k\in [1,n],i_k\ne i$$
  记 $D_n$ 表示 { 1 , 2 , 3 , ⋯   , n } \{1,2,3,\cdots,n\} {1,2,3,⋯,n} 的错位排列的数目
• 定理：对于 $n\ge 1$，我们有：
  $$D_n=n!(1-\frac{1}{1!}+\frac{1}{2!}-\frac{1}{3!}+\cdots +(-1{)}^n\frac{1}{n!})$$
  证明：设 $P_j$ 性质表示 $i_j=j$ 的时候的一个排列
  错位排列就是 $|\cap \overline{A_i}|$
  容易想到，
  $$\begin{gathered} |A_i|=(n-1)! \\ |A_i\cap A_j|=(n-2)! \\ \cdots \\ |A_{i_1}\cap A_{i_2}\cap \cdots \cap A_{i_k}|=(n-k)! \end{gathered}$$
  直接应用容斥定理，我们得到：
  $$\begin{array}{rl}|\cap \overline{A_i}|& =|S|-\sum |A_i|+\sum |A_i\cap A_j|+\cdots \\ & =n!-C_n^1(n-1)!+C_n^2(n-2)!+\cdots \\ & =n!-\frac{n!}{1!}+\frac{n!}{2!}-\frac{n!}{3!}+\cdots +(-1{)}^n\frac{n!}{n!}\\ & =n!(1-\frac{1}{1!}+\frac{1}{2!}-\frac{1}{1!}+\cdots +(-1{)}^n\frac{1}{n!})\end{array}$$
• 我们注意到：$e^{-1}\approx \frac{D_n}{n!}$
  关系式：$D_n=(n-1)(D_{n-2}+D_{n-1})$
  或者等价的 $D_n=n{D}_{n-1}+(-1{)}^n$
  在下一章讨论求解常系数线性递推关系

带有禁止位置的排列

• 例子：在 $n\times n$ 的网格中，令 $X_i$ 表示第 $i$ 行不能放的所有列的集合
  现在问你，放 $n$ 个非攻击型车在能放的位置的方案数
  解：令 $P_i$ 表示第 $i$ 行的车是属于 $X_i$ 的列上的性质
  我们最终求的就是：
  $$|\cap \overline{A_i}|$$
  容斥原理带入，我们需要求出：
  $$\begin{gathered} |S|=n! \\ |A_1|=|X_1|(n-1) \\ |A_i|=|X_i|(n-1)! \end{gathered}$$
  令 $r_1=\sum |X_i|$，我们可以得到：
  $$\sum |A_i|=r_1(n-1)!$$
  设 $r_2$ 表示把两个非攻击型车放在禁止位置上的方案数，则：
  $$\sum |A_i\cap A_j|=r_2(n-2)!$$
  顺理推到 $r_k$ 表示把 $k$ 个非攻击型车放在禁止位置上的方案数，则：
  $$|\cap \overline{A_i}|=\sum _{i=0}^n(-1{)}^i{r}_i(n-i)!$$
  注意到，如果禁止格子特别多的话，那么 $r_k$ 反而更难计算。

莫比乌斯反演

• 容斥原理是莫比乌斯反演在有限偏序集上的一个实例
  跳过大量证明…
  具体的 $ACM$ 竞赛内容可以看我之前的博客
• 例子：计算 $k$ 个不同元素的循环 $n$ 排列的个数，每个元素都可以用任意多次

---

更新了自己的理解

---

• 解：我们设 $d$ 表示一个循环排列的周期（每个元素都绕中心旋转最少的次数让出现两次相同排列）
  一个循环 $n$ 排列的周期 $d$ 满足 $1\le d\le n$ 且 $d|n$
  设 $h(n)$ 是可能使用这 $k$ 个元素的循环 $n$ 排列的个数
  设 $f(d)$ 等于长度为 $n$ 且可能使用这 $k$ 个元素的串的个数，每个串有一个循环周期 $d$，满足 $d|n$
  因此下式成立：
  $$h(n)=\sum _{d|n}\frac{f(d)}{d}$$
  因为循环周期为 $d$，所以有 $d$ 次重复计数。
  设 $g(n)$，表示：
  $$g(n)=\sum _{e|n}f(e)$$
  于是，$g(n)$ 是长度为 $n$ 的串的总数，因为一个串的循环周期一定满足 $d|n$ ，且 $g(n)=n^k$
  根据经典的莫比乌斯反演，我们得到：
  $$f(d)=\sum _{e|d}\mu (\frac{d}{e})g(e)$$
  在第一个式子中我们使用上述式子，得到：
  $$\begin{array}{rl}h(n)& =\sum _{d|n}\frac{f(d)}{d}\\ & =\sum _{d|n}\frac{1}{d}\sum _{e|d}\mu (\frac{d}{e})k^e\\ & =\sum _{e|n}(\sum _{m|\frac{n}{e}}\frac{1}{me}\mu (m))k^e\\ & =\sum _{e|n}(\sum _{r|\frac{n}{e}}\frac{r}{n}\mu (\frac{n}{er}))k^e\\ & =\sum _{e|n}\frac{\phi (n/e)}{n}{k}^e\\ & =\frac{1}{n}\sum _{e|n}\phi (\frac{n}{e})k^e\end{array}$$
  注意到，第二行到第三行的转换：因为 $d|n$ 且 $e|d$ ，所以相当于 $d=me$ ，从而 $me|n$，从而 $m|(n/e)$
  第三行到第四行的转换：令 $n=rd=mer$，所以得到了 $r=\frac{n}{me}$，那自然 $r|\frac{n}{e}$
  第四行到第五行的转换：因为我们有
  $$\phi (n)=\sum _{d|n}\mu (\frac{n}{d})d$$
  注意到，这个问题就是珠子染色问题，但是我们没有使用 $Polya$ 定理或者 $BurnSide$ 引理，因此化简会非常复杂
  更简单的化简方式这篇博客有：【算法讲27：Polya定理】BurnSide 引理 Polya定理 | 洛谷 P4980