组合数学——生成函数

最新推荐文章于 2022-07-19 15:10:48 发布

爱寂寞的时光

最新推荐文章于 2022-07-19 15:10:48 发布

阅读量3.9k

点赞数 3

分类专栏： # 组合数学、概率论、博弈论与数理统计等基础数学文章标签：算法

本文链接：https://blog.csdn.net/jiahonghao2002/article/details/115468694

版权

组合数学、概率论、博弈论与数理统计等基础数学专栏收录该内容

40 篇文章 2 订阅

订阅专栏

生成函数

标题

生成函数

生成函数（generating function），又称母函数，是一种形式幂级数，其每一项的系数可以提供关于这个序列的信息。

生成函数有许多不同的种类，但大多可以表示为单一的形式：

$\sum_{n=0}^{\infty}a_{n}k_{n}(x)$

其中 $a_{n}$ 是数列的第 $n$ 项， $k_{n}(x)$ 被称为核函数，不同的核函数决定了不同类型的生成函数。生成函数既可以是有限的，也可以是无限的。

其中，根据核函数的不同，生成函数大致分为：

普通型生成函数： $k_{n}(x)=x^{n}$
指数型生成函数： $k_{n}(x)=\frac{x^{n}}{n!}$
狄利克雷生成函数： $k_{n}(x)=\frac{1}{x^{n}}$

其中，最常用的为普通型生成函数，本文围绕普通型生成函数展开讨论。

普通型生成函数

普通型生成函数定义为幂级数的形式：

$\sum_{n=0}^{\infty}a_{n}x^{n}$

通常假设数列以 $a_{0}$ 开始。

给出几个常见的生成函数：

有限项生成函数： $a=\{1,2,3\}$ ， $F(x)=1 + 2x+3x^{2}$
$1$ 序列生成函数： $a=\{1,1,1,\ldots\}$ ， $F(x)=\sum_{n=0}^{\infty}x^{n}$
$2$ 的幂次生成函数： $a=\{1,2,4,,8\ldots\}$ ， $F(x)=\sum_{n=0}^{\infty}2^{n}x^{n}$
奇数生成函数： $a=\{1,3,5,,7\ldots\}$ ， $F(x)=\sum_{n=0}^{\infty}(2n+1)x^{n}$

换句话说，如果 $a_{n}$ 有通项公式，那么幂级数的形式很容易写出。

基本运算

加法运算

考虑两个序列 $a, b$ 的普通生成函数，分别为 $F (x), G (x)$ 。那么有

$F(x)\pm G(x)=\sum_{n=0}^{\infty} (a_n\pm b_n)x^n$

因此 $F(x)\pm G(x)$ 是序列 $\langle a_n\pm b_n\rangle$ 的普通生成函数。

乘法运算（卷积）

考虑乘法运算，也就是卷积：

$F(x)G(x)=\sum_{n=0}^{\infty} x^n \sum_{i=0}^na_ib_{n-i}$

因此 $F (x) G (x)$ 是序列 $\langle \sum_{i=0}^n a_ib_{n-i} \rangle$ 的普通生成函数。

注意：序列 $\langle \sum_{i=0}^n a_ib_{n-i} \rangle$ 满足一般乘法，加法原理。这是以后卷积运算在排列组合中的内部原理。

幂级数的封闭

有时，我们做两个无穷幂级数的卷积的时候，对于无穷多项多项式做卷积是一件很困难的事，我们可以用幂级数的收敛值来代替幂级数做卷积运算。

这种方法称为幂级数的封闭。

例如 $\langle 1,1,1,\cdots\rangle$ 的普通生成函数 $F(x)=\sum_{n=0}^{\infty}x^n$ ，我们可以发现

$F (x) x + 1 = F (x)$

那么解这个方程得到

$F(x)=\frac{1}{1-x}$

这就是 $\sum_{n = 0}^{\infty}x^n$ 的封闭形式。

注意：上式等号仅在收敛域内成立。

考虑等比数列 $\langle 1,p,p^2,p^3,p^4,\cdots\rangle$ 的生成函数 $F(x)=\sum_{n\ge 0}p^nx^n$ ，有

$\\ F(x) =\frac{1}{1-px}$

等比数列的封闭形式与展开形式是常用的变换手段。

除了上述两个下列给出几种常见的封闭运算：

求导型

$a=\langle 1,2,3,4,\cdots \rangle$ 。

第三个（求导）：

$\begin{aligned}F(x)&=\sum_{n\ge 0}(n+1)x^n\\&=\sum_{n\ge 1}nx^{n-1}\\&=\sum_{n\ge 0}(x^n)'\\&=\left(\frac{1}{1-x}\right)'\\&=\frac{1}{(1-x)^2}\end{aligned}$

