组合数学学习笔记(二):鸽巢原理、容斥及反演

鸽巢原理

i=1npin+1 个物品放入 n 个盒子,一定存在一个盒子 i,使得第 i 个盒子至少装了 pi 个物品。

证明:考虑反证法。

假设每个盒子 i 都装了小于 pi 个物品,那么物品数量最大当且仅当每个盒子 i 都装了 pi1 个物品,那么物品数就是 (p11)+(p21)+(p31)++(pn1)=i=1npin,此时物品数小于物品总数,矛盾,因此 一定存在一个盒子 i,使得第 i 个盒子至少装了 pi 个物品。

推论: 如果把 n(r1)+1 个物品放入 n 个盒子中,那么一定有一个盒子中至少装了 r 个物品。

证明:把所有的 pi 取成一样就可以了。

练习一

有十个数 a1,a2,a10 满足 1i101ai60,证明能够从 ai 中挑出两个交为空的子集,使得它们的和相等。

证明:首先,由于 1ai60,那么这个子集的和就有可能是 0600601 个数中的一个,而可能的子集数为 210=1024,根据鸽巢原理,一定有两个子集的和相等。

如果这两个集合没有交集,那么已经满足题目的要求;如果这两个集合有交集,那就将这两个集合同时减掉交集,此时两个集合的和依然相等,而且仍然为全集的子集,原命题证明成立。

练习二

证明一张有超过 1 个点的简单无向图必定有两点度数相等。

解:考虑反证法。

一个 n 个点的简单无向图,每个点的度数是 0n1 中的一个,假设这 n 个点的度数恰好是 0n1 各一个,但由于度数为 n1 的点一定向其它 n1 个点个连了一条边,因此不存在度数为 0 的点,矛盾。于是度数 0n1 不能同时出现,因此度数的选择只有 n1 种。

根据鸽巢原理,这 n 个点中一定有两个点度数相等,原命题证明成立。

练习三

前置知识:线性代数学习笔记(一)(2024.7.24)

证明能从任意 11 个实数中挑选出 4 个数 a,b,c,d 满足:

(ac+bd)212(a2+b2)(c2+d2)

证明:考虑两个向量 u=[ab],v=[cd],它们的点积为 [ab][cd]=ac+bd,那么右边又可以改写成 (uv)2 这是左边的几何含义。

这两个向量的模长分别是 a2+b2c2+b2,那么右边又可以改写成 12|u|2|v|2,这是右边的几何含义。

考虑到点积的几何意义是将一个向量 u 投影到另一个向量 v 上,再用投影长度乘以 v 的模长,设这两个向量的夹角为 θ, 那么点积可以表示成 |u||v|cosθ,那么左边又可以写成 |u|2|v|2cos2θ

这时原先的不等式变成了 |u|2|v|2cos2θ12|u|2|v|2,消去相同的项,可以得到 cos2θ12,即 cosθ22cosθ22,即 θ[0,π4][3π4,π]

现在考虑将平面每旋转 π4 划分一次,将平面划分成 8 个部分,由于有 11 个数字,根据鸽巢原理,一定有 2 个向量位于同一个部分内,那 θ 一定满足要求,原命题证明成立。

练习四

证明由任意 n2+1 个互不相同的实数构成的序列 a1,a2,,an2+1,要么含有长度为 n+1 的递增子序列,要么含有长度为 n+1 的递减子序列。

证明:先假设原序列不存在长度为 n+1 的递增子序列,来证明原序列存在长度为 n+1 的递减子序列。

mk 表示从 ak 开始的最长递增子序列的长度,那么根据假设,k[1,n2+1],mkn,又因为 mk1,所以 m1,m2,,mn2+11n 之间的整数,根据鸽巢原理,这 n2+1 个数中必然有 n+1 个是相同的数。

假设相同的 m 分别是 mk1=mk2==mkn+1,其中 1k1<k2<<kn+1n2+1,再假设 i[1,n],aki<aki+1,由于 ki<ki+1,那么可以在 aki+1 的最长递增子序列前面加上 aki 使得 mki>mki+1,这与假设矛盾,因此不存在 aki<aki+1,于是 ak1>ak2>>akn+1,此时这是一个长度为 n+1 的递减子序列,原命题证明成立。

