从组合数学开始的数学相关知识总结（适合复习用）

Belief_yfly

已于 2023-10-13 12:42:09 修改

阅读量349

点赞数 1

分类专栏：总结文章标签：组合数学总结

于 2023-07-15 16:00:02 首次发布

本文链接：https://blog.csdn.net/weixin_44346868/article/details/131738651

版权

总结专栏收录该内容

7 篇文章 1 订阅

订阅专栏

全排列

例子： $n$ 个数取 $m$ 个数有序排放

通项公式： $A_n^m(P_n^m)=n*(n-1)*(n-2)*···*(n-m+1)= \frac{n!}{(n-m)!}$

因此在写公式时也可以直接把 $A_n^n$ 写成 $n!$

组合数

例子： $n$ 个数取 $m$ 个数，杨辉三角

递推公式： $C_i^0=1,\space C_i^j=C_{i-1}^j+C_{i-1}^{j-1}$

通项公式： $C_n^m=\frac{n!}{m!(n-m)!}$

组合数是全排列去掉顺序，所以也有 $C_n^m=\frac{A_n^m}{m!}$

性质： $C_n^m=C_n^{n-m},\space C_n^0+C_n^1+C_n^2+···+C_n^n=2^n$

需要注意 $n < m$ 时， $C_n^m=0$

杨辉三角

最简单的杨辉三角与组合数的递推公式相同

组合数的性质在杨辉三角中体现为对称以及第 $i$ 行所有数的和为 $2^{i-1}$

不同的是杨辉三角可以带有系数 $(a, b)$ ，即将杨辉三角第二行赋值为 $a$ 和 $b$

而杨辉三角的第 $i$ 行就是 $(a, b)$ 的 $i - 1$ 次方

我们可以通过这种方法求 $ax+b)^{i-1}$ 展开后每一项前的系数

p.s.我记得这玩意好像也可以用矩阵快速幂求，但是我不会，而且可能也没啥用，就不拓展了

递推公式为第 $i$ 行第 $j$ 个数 $x_{i,j}=a*x_{i-1,j}+b*x_{i-1,j-1}$

卡特兰数

例子：括号匹配，走楼梯，出栈顺序，二叉树计数

小数据：1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796

递推公式： $h_0=h_1=1,\space h_n=\sum_{i=0}^{n-1}h_i*h_{n-i-1}=\frac{h_{n-1}*(4n-2)}{n+1}$

通项公式： $h_n=C_{2n}^n-C_{2n}^{n-1}=\frac{C_{2n}^{n}}{n + 1}$

对于通项公式的拓展（以走楼梯为例）：

走楼梯问题可以抽象的认为是从 $(0, 0)$ 走到 $(n, n)$ 并且不超过直线 $y = x$ （只能向右或向上走）

如图为从 $(0, 0)$ 走到 $(8, 8)$ 并且不超过直线 $y = x$ 的一种走法
在这里插入图片描述
而有时我们需要求的并不是从 $(0, 0)$ 走到 $(n, n)$ ，假定我们需要求的是从 $x_a,y_a)$ 走到 $x_b,y_b)$

再具体的，我们假定从 $(3, 0)$ 走到 $(6, 5)$
在这里插入图片描述
接下来以这个例子来解释做法，考虑容斥原理

如果不考虑不超过 $y = x$ ，那么从 $(3, 0)$ 走到 $(6, 5)$ 总共有 $C_{6-3+5-0}^{6-3}$ 种

不超过 $y = x$ 等价于不经过 $y = x + 1$ ，我们将点 $(6, 5)$ 对称过去得到点 $(4, 7)$
在这里插入图片描述
从 $(3, 0)$ 走到 $(4, 7)$ 的方案数与从 $(3, 0)$ 走到 $(6, 5)$ 的不合法的方案是相等

解释如图，从 $(3, 0)$ 走到 $(4, 7)$ 的路线超出 $y = x + 1$ 的部分对称回去即可走到 $(6, 5)$
在这里插入图片描述
从 $(3, 0)$ 走到 $(4, 7)$ 总共有 $C_{4-3+7-0}^{4-3}$ 种

因此从 $(3, 0)$ 走到 $(6, 5)$ 的合法的方案即为 $C_{6-3+5-0}^{6-3}-C_{4-3+7-0}^{4-3}$