牛顿二项式定理型

$a_n=\binom{m}{n}$ （ $m$ 是常数， $n\ge 0$ ）。

$F(x)=\sum_{n\ge 0}\binom{m}{n}x^n=(1+x)^m$

倒数型

$a_n=\binom{m+n}{n}$ （ $m$ 是常数， $n\ge 0$ ）。

$F(x)=\sum_{n\ge 0}\binom{m+n}{n}x^n=\frac{1}{(1-x)^{m+1}}$

积分型

不常用。

幂级数的展开

幂级数在封闭形式下运算之后，得到的函数应该再次展开成幂级数的形式，才能获取到新的数列信息。

常用的展开方法就是封闭的逆运算，以及适当的换元，求导，积分，泰勒级数等等。

应用

承德汉堡：偶数个
可乐：0 个或 1 个
鸡腿：0 个，1 个或 2 个
蜜桃多：奇数个
鸡块：4 的倍数个
包子：0 个，1 个，2 个或 3 个
土豆片炒肉：不超过一个。
面包：3 的倍数个

每种食物都是以“个”为单位，只要总数加起来是 $n$ 就算一种方案。对于给出的 $n$ 你需要计算出方案数，对 $10007$ 取模。

这是一道经典的生成函数题。对于一种食物，我们可以设 $a_n$ 表示这种食物选 $n$ 个的方案数，并求出它的生成函数。而两种食物一共选 $n$ 个的方案数的生成函数，就是它们生成函数的卷积。多种食物选 $n$ 个的方案数的生成函数也是它们生成函数的卷积。

在理解了方案数可以用卷积表示以后，我们就可以构造生成函数（标号对应题目中食物的标号）：

$\displaystyle\sum_{n\ge 0}x^{2n}=\dfrac{1}{1-x^2}$ 。
$1 + x$ 。
$1+x+x^2=\dfrac{1-x^3}{1-x}$ 。
$\dfrac{x}{1-x^2}$ 。
$\displaystyle \sum_{n\ge 0}x^{4n}=\dfrac{1}{1-x^4}$ 。
$1+x+x^2+x^3=\dfrac{1-x^4}{1-x}$ 。
$1 + x$ 。
$\dfrac{1}{1-x^3}$ 。

那么全部乘起来，得到答案的生成函数：

$F(x)=\frac{(1+x)(1-x^3)x(1-x^4)(1+x)}{(1-x^2)(1-x)(1-x^2)(1-x^4)(1-x)(1-x^3)} =\frac{x}{(1-x)^4}$

然后将它转化为展开形式：

$F(x)=\sum_{n\ge 1}\binom{n+2}{n-1}x^n$

因此答案就是 $\dbinom{n+2}{n-1}=\dbinom{n+2}{3}$ 。

例题

CF1548C

求解 $\sum_{i = 1}^n \binom{3i}{x}$ 。

考虑构造一个系数是 $\sum_{i = 1}^n \binom{3i}{t}$ 关于 $x^t$ 的一个多项式生成函数。

固定 $3 i$ 那么 $\binom{3i}{t} = [x^t](1+x)^{3i}$ ，那么所求多项式为 $\sum_{i = 1}^n \binom{3i}{t} = \sum_{i = 1}^n[x^t](1+x)^{3i} = [x^t]\sum_{i = 1}^n(1+x)^{3i}$ 。

考虑对 $\sum_{i = 1}^n(1+x)^{3i}$ 进行化简，套用等比数列求和公式得到：

$\frac {(x+1)^{3n+3}-(x+1)^3}{(x+1)^3-1}$

考虑多项式除法即可在 $O (n)$ 时间内算出所有的系数。

总结

生成函数是连接离散数学和连续数学的桥梁，离散数学的加法、乘法原理对应连续数学中的卷积运算，使得复杂的加法、乘法原理变成了相对简单的卷积运算，引入卷积之后，连续数学所有的求导、极限、积分理论均可使用。还有加快多项式卷积的FFT也可以使用。

爱寂寞的时光

关注

3
点赞
踩
15

收藏

觉得还不错? 一键收藏
0
评论
组合数学——生成函数

生成函数生成函数（generating function），又称母函数，是一种形式幂级数，其每一项的系数可以提供关于这个序列的信息。生成函数有许多不同的种类，但大多可以表示为单一的形式：F(x)=∑n=0∞ankn(x)F(x) = \sum_{n=0}^{\infty}a_{n}k_{n}(x)F(x)=n=0∑∞ankn(x)其中ana_{n}an是数列的第nnn项，kn(x)k_{n}(x)kn(x)被称为核函数，不同的核函数决定了不同类型的生成函数。生成函数既可以是有限的，也可
复制链接

扫一扫

专栏目录