组合数学总结

本篇文章是为了总结一下最近做的组合数学的题目以及涉及到的知识点,以后可能会不定期补充。同时也参考了大佬 lyt 的博客。同时如果想要更深入地了解组合数学中的推式子技巧,可以看下 pjw 的博客。地址1 地址2

基本知识

排列数 Anm 表示从 n 个不同的物品中选出 m 个组成的不同排列的数量。Anm=n!(nm)!

组合数(高中常用 Cnm 来表示,但是目前通常用 (nm) 表示),表示从 n 个不同物品中选出 m 个物品的方案数。(nm)=n!m!(nm)!

二项式定理:(x+y)n=i=0n(ni)xiyni

一些比较重要的组合数等式

  1. (nm)=(nnm)
  2. i=0n(ni)=2n。二项式定理当 x=y=1 时的情况。
  3. i=0n(1)i(ni)=[n=0]。二项式定理当 x=1,y=1 时的情况。
  4. (nm)(mk)=(nk)(nkmk)。从组合意义上即可证明。
  5. i=0k(ni)(mki)=(n+mk)。范德蒙德卷积,组合意义上证明即可。
  6. i=0n(im)=(n+1m+1)
  7. (nm)=(n1m1)+(n1m)

二项式反演

fn 为恰好 n 个元素组成特定结构的方案数,gn 为从 n 个元素中选出 i 个元素组成特定结构的方案数,显然 gn=i=0n(ni)fi。则 fn=i=0n(1)ni(ni)gi

fi 恰好,gi 至多)


证明如下(直接从十二重计数法那搬过来的)

i=0n(ni)(1)nif(i)=i=0n(ni)(1)nij=0i(ij)g(j)=i=0nj=0i(ni)(ij)(1)nig(j)=i=0nj=0i(nj)(njij)(1)nig(j)=j=0n(nj)g(j)i=jn(njij)(1)ni=j=0n(nj)g(j)i=0nj(nji)(1)nij=j=0n(nj)g(j)×0nj=g(n)

据说也可以用 EGF 证,但是我不会(


二项式反演的另一种常见形式:

fi 为恰好 i 个条件满足的方案数,gi 为至少 i 个条件满足的方案数(也就是强制让 i 个条件满足后其他的随便),则显然 gk=i=kn(ik)fi,那么 fk=i=kn(1)ik(ik)gi

fi 恰好,gi 至少)

同时,当 k=0 时,f0=i=0n(1)igi。因此容斥原理其实是二项式反演的一种特殊情况。

多重集排列数

设多重集 S={n1x1,n2x2,,nkxk},则该多重集的全排列数量为 n!n1!n2!n3!nk!(n=ni)

多重集组合数

从集合 S={n1x1,n2x2,,nkxk} 中选出 m 个数的方案数。

我们可以把问题转化成另一个形式:n1+n2++nk=m,求正整数解的个数。

多重集组合数分两种情况,没有上界和有上界。

先分析没有上界的情况,也就是 ni= 时,我们可以把它转化成插板法,方案数即为 (m+k1m1)

再分析有上界的情况。

i[1,k],niri 时,我们可以用容斥把约束干掉。依次考虑哪些约束被选中,然后删去这些约束后用无上界的情况解决,最后容斥原理即可。

公式太长了不放了

卡特兰数

原始定义:Catn=x+y=n1CatxCaty

递推公式:Catn=4n2n+1Catn1

通项公式:Catn=(2nn)(2nn+1)

以下问题的解均为卡特兰数:

  1. 大小为 n 的二叉树数量。
  2. n+2 边形的三角剖分数量。
  3. n 个元素进栈出栈的排列数量。
  4. (0,0) 走到 (n,n),每次只能向右或向上走,且不能越过 y=x 这条线的方案数。

如果想要更详细的卡特兰数讲解,可以看 pjw 的博客。这里不过多讲解。

错位排列