$x_b,y_b)$ 关于 $y = x + 1$ 对称即为 $y_b-1,x_b+1)$ ，不合法方案数 $C_{y_b-1-x_a+x_b+1-y_a}^{y_b-1-x_a}=C_{x_b-x_a+y_b-y_a}^{y_b-x_a-1}$

由此可以推导得出，从 $x_a,y_a)$ 走到 $x_b,y_b)$ 的方案数为 $C_{x_b-x_a+y_b-y_a}^{x_b-x_a}-C_{x_b-x_a+y_b-y_a}^{y_b-x_a-1}$

例题：登山

Q：如果出现 $y_b-1<x_a$ 的情况怎么办？（可以证明不存在 $x_b+1<y_a$ 的情况）

A：可以发现这种情况不存在不合法方案，所以方案数即为 $C_{x_b-x_a+y_b-y_a}^{x_b-x_a}$

第一类斯特林数

定义：用 $S1_{n,m}$ 表示将 $n$ 个不同元素构成 $m$ 个圆排列的数目

递推公式及证明：

边界情况，显然 $S1_{0,0}=1,S1_{n,0}=0$

容易发现 $S1_{n,m}$ 可以从两个状态转移过来

一个是 $S1_{n-1,m-1}$ ，这种情况下第 $n$ 个元素单独构成第 $m$ 个圆，产生 $S1_{n-1,m-1}$ 的贡献

另一个是 $S1_{n-1,m}$ ，这种情况下第 $n$ 个元素可以放在前 $n - 1$ 个元素的任意一个的前面，产生 $S1_{n-1,m}*(n-1)$ 的贡献

性质：

$S1_{n,1}=(n-1)!,\space S1_{n,2}=(n-1)!*\sum_{i=1}^{n-1}\frac{1}{i}$

$S1_{n,n}=1,\space S1_{n,n-1}=C_n^2, \space S1_{n,n-2}=2*C_n^3+3*C_n^4$

$\sum_{i=0}^nS1_{n,i}=n!$

第二类斯特林数

定义：用 $S2_{n,m}$ 表示把 $n$ 个不同的元素划分到 $m$ 个集合的方案数（集合不能为空）

很容易联想到放球模型，描述为：将 $n$ 个不同的小球放入 $m$ 个相同的盒子中，盒子不能为空，有几种方案

证明过程与第一类类似，不再赘述

递推公式： $S2_{0,0}=1,S2_{n,0}=0,S2_{n,m}=S2_{n-1,m-1}+S2_{n-1,m}*m$

通项公式： $S2_{n,m}=\frac{1}{m!}\sum_{i=0}^m(-1)^i*C_m^i*(m-i)^n$

性质：

$S2_{n,1}=1,\space S2_{n,2}=2^{n-1}-1,\space S2_{n,3}=\frac{3^{n-1}+1}{2}-2^{n-1}$

$S2_{n,n}=1,\space S2_{n,n-1}=C_n^2, \space S2_{n,n-2}=C_n^3+3*C_n^4,\space S2_{n,n-3}=C_n^4 +10*C_n^5+15*C_n^6$

$\sum_{i=0}^nS2_{n,i}=B_n$ （ $B_n$ 为倍尔数，下面会简单介绍）

两类斯特林数之间的递推式和实际含义很类似，他们之间存在一个~~没啥用的~~ 互为转置的转化关系：

$\sum_{k=0}^nS1_{n,k}*S2_{k,m}=\sum_{k=0}^nS2_{n,k}*S1_{k,m}$

倍尔数（贝尔数）

~~听说倍尔数与诗词有着奇妙的联系，应用倍尔数可以算出诗词的各种押韵方式~~

扯远了，回归正题

定义：用 $B_n$ 表示将 $n$ 个不同元素划分成若干个不相交集合的方案数

小数据：1, 1, 2, 5, 15, 52, 203, 877, 4140

规律：如图为倍尔三角形
在这里插入图片描述
我们可以发现第一竖列和右边斜行都是倍尔数

它有两条规律：每排的最后一个数都是下一排的第一个数；其他任何一个数等于它左边相邻数加左边相邻数上面的一个数

放球模型

将 $n$ 个小球放入 $m$ 个盒子，根据小球是否相同，盒子是否相同，盒子能否为空三个条件可分为八类模型

① 小球相同，盒子相同，盒子不可以为空

考虑dp，用 $f_{i,j}$ 表示把 $i$ 个小球放入至多 $j$ 个盒子中的方案数

边界情况： $f_{0,j}=1$