容斥原理

其实这个东西我们在小学时就学过了,比如题目,其实这就是容斥原理。一般,如果我们要不重不漏地统计至少具有一种特性的某些物品的系数,就需要考虑如何减掉重复的,如何加上漏掉的,这时就需要用容斥原理来统计。用计算式来表达就是:

|i=1nAi|=j=1n(1)j11k1<k2<k3<<kjn|l=1jAkl|

用人话说就是如果你想求出一些集合的并集的大小,你只需要加上这些集合的大小,减掉任意两个集合的交集,加上任意 3 个集合的交集,重复这种奇数个集合就加,偶数个集合就减的方式,就可以得到,下面来证明为什么。

定理2.1

P1,P2,,Pm 为研究对象的 m 种性质,Am 为具有 Pm 性质的对象所构成的集合,AiAjAk 表示同时具有 Pi,Pj,Pk 性质的对象所具有的,A1AjAm,表示不具有任何性质的对象所构成的集合。

那么不具有性质 P1,P2,,Pm 的对象个数由下面的交错表达式给出 |AiAjAm|=|S||Ai|+|AiAj||AiAjAk|++(1)m|AiAjAm|

证明:先考虑只有两个性质的情况,那么 |A1A2|=|S||A1||A2|+|A1A2|,现在我们只需要证明,一个既不具有 P1 性质,也不具有 P2 性质的元素对于右边式子的贡献为 1,而具有 P1P2 性质的元素对于右边式子的贡献为 0

假设元素 x 既不具有 P1 性质,也不具有 P2 性质,那么它只被包含在 S 中,对于右边式子的贡献为 100+0=1

假设元素 x 只具有 P1 性质,那么它被包含在 S,A1 中,对于右边式子的贡献为 110+0=0

假设元素 x 只具有 P2 性质,那么它被包含在 S,A2 中,对于右边式子的贡献为 101+0=0

假设元素 x 既具有 P1 性质,也具有 P2 性质,那么它被包含在 S,A1,A2,A1A2 中,对于右边式子的贡献为 111+1=0

那么右边等式计算出的结果就是既不具有性质 P1,也不具有性质 P2 的元素的个数。

现在考虑拓展到一般,假设元素 x 不具有任意一个性质,那么他对右边的贡献就是 10+00++(1)m0=1

考虑一个具有 n(nm) 个性质的元素 y,我们可以从这 n 个性质中任意选出 k 个来,那么 A1A2Ak 这个集合中,而一共有 (nk) 个包含 y 且大小为 k 的集合你,那么答案就是 k=1m(1)k(nk)=k=0n(1)k(nk),套用二项式定理,这个式子 =[1+(1)]n=0,于是 y 对右边式子的贡献为 0

现在,只有不具有任何性质的元素对等式右边有贡献,定理 2.1 证明成立。

感性理解就是减多了就加,加多了就减。

定理2.2(容斥原理)

|A1A2Am|=|Ai||AiAj|+|AiAjAk|(1)m+1|AiAjAm|

证明:将定理1.1两边同时用 |S| 减即可证明。

容斥系数

例题三

如何现在有一个序列 an0anm,请问有多少个序列满足:不存在两个连续的0?

例题四

问对于所有长为 n 的排列,有多少排列存在一个连续上升段 k。对所有 k 回答,对大质数取模。

例题五

反射容斥

例题六

例题七

反演

广义的反演,就是两个函数之间互相推倒的关系,比如 f(x)=xg(x),那么 g(x)=f(x)x,放到离散意义下,就可以表示两个数列之间之间互相推倒的关系,比如 F(n)=i=1nG(i)(前缀和),那么 G(n)=F(n)F(n1)(差分),这就是信息学中常见的反演。

考虑如果将第二个式子带入第一个式子,那么会得到 F(n)=F(n) 这样一个无聊的等式,相当于是什么也没有做,这和线性代数中的线性变换与逆变换极其类似,这也可以从矩阵的角度分析,假设 F(n)=AG(n),那么一定有 G(n)=A1F(n),若有多个反演叠加,将它们的转移矩阵全部乘在一起就行了。