设满足长度为 n 且满足 i[1,n],pii 的排列数量为 Dn,则 Dn=(n1)(Dn1+Dn2),或者 Dn=nDn1+(1)n.

我们设 fn 为恰好 n 个位置错开,gn 为至多 n 个位置错开。易知 gn=n!。根据二项式反演可以得到 fn=i=0n(1)ni(ni)i!

斯特林数

第二类斯特林数

第二类斯特林数 {nm} 表示将 n 个不同元素划分成 m 个相同的非空集合的方案数。

递推公式为 {nm}={n1m1}+m{n1m}。表示枚举将最后加入的一个元素另开一个集合还是加入原来的集合。

通项公式为:{nm}=i=0m(1)miini!(mi)!。直接用下面的等式中第一个二项式反演一下即可。

mn=i=0m or n{ni}i!(mi)

可以类比十二重计数法中的情况,mn 表示将 n 个不同的小球放入 m 个不同的盒子中,我们可以枚举有 i 个盒子装了小球,斯特林数为装起来的方案数,再考虑盒子的排列 i!,同时再考虑选出哪些盒子。就能得到这个等式。这个等式在化简一些 xk 的式子中十分有用。

i=0nik=i=0nj=0k{kj}j!(ij)=j=0k{kj}j!i=0n(ij)=j=0k{kj}j!(n+1j+1)

我们也可以得到自然数幂求和的一个结论:

i=0nik=i=0nj=0k{kj}j!(ij)=j=0k{kj}j!i=0n(ij)=j=0k{kj}j!(n+1j+1)

第一类斯特林数

第一类斯特林数 [nm] 为将 n 个不同元素划分为 m 个不区分的环的方案数。(环的意思就是可以轮换,比如 1,2,3,4,55,1,2,3,4 就是一个环)

递推公式:[nm]=[n1m1]+(n1)[n1m]。枚举将一个新元素放入原有的环或者新开一个环。

第一类斯特林数没有实用的通项公式。

本来这里想咕掉不写斯特林的一些科技的,直到我发现 2020 年省选考了下降幂

上升幂与下降幂

下降幂:xn=x!(xn)!=i=0n1(xi)

上升幂:xn¯=i=0n1(x+i)

下降幂性质:mk(nm)=(nkmk)nk

证明:

mk(nm)=m!(mk)!n!m!(nm)!=(nk)!(mk)!(nm)!n!(nk)!=(nkmk)nk

下降幂上升幂互换:xn=(1)n(x)n¯

普通幂转下降幂:xn=i=0n{ni}xi。下降幂转普通幂:xn=i=0n(1)ni[ni]xi

普通幂转上升幂:xn=i=0n(1)ni{ni}xi¯。上升幂转下降幂:xn¯=i=0n[ni]xi


对第二个式子和第四个式子的证明:

第二个式子:

考虑数学归纳法

xn+1=(xn)xn=(xn)i=0n(1)ni[ni]xi=xi=0n(1)ni[ni]xini=0n(1)ni[ni]xi=i=0n(1)ni[ni]xi+1ni=0n+1(1)ni[ni]xi=i=1n+1(1)ni+1[ni1]xi+ni=1n+1(1)ni+1[ni]xi=i=1n+1([ni1]+n[ni])(1)ni+1xi=i=1n+1[n+1i](1)n+1ixi=i=0n+1[n+1i](1)n+1ixi

第四个式子可以把上升幂转下降幂,再用第二个式子的结论。

xn¯=(1)n(x)n=i=0n[ni](1)n(1)ni(x)i=i=0n[ni]xi


斯特林反演

gn=i=0n{ni}fi,则 fn=i=0n(1)ni[ni]gi

斯特林反演不太常用,可以暂时不掌握。

题目

组合数学不是会数数就行吗
我觉得组合数学它都不能叫做数学
你高中三年数学不学都会组合数学
------------------------- 白宇

基础组合数学

[SDOI2016] 排列计数