$f_{i,j}$ 会由两种状态转移过来，一种是有空盒，那么就至少有一个为空，方案数等于 $f_{i,j-1}$

另一种是没有空盒，那么 $j$ 个盒子每个都必须有一个小球，方案数等于 $f_{i-j,j}$

转移方程： $f_{i,j}=f_{i,j-1}+f_{i-j,j}$

而题目要求是没有空盒，所以答案是 $f_{n,m}-f_{n,m-1}$ ，即 $f_{n-m,m}$

② 小球相同，盒子相同，盒子可以为空

与第一种一样，答案是 $f_{n,m}$

③ 小球相同，盒子不相同，盒子不可以为空

搁板法（插空法），要形成 $m$ 个盒子，就需要 $m - 1$ 块板， $n$ 个小球会有 $n - 1$ 个空给这些板放置

盒子不可以为空，因此每个空最多只能放一块板，所以答案是 $C_{n-1}^{m-1}$

④ 小球相同，盒子不相同，盒子可以为空

与第三种相似，但因为盒子可以为空，所以可以假设现在每个盒子都放了一个

那么总共就有 $n + m$ 个球，所以答案是 $C_{n+m-1}^{m-1}$

⑤ 小球不相同，盒子相同，盒子不可以为空

即第二类斯特林数，答案是 $S2_{n,m}$

⑥ 小球不相同，盒子相同，盒子可以为空

与第五种相似，因为盒子可以为空，所以统计 $n$ 个球放入 $1, 2, 3, ..., m$ 个盒子的方案数之和即可，答案是 $\sum_{i=1}^mS2_{n,i}$

⑦ 小球不相同，盒子不相同，盒子不可以为空

与第五种相似，只需要多处理一个盒子的排列问题，答案是 $m!*S2_{n,m}$

⑧ 小球不相同，盒子不相同，盒子可以为空

这种情况下各个小球放入盒子互为独立事件，每个事件有 $m$ 种情况，最后答案便是 $m^n$

补充：因为以上都只考虑了 $n\ge m$ 的情况，感觉不是很完善，再来补充一点

盒子不可以为空时， $n < m$ 的情况输出 $0$

盒子可以为空时，如果盒子相同， $n < m$ 的情况与 $n = m$ 时相同；如果盒子不相同，处理方式与 $n\ge m$ 时一样

这个其实蛮好想的，就不再证明啦

整数划分

将一个正整数 $n$ 划分成 $m$ 个正整数，与放球模型第一种很相似

但我们要求的往往不是这种形式

先来点定义：

若 $n=a_1+a_2+…+a_i$ ，则称 ${a_1,a_2,...,a_m\}$ 为 $n$ 的一个划分

若 $max\{a_1,a_2,...,a_m\}\le m$ ，则称其为 $n$ 的一个 $m$ 划分

记 $n$ 的 $m$ 划分个数为 $f_{n,m}$

我们往往要求的是 $n$ 的所有划分个数，即 $f_{n,n}$

我们考虑如何求 $f_{n,m}$ （保证 $n, m > 0$ ）

边界情况： $f_{n,1}=f_{1,m}=1$

$n\ge m$ 时，可以从两种状态转移过来，一种是划分中没有 $m$ ，那么方案数就是 $f_{n,m-1}$

另一种是划分中有 $m$ ，可以假设 $a_m=m$ ，那么显然 $a_1+a_2+…+a_{m-1}=n-m$ ，可能还有 $m$ ，那么方案数就是 $f_{n-m,m}$

因此 $f_{n,m}=f_{n,m-1}+f_{n-m,m}$

很凑巧，它和放球模型第一种讲述的dp递推式一样，可以浅记一下

逆元

当组合数过大时不适宜 $O(n^2)$ 递推求解，这时我们就需要想另外办法处理

考虑到因为数字过大，很有可能会取模，那么就有可能会用到逆元

这里有一点可以记一下，逆元需要模数是质数，模数非质数大概率需要用递推求解

定义：当 $ax\equiv1(mod\space b)$ ， $x$ 即为 $a$ 在 $mod\space b$ 意义下的逆元

特殊的， $0$ 的逆元为 $1$

逆元的数学符号为inv

由定义可以发现在取模的意义下除以一个数就等于乘这个数的逆元

因此对于 $C_n^m=\frac{n!}{m!(n-m)!}$ 可以 $O (n)$ 预处理阶乘和其逆元然后 $O (1)$ 求解组合数