信息学竞赛中常见的反演有以下几种:

  • 二项式反演:与容斥原理有关,一般解决“恰好”与“至多(少)”、“给定”与“任意”之间的转换关系;

  • Min-Max 反演:与容斥原理有关,将 minmax 表示,或将 maxmin 表示,在某些题中可以改变枚举顺序或枚举类型;

  • 单位根反演:与多项式有关,DFT 对单位根多点求值后,用单位根反演重新插值出原函数;

  • 莫比乌斯反演:与数论有关,一般用于含有 gcd 的式子的简化;

  • 斯特林反演:与组合数学有关,由于第二类斯特林数是好求的,运用斯特林反演之后可以计算出第一类斯特林数,在上升幂、普通幂、下降幂之间的转化有用处。

  • 子集反演:与高维前缀和有关系,一般结合快速莫比乌斯变换或快速沃尔什变换可以处理高位空间中特殊的数点问题;

二项式反演

先来看一下最基本的式子:

F(n)=i=mn(ni)G(i)G(n)=i=mn(1)ni(ni)G(i)

这个式子有一个等价形式,就是

F(n)=i=mn(1)i(ni)G(i)G(n)=i=mn(1)i(ni)G(i)

这个式子可以说是非常的对称,根据反演可以看作矩阵乘法,那么该矩阵 A×A=I

我们可以先从线性代数的角度来证明一下。

证:

或者可以最暴力地带入证明

证:F(n)=i=mn(ni)j=mi(1)ij(ij)F(j)

考虑改变求和顺序,那么原式

=j=mnF(j)i=jn(ni)(ij)(1)ij

套用三项式系数恒等式,那么原式

=j=mnF(j)i=jn(nj)(njij)(1)ij

再次调整计算顺序,那么原式

=j=mn(nj)F(j)i=jn(njij)(1)ij

将枚举 i 替换为枚举 ij,那么原式

=j=mn(nj)F(j)i=0nj(nji)(1)i

发现内层求和是一个二项式定理的形式,可以写成 [1+(1)]nj 形式,而这个式子只有在 nj=0 时才有值 1,于是原式最终写成

j=mn(nj)F(j)[n=j]=F(n)

F(n)=i=mn(im)G(i)G(n)=i=mn(1)im(im)F(i)

第二种形式证明类似。

练习(洛谷什么时候把这道题传上去的)

n 个元素,问有多少种选择若干个子集的方案,使得选出的子集的交集大小恰好为 k

0<kn106

解:直接钦定 i 个的方案数是 F(i)=(ni)(22ni1),设答案为 G(i),观察到

F(i)=j=in(ji)G(j)

那么反演即可得到:

G(k)=i=kn(ik)F(i)

=i=kn(1)ik(ik)(ni)(22ni1)

这样之后直接做就行了。

例题一

一句话题意:有两个序列 ai,bi 保证所有元素互不相同。你需要重排 b 序列,使得恰好有 ki 满足 ai>bi

0<kn2000

解:先将 a 序列排序。