求有多少种 1n 的排列满足恰好有 m 个位置满足 ai=i

数据范围:T5×105,n,m1×106


原排列除去这 m 个位置就是一个错位排列,直接用递推公式求。而我们要从 n 个位置中选出 m 个位置不动,总方案数就为 Dnm(nm)。预处理即可。时间复杂度 O(T+n)

[NOI Online #2 入门组] 建设城市

求满足以下条件的长度为 2n 的序列 a 的个数:最大值不超过 m,前 n 个数单调不降,后 n 个数单调不增,且 ax=ay

数据范围:n2×105


我们可以枚举 x 位置的高度,假设当前枚举到的高度为 i,之后分情况讨论:

xny 时,整个序列被 x,y 分成了四段。我们只拿出 [1,x1] 这一段来讨论,其他段同理。我们可以把 [1,x1] 这些城市看成小球,每个高度看成不同的盒子。因为最后都要保持有序,我们不妨把小球都看成相同的。我们就转化成了一个经典模型,插板法即可,方案数为 (x1+i1i1)

x<y<nn<x<y 时,我们可以把 xy 这些建筑看成一个建筑,然后像上面的方法计算即可。

时间复杂度 O(n)

[HAOI2011]Problem c

n 个人,每个人有一个座位 ai。现在从 1n 依次入座,如果他要坐的地方已经有人了,他就会依次尝试 ai+1,ai+2,依此类推。如果一直尝试到 n 都没有座位,那么该方案就不合法。现在已经确定了 m 个人的座位,求合法方案数。

数据范围:n300


首先考虑无解的情况,我们设 sumi 表示已经确定的编号大于等于 i 的人数。当 sumi>ni+1 时无解。再考虑有解的情况。我们设 fi,j 为剩余编号大于等于 i 的人中已经确定了 j 个人的座位。我们枚举 k 表示当前选 k 个人的编号为 i。我们就有转移方程 fi,j=k=0j(jk)fi+1,jk(0jni+1sumi)。答案即为 f1,nm

[SHOI2015]超能粒子炮·改

i=0k(ni)mod2333


fn,k=i=0k(ni)mod2333

开始大力推式子

fn,k=i=0k(ni)modp=i=0k(n/pi/p)(nmodpimodp)=(n/p0)i=0p1(nmodpi)+(n/p1)i=0p1(nmodpi)++(n/pk/p1)i=0p1(nmodpi)+(n/pk/p)i=0kmodp(nmodpi)=i=0p1(nmodpi)×i=0k/p1(n/pi)+(n/pk/p)i=0kmodp(nmodpi)=fnmodp,p1×fn/p,k/p1+(n/pk/p)fnmodp,kmodp

我们就可以递推出前 2333 项的 f,然后用类似卢卡斯定理的方式求解即可。

[CQOI2018]交错序列

交错序列的定义是:只由 01 组成的一个序列,且序列中所有 1 都互不相邻。一个交错序列的权值为:设 0 的个数为 x 个,1 的个数为 y 个,那么权值即为 xayba,b 为给定的整数)。求长度为 n 的所有交错序列的权值和。

数据范围:n107


我们枚举 1 的个数,考虑有 x1 的交错序列的个数。我们可以先把 nx0 放入序列,再在 nx+1 个空隙中放入 x1,方案数即为 (nx+1x)。我们再线性筛出 iaib。时间复杂度 O(n)

Anton and School - 2

给定一个长度为 n 的由小括号组成的序列,求有多少子序列满足左半部分为 (,右半部分为 )

数据范围:n2×105