逆元一般有费马小定理和扩展欧几里德两种求法，这里仅介绍费马小定理，扩展欧几里德留到之后

费马小定理提出如果 $p$ 是一个质数，而整数 $a$ 不是 $p$ 的倍数，则有 $a^{p-1}\equiv1(mod\space p)$

因此 $p$ 一般是大质数比如 $998244353$ 和 $1 e 9 + 7$

不难看出 $a$ 和 $p$ 是互质的，所以式子两侧可以同除以 $a$ 得到 $a^{p-2}\equiv\frac{1}{a}(mod\space p)$

这个可以用快速幂 $O (l o g n)$ 求出

可是每个都求逆元的话时间复杂度是 $O (n l o g n)$ ，容易超时，考虑优化

$n!$ 的逆元可以看成 $\frac{1}{n!}$ ，因此乘 $n$ 就可以求出 $(n - 1)!$ 的逆元

因此也只需要 $O (n)$ 就可以求出所有阶乘的逆元了

扩展欧几里德

先引入一个常理，给予二整数 $a$ 与 $b$ , 必存在有整数 $x$ 与 $y$ 使得 $ax+by=gcd(a,\space b)$

扩欧是在欧几里德算法（辗转相除法）的基础上进行的拓展，目的就是在求出 $gcd(a,\space b)$ 的同时找到 $x$ 和 $y$

可以将代码当做板子记

int exgcd(int a, int b, int &x, int &y) {
    if (b == 0) {
    	x = 1; y = 0;
    	return a;
    }
    int res = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return res;
}

扩欧可以用来求逆元

假设 $a,\space x$ 在 $mod\space b$ 意义下互为逆元，即 $ax\equiv1(mod\space b)$

这就相当于可以找到一个 $y$ 使得 $a x - 1 = b y$ ，即 $a x - b y = 1$

而因为 $y$ 是自己取的可以直接取负，就可以化成经典的扩欧形式 $a x + b y = 1$

用扩欧求出来的 $x$ 就是 $a$ 在 $mod\space b$ 意义下的逆元

卢卡斯定理

卢卡斯定理一般用于解决组合数 $C_n^m$ 中 $n$ 和 $m$ 很大并且模数 $p$ 不大且为质数的问题

公式： $C_n^m=C_{n/p}^{m/p}*C_{n\space mod\space p}^{m\space mod\space p}(mod\space p)$

因为证明比较麻烦，不作证明，仅当作结论记

一般可以单独写出函数

int Lucas(int n, int m) {return n < MOD && m < MOD ? C(n, m) : Mul(Lucas(n / MOD, m / MOD), C(n % MOD, m % MOD));}

斐波那契

接下来的内容可能跟组合数学关系不大，仅做类比整理

正如组合数学中出现了很多经典的递推式，斐波那契也是一个非常经典的递推式

这里需要引出一个经典的问题：

一般而言，兔子在出生两个月后，就有繁殖能力，一对兔子每个月能生出一对小兔子来

如果我们在第一个月买回来一对幼兔，假设兔子不会死亡，那么一年后会有多少对兔子

列举每个月的兔子数量为 $1,\space1,\space2,\space3,\space5,\space8,\space13,\space21,\space34,\space55,\space89,\space144$

可以发现此数列中的某一项恰好会等于前两项相加，由此引导出了斐波那契的公式

递推公式： $f_1=f_2=1,\space f_i=f_{i-1}+f_{i-2}$

通项公式： $f_i=\frac{1}{\sqrt5}[(\frac{1+\sqrt5}{2})^n-(\frac{1-\sqrt5}{2})^n]$

斐波那契数列还有两个有趣的性质：

① 斐波那契数列中任一项的平方数都等于跟它相邻的前后两项的乘积相差 $1$

② 任取相邻的四个斐波那契数，中间两数之积（内积）与两边两数之积（外积）相差 $1$

由斐波那契序列中的质数组成的序列叫斐波那契质数序列

小数据： $2,\space3,\space5,\space13,\space89,\space233,\space1597,\space28657,\space514229,\space433494437,\space2971215073$

自然对数

求自然对数的过程中可以产生一个序列 ${a_n}=\sum_{i=1}^n\frac{1}{i}$

小数据： $1,\space1.5,\space1.83333,\space2.08333,\space2.28333,\space2.45,\space2.30714,...\space2.71828$

打算补充：概率期望，生成函数，多项式