考虑设 dpi,j 表示考虑了前 i 对数,恰有 jax>bx,发现完全没法转移(

主要是你配大的和配小的总有一个转不了,所以你考虑只算其中一个,也就是钦定有 j 对数满足 ax>bx,发现这样就可以转了。

dpi,j=dpi1,j+dpi1,j1×(cnt(ai)j+1)

cnt(ai) 表示 b 序列中比 ai 小的数的个数。

然后设 f(i) 表示钦定了 i 对的方案数,g(i) 表示恰好 i 对的方案数。观察可得:

(ni)!dpn,i=f(i)=j=in(ji)g(j)

g(k)=i=kn(1)ik(ik)(ni)!dpn,i

就可以直接做了。

例题二

一句话题意:有一个 n×n 的矩阵,将其三染色,使得至少有一行或者一列同色,问方案数。

n106

解:直接钦定 ij 列同色:

F(i,0)=3(ni)n+i(ni)

F(0,j)=3(nj)n+j(nj)

F(i,j)=3(ni)(nj)+1(ni)(nj)

考虑恰好 ij 列同色:

G(0,0)=3n2+i=1n(1)ni(ni)(F(0,i)+F(i,0))+i=1nj=1n(1)2nij(i0)(j0)F(i,j)

后面那坨东西带进去:

i=1nj=1n3n2ninj+ij+1(1)2nij(ni)(nj)

=3n2+1i=1n3ni(1)ni(ni)j=1n(3ni)j(1)nj(nj)

=3n2+1i=1n3ni(1)ni((3ni1)n(1)n)

然后就做完了。

min / max反演(容斥)

maxS=TS(1)|T|1minTmaxkthS=TS(1)|T|k(|T|1k1)minT

min/max 反过来也成立。

证:maxkthS=TS(1)|T|k(|T|1k1)minT

=xSxxTS(1)|T|k(|T|1k1)[minT=x]

f(x) 表示 S 中大于 x 的元素构成的集合。

=xSxxTf(x)(1)|T|k(|T|1k1)

枚举 |T|:

=xSxl=1|f(x)|(1)lk(|f(x)|1l1)(l1k1)

=xSxl=1|f(x)|(1)lk(|f(x)|1k1)(|f(x)|klk)

=xSx(|f(x)|1k1)l=1|f(x)|(1)lk(|f(x)|klk)

易知 |f(x)|<k 时无贡献。

=xSx(|f(x)|1k1)l=0|f(x)|k(1)l(|f(x)|kl)

=xSx(|f(x)|1k1)[|f(x)|=k]

=maxkthS

练习

给定三个序列 ai,bi,ci,求

1i<jnmaxai+aj,bi+bj,ci+cjminai+aj,bi+bj,ci+cj

n2×105

解:暴力拆开那个 max :

maxai+aj,bi+bj,ci+cj

=minai+aj+minbi+bj+minci+cjminai+aj,bi+bjminai+aj,ci+cjminbi+bj,ci+cj+minai+aj,bi+bj,ci+cj

抵消掉最后那一项,剩下的项中,只有一个的是平凡的,有两个的可以二维偏序,总复杂度就是 O(nlogn)

例题三

一句话题意:给 n 个元素,每次会随机选择一个,有 pi 的概率选择第 i 个,问第一次所有元素都被选择过的期望时间。

1n20

解:因为 min/max 容斥都是线性运算,且期望具有线性性,所以可以直接套上期望:

E(maxS)=TS(1)|T|1E(minT)

E(minT) 的含义实际上就是第一次选到 T 中元素的期望时间,可以把 TT 的补集看做两个整体,那么一次选中的概率就是 1xTpx,期望时间就是其倒数 xTpx

然后直接带进式子里计算即可,O(n2n)

例题四

一句话题意:给 n 个元素,每次会随机选择一个,有 Mpi 的概率选择第 i 个,问第一次有 k 元素被选择过的期望时间。

1ln103,n10kn,i=1npi=M104

解:相当于求期望出现时间第 nk 大,令 knk,一样的期望式子列出来:

E(maxkthS)=TS(1)|T|1(|T|1k1)mxTpx

直接求肯定 G。因为 M104,所以考虑计算每一种 xTpx 作为分母的项的系数之和,再算答案。

考虑设 dpi,j 表示前 i 个,分母和为 j 的项的系数和。

不选很好转移,选的话中间那坨组合数就不太能转的动。

考虑到 (|T|1k1)=|T|1|T|k(|T|2k1),所以我们可以在方程里及一个 l=|T|,那就可以转了:

dpi,j,l=dpi1,j,ll1lkdpi1,jpi,l1

但是这个方程是 O(n2m) 的,显然过不了,且没用到 k10 的性质。

又考虑到 (|T|1k1)=(|T|2k1)+(|T|2k2)

那么只需要在方程中记 l=k 即可转移:

dpi,j,l=dpi1,j,l+dpi1,jpi,l1dpi1,jpi,l

就可以 O(nmk) 了。

子集反演

快速沃尔什变换(FWT)

快速莫比乌斯变换(FMT)

参考资料

  • Introductory Combinatorics(Fifth Edition) Richard A.Brualdi

  • Discrete Mathematics and Its Applications(Eighth Edition) Kenneth H.Rosen

  • 炫酷反演魔术 command_block

posted @   JPGOJCZX  阅读(17)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 开源Multi-agent AI智能体框架aevatar.ai,欢迎大家贡献代码
· Manus重磅发布:全球首款通用AI代理技术深度解析与实战指南
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· AI技术革命,工作效率10个最佳AI工具
点击右上角即可分享
微信分享提示