我们枚举每一个 (,在它左边有 a((包括它自己),右边有 b)。那么它对答案的贡献为 i=0min(a1,b1)(a1i)(bi+1)。我们发现这玩意和范德蒙德卷积 i=0k(ni)(mki)=(n+mk) 差不多。我们考虑化简这个式子。

我们把上面的 min 拆掉,那么当 ab

i=0a1(a1i)(bi+1)=i=0a1(a1a1i)(bi+1)=i=0a(a1a1i)(bi+1)=(a+b1a)

ba 时,

i=0b1(a1i)(bi+1)=i=0b1(a1i)(bb1i)=(a+b1b1)=(a+b1a)

因此对答案的贡献即为 (a+b1a)。直接 O(n) 即可

梦现时刻

fa,b=i=0b(bi)(nia),求 a=1mb=1mfa,bmod998244353

其中 为异或运算。

数据范围:1n100,m5000


比较牛逼的推式子题。

fa,b=i=0b(bi)(nia)=i=0b(b1i)(nia)+i=0b(b1i1)(nia)=fa,b1+i=0b1(b1i)(ni1a)=fa,b1+i=0b1(b1i)(nia)i=0b1(b1i)(ni1a1)=2fa,b1i=0b1(b1i)(ni1a1)=2fa,b1(i=0b1(bi+1)(ni1a1)i=0b1(b1i+1)(ni1a1))=2fa,b1i=0b(bi)(nia1)+i=0b(b1i)(nia1)=2fa,b1fa1,b+fa1,b1

直接 O(n2) 递推。

[AGC002F] Leftmost Ball

给你 n 种颜色的球,每种颜色的球有 k 个,把这 n×k 个球排成一排,把每一种颜色的最左边出现的球涂成白色(初始球不包含白色),求有多少种不同的颜色序列。

数据范围:1n,k2000


合法的情况为:在任意前缀中,白球的种类都大于等于其他颜色的种类。

fi,j 为已经放了 i 个白球和 j 种颜色的方案数。

为了防止方案重复,我们规定,每次放球时,必须将球放到第一个空位上。

考虑最左边的第一个空位放什么球。当放白球时,可以从 fi1,j 转移过来。

当放其他颜色的球时,可以从 fi,j+1 转移过来。我们考虑这种情况的转移:先从剩余的 nj+1 种颜色中选出一种,再让一个球放到最左边的空位中,剩下 k2 个球在剩下的空位中随便放。因此转移方程即为 fi,j=fi1,j+fi,j1×(nj+1)×(nki(j1)(k1)1k2)

时间复杂度为 O(n2)

容斥原理 二项式反演

硬币购物

有四种面值为 c1,c2,c3,c4T 次询问,每次询问给出 lim1,lim2,lim3,lim4,tot,求出在第 i 种硬币不超过 limi 的情况下凑出 tot 面值的方案数。

数据范围:T1000


先做一遍完全背包,求出没有限制时凑出面值 t 的方案。再用类似多重集组合数的思想,容斥一下把限制消掉,即可得出答案。时间复杂度 O(16T)

[CSP-S2019] Emiya 家今天的饭

n 种烹饪方法,m 种主要食材,能用第 i 种烹饪方法和第 j 种食材做出 ai,j 种菜。同时有 3 条限制:做至少一道菜,每道菜烹饪方法互不相同,每种食材至多在一半的菜中被使用。求合法方案数。

数据范围:n100,m2×103


考虑容斥,让一种菜超出一半的限制,再用总方案减去即可。考虑求出总方案。设 fi,j 为前 i 种菜选 j 种食材的方案数,sij=1maj。转移即为 fi,j=fi1,j+si×fi1,j1,总方案即为 i=1nfn,i。再减去不符合限制的。设当前枚举到的食材为 p,设 gi,j,k 表示前 i 个食材中选了 p 这种食材 j 次,其他食材 k 次的方案,根据选不选 p 食材有下面的转移:gi,j,k=gi1,j,k+ai,pgi1,j1,k+(siai,p)gi1,j,k1,但是这样复杂度为 O(mn3),考虑优化。注意到我们其实只关心 jk 的相对大小,那我们就设 gi,j 为前 i 个食材中选了 p 这种食材的次数再减去其他食材选的次数为 j 的方案。那么我们就有转移:gi,j=gi1,j+ai,pgi1,j1+(siai,p)gi1,j+1。时间复杂度 O(mn2)

[HNOI2011] 卡农

从集合 S={1,2,3n} 中选出 m 个非空子集,满足所有元素在这些集合中出现的总次数都为偶数,求方案数。

数据范围:n106


我们可以把题意转化为:有 2n1 个二进制数,求从其中选出 m 个数异或为 0 的方案数。我们考虑先求出选择的排列数,最后再除以 m! 即可。我们设 fi 为从这些数中选出 i 个数的排列使得异或为 0 的方案。总的方案数为 A2n1i1,也就是先选出 i1 个数,这些数异或为 x,再将 x 选上,再减去不合法的方案数。考虑有哪些情况是不合法的,第一种是当前选择 x=0 的情况,这时前 i1 个数已经构成了合法方案,需要减去。再考虑重复的问题,假设 x 与前面选择的数有重复,那么 x 和这个数异或后为 0,剩下 i2 个数也是合法的。x2n1(i2) 种取值,再考虑第几个数和 x 重复,方案就有 i1 种。因此总的递推式即为 fi=A2n1i1fi1fi2×(2ni+1)×(i1)。答案就为 fmm!

[HAOI2018]染色

有一个长度为 n 的序列,有 m 种颜色给这个序列染色。我们只关心出现次数恰好为 s 的颜色种类。当有 k 种颜色恰好出现了 s 次时,愉悦度为 wk。求所有方案的愉悦度之和模 1004535809 的值。


1004535809 为质数,原根为 3,常用作 NTT 模数。

我们设 fi 为恰好出现 s 次的颜色恰好有 i 种的方案数,gi 为恰好出现 s 次的颜色至少有 i 种的方案数。我们考虑先求出 gi。我们先从 m 种颜色中选出 i 种颜色,再考虑在序列中如何放置这些颜色。我们把空着的 nis 位置都看成相同的颜色,这就是一个多重集组合数,方案为 n!(s!)i(nis)!。再把剩下的 mi 种颜色放到剩下 nis 个位置上,方案数为 (mi)nis。因此 gi=(mi)n!(s!)i(nis)!(mi)nis

再考虑求出 fk。设 up=min(m,ns)

fk=i=kup(1)ik(ik)gi=i=kup(1)iki!(ik)!k!gi=1k!i=kup(1)ik(ik)!(i!gi)

F(x)=i=0up(1)ii!xi,G(x)=i=0upi!×gi×xi。则 fk=1k!i=kupFikGi。这就是一个常见的差卷积形式,翻转一下即可转变为和卷积形式,直接 NTT 即可。时间复杂度 O(mlogm)

[JSOI2011]分特产

m 种不同的特产,每种特产有 ai 个,要把这些特产分给 n 个人,要求每个人必须有至少一件特产,求方案数。

数据范围:n,m,ai1000


gi 为至少 i 个人没有特产的方案数,先选出没有拿到特产的人,再把剩下的特产每种分别分给剩下的人,直接插板法即可,则 gi=(ni)j=1m(aj+ni1ni1)。最后二项式反演 ans=i=0n(1)igi。时间复杂度 O(nm)

[TJOI2019]唱、跳、rap和篮球

现在有 aAbBcCdD,要求把这些字母中的 n 个填到长度为 n 的序列中,要求序列中不能按顺序出现相邻的 ABCD。求合法方案数对 998244353 取模的值。

数据范围:n1000,a,b,c,d500


考虑二项式反演,设 gi 为序列中至少出现了 iABCD 的方案,我们先把这 i 个位置填上,再在剩下 n4i 个位置多重集组合数下,就可以求出 gi。考虑写出 gi 的表达式,设 hn,a,b,c,d 为有 aAbBcCdD 的情况下在 n 个位置中随便放的方案数,则 gi=(n3ii)hn4i,ai,bi,ci,di。再考虑求出 h。根据多重集组合数,hn,a,b,c,d=p=0aq=0bs=0ct=0dn!p!q!s!t![p+q+s+t=n],这就是一个卷积形式了,直接 NTT 即可。时间复杂度 O(n2logn)

[ARC160D] Mahjong

统计满足经过以下两种操作可以使得元素全为 0 的长度为 n,元素和为 mA 序列的个数:1. 选取一个元素使其减去 k。2. 选取长度为 k 的子串,使其全部都减去 1

数据范围:k,n2000,m1018


有意思的计数题。注意到无解当且仅当 km。再考虑有解如何统计。我们可以把操作也看成一个序列 BB 序列的长度为 2nk+1B 序列的前 n 个数代表的是第一种操作。Bi=p 就意味着在 i 这个位置进行了 p1 操作。B 序列的后 nk+1 个数表示的是第二种操作,Bi+n=p 就意味着在以 i 开头的子串中进行了 p 次操作。同时为了防止重复,我们规定 Bi<k(n<i2nk+1)。这样我们就能建立起操作和原序列的一一对应关系了。我们再考虑 B 序列的计数。首先 i=12nk+1Bi=mk,其中只有后面的 Bi 是有限制的,这就转化成了多重集组合数问题,直接容斥。设 gi 为至少有 iBi 不满足限制的方案数,先选出哪些位置不满足要求,然后直接插板即可,因此 gi=(nk+1i)(mki×k+2nk2nk)。最后二项式反演即可得出答案。

[ABC214G] Three Permutations

给定长度为 n 的两个排列 p,q,求满足 i[1,n],ripi,riqir 排列的方案数。

数据范围:n3000


考虑二项式反演,求出至少有 i 个限制不满足的方案 hi。我们从 piqi 连边,这时我们把这条边称作第 i 条边,第 i 个限制不满足就等价于第 i 条边对应它的起点或终点。把这些边都连上之后会形成一些环,我们把没有对应自己起点或终点的边删去。我们先考虑一个环,设 dpi,j 为把长度为 i 的链删去 j 条边之后对应的方案数,易知大小为 i 的链对应的方案数为 i,则 dpn,m=i=1nidpni,m1。再考虑环,设 fi,j 为把大小为 i 的环删去 j 条边之后对应的方案数,则 fn,m=i=1ni2fni,m1。再考虑所有的环,设 gi,j 为前 i 个环共删去 j 条边的对应方案数,则 gn,m=i=1vngn1,mfvn,i。最后 hi 就等于 gcnt,ni

随机立方体

有一个 n×m×l 的立方体,立方体中每个格子上都有一个数,如果某个格子上的数比三维坐标至少有一维相同的其他格子上的数都要大的话,我们就称它是极大的。现在将 1n×m×ln×m×l 个数等概率随机填入 n×m×l 个格子(即任意数字出现在任意格子上的概率均相等),使得每个数恰出现一次,求恰有 k 个极大的数的概率。

数据范围:n,m,l5×106,k100


看到恰好,直接二项式反演。我们只需要考虑至少有 k 个极大的数的情况。我们把这种情况的方案数记作 f(k)。考虑如何求出 f(i)

我们钦定完这 i 个数后,剩下的 (ni)(mi)(li) 个位置都是可以随便放置的。然后我们再按顺序选择出来这 i 个坐标不相同的数。总的方案数我们可以写作:

fi=(nimili)((nml)(ni)(mi)(li))gi

第一项为选择这 i 个坐标不同的数,第二项为放置剩下的位置,第三项为被这些极大点所限制的数的填法。我们再考虑求出 gi

由于每个点可能会受到多个极大点的限制,因此我们需要从大到小来考虑填数,因为小的极大点的限制更强。我们结合下面这张图来推式子。

在图中,粉色部分为无限制区域,一半红的是放了关键点的区域,与关键点另一半相同颜色的区域是受该关键点控制的区域(这里为了方便,将关键点从左上到右下按从大到小排序)。这时我们可以计算一下:g(3)=297×219×1111

因此我们可以得到 gi 的公式:

h(x)=(nx)(mx)(lx),则

gi=j=1i(h(0)h(j)1)h(j1)h(j)1=j=1i(h(0)h(j)1)!(h(0)h(j1))!=(h(0)h(i)1)!j=1i11(h(0)h(j))!

再带回 fi 中,即可得到:

fi=(nimili)((nml)(ni)(mi)(li))gi=(nimili)((nml)(ni)(mi)(li))(h(0)h(i)1)!j=1i11(h(0)h(j))!=(nimili)(nml)!(nmlh(i))(nmlh(i)1)!j=1i11nmlh(j)=(nimili)(nml)!j=1i1nmlh(j)

由于最终还要除掉 (nml)!,则我们直接把它去掉。

f(i)=(nimili)j=1i1nmlh(j)ans=i=k(1)ik(ik)f(i)

Ban Permutation

问长度为 n,且 i,|ipi|X 的排列的数量。

数据范围:n100,X5


我们先用二项式反演将题目的限制变为 |ipi|<X。我们考虑钦定 i 个位置满足的方案。

我们设 gi,j,s 为用前 i 个数,已经钦定了 j 个数,并且 [iX+1,i+X1] 这段区间的占用情况为 s 时的方案数。转移就按照是否钦定当前数来转移。具体的说:

gi,j,T=gi1,j,s+gi1,j1,s

这样就能求出钦定 i 个数时的方案数,二项式反演即可得到结果。

时间复杂度 O(4Xn2)

Lamps and Buttons

n 盏灯,其中前 m 盏灯是亮着的。现在随机生成一个排列 p,第 i 个开关能被按下当且仅当第 i 盏灯亮着,按下第 i 个开关会使得第 pi 盏灯的状态改变。现在一个人来按开关,这个人不知道 p 排列。当所有灯亮起时他赢,如果无法操作时他输。假设这个人绝顶聪明,问对于所有的排列 p,这个人的胜率之和。

数据范围:n107,m5000


我们把 ipi 连边。我们先考虑什么情况会输:按到自环或者在没亮的灯中出现了环。这样无论如何也无法使得所有灯亮起来。由于这个人不知道排列 p,因此他最好的顺序就是没有顺序,从左往右按就完事了。我们的问题就转化为了:把灯从左往右按,在碰到第一个自环前把所有的灯按亮的方案数。

我们枚举第一个自环的位置 i,我们发现在它前面没有自环的方案数很难求,我们可以用二项式反演求出在它前面有 x 个自环的方案数,同时我们还要满足后面 nm 个数都能包在前面的环中。我们可以让后面这些数插入到前面 ix 个位置所成的环中,因此总方案数为(下面的公式中 a=ix1,b=nm,c=mi

1×2××a×a×(a+1)××(a+b1)×(a+b+1)×(a+b+c)=a(a+b+c)(a+b)

总时间复杂度 O(n+m2)

推式子

Team Work

给定 n,k,求:i=1n(ni)×ik

数据范围:1k5000,1n109


i=1n(ni)ik=i=0n(ni)j=0k{kj}(ij)j!=j=0k{kj}j!i=0n(ni)(ij)=j=0k{kj}j!i=0n(nj)(njij)=j=0k{kj}j!(nj)i=0n(njij)=j=0k{kj}n!(nj)!2nj

直接 O(n2) 递推斯特林数即可。

[HEOI2016/TJOI2016]求和

i=0nj=0i{ij}2j×j!

数据范围:n105


i=0nj=0i{ij}2j×j!=j=0n2jj!i=0n{ij}=j=0n2jj!i=0nk=0j(1)jkkik!(jk)!=j=0n2jj!k=0j(1)jkk!(jk)!i=0nki=j=0n2jj!k=0j(1)jk(jk)!(kn+11)(k!)(k1)

推出卷积形式直接 NTT 即可。

[国家集训队] Crash 的文明世界

给定一棵树,对于每个点 i,求出 j=1ndisti,jk

数据范围:n50000,k150


i=1ndistu,ik=i=1nj=0k{kj}j!(distu,ij)=j=0k{kj}j!i=1n(distu,ij)

因此我们只需要求出对于每个点的 i=1n(distu,ij),就可求出答案。考虑树上 dp,设 ansu,ju 子树内 (distu,ij) 之和,考虑 (nm)=(n1m)+(n1m1),那么 ansu,j=vsonuansv,j+ansv,j1。最后换根 dp 就可以求出所有点的答案了。

CF1278F Cards

m 张牌,其中有一张为王牌。现在执行 n 次洗牌后选出最上面的一张牌的操作。设 x 为第一张为王牌的次数,求 xk 的期望值。

数据范围:n,m<998244353,k5000


ans=i=0n(ni)(1m)i(m1m)niik=i=0n(ni)(1m)i(m1m)nij=0k(ij)j!{kj}=1mnj=0k{kj}j!i=0n(m1)ni(ni)(ij)=1mnj=0k{kj}j!i=0n(m1)ni(nj)(njij)=1mnj=0k{kj}j!(nj)i=0nj(nji)(m1)nij×1i=1mnj=0k{kj}j!(nj)mnj=j=0k{kj}n!(nj)!×1mj

O(n2) 递推即可。

Partitions

n 个物品,每个物品有一个权值 wi。一个集合的权值为该集合的大小再乘上其内部物品的权值和,即 W(S)=|S|wi。一个划分的权值为所有集合的权值和。求 n 个物品划分成 k 个集合的所有方案的权值和。

数据范围:n,k2×105


易知每个物品在所有划分方案中贡献的次数是相同的,记为 cnt,开始推式子。

cnt=s=1n(s×(n1s1)×{nsk1})=s=1ns(n1s1)i=0k1(1)i(k1i)nsi!(k1i)!=i=1k1(1)ii!(k1i)!s=1ns(n1s1)(k1i)ns

我们再单独考虑后面的和式。

有一个公式 m(nm)=n(n1m1)

tmp=s=1n(n1s1)(k1i)ns+s=1n(s1)(n1s1)(k1i)ns=s=1n(n1s1)(k1i)ns+(n1)s=1n(n2s2)(k1i)ns

然后二项式定理即可。

tmp=(ki)n1+(n1)(ki)n2=(ki)n2(ki+n1)

因此最终答案为 i=1k1(1)i(ki)n2(ki+n1)i!(k1i)!

[省选联考 2020 A 卷] 组合数问题

给定一个 m 次多项式 f(x)=a0+a1x++amxm,求 k=0nf(k)×xk×(nk)

数据范围:n109,m1000


本题要用新科技:下降幂。本来不想写这题的发现有不会的新科技才写的

我们先把 f(x) 转化为下降幂多项式 f(x)=b0+b1x+bmxm,再开始推式子。

i=0nf(i)xi(ni)=i=0nj=0mbj×ij(ni)xi=i=0nj=0mbj×nj(njij)xi=i=0nxij=0mbj×nj(njij)=j=0mbjnji=0n(njij)xi=j=0mbjnji=0nj(nji)xi+j=j=0mbjnjxji=0nj(nji)xi=j=0mbjnjxj(x+1)nj

再回顾一开始普通幂转下降幂的式子,则有

i=0maixi=i=0maij=0i{ij}xj=i=0mxij=im{ji}aj

因此 bj=j=im{ji}aj。直接 O(n2) 递推即可。

The End

组合数学博客也肝了一段时间了,从开始做到肝完博客也花了十天左右,以后可能还会加点有意思的组合数学题(也感觉自己数数能力好了点)。

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