【笔记】生成函数/母函数在通项公式上的应用

前言

在看了好多篇博客，翻阅了几本书后，终于对生成函数有了一点点理解，还请各位和我一样刚入门的同学一起静下心来仔细思考，最好在草稿纸上演算一下

本文面向oier，并没有运用严格的数学语言

---

定义

对于序列$\{a_i\}$，它对应的生成函数为

$$G(x)=\sum_{i=0}^{+\infty}a_ix^i=a_0+a_1x+a_2x^2+a_3x^3+…$$

相应的，一个生成函数与一个序列之间应该是唯一对应的，即已知一个序列$\{a_i\}$可以得到与其对应的唯一的生成函数$G(x)$，而已知一个生成函数$G(x)$可以得到对应的唯一的序列$\{a_i\}$

而我们可以利用这个生成函数中多项式的性质解决一些问题，一般是比较复杂的组合问题

---

前置技能

在此之前先介绍一个式子：

$\frac 1{1-x}=1+x+x^2+x^3+…$

即$\frac 1{1-x}=\sum_{i=0}^{+\infty}x^i$

证明：

将分母乘过去：
$1=(1-x)\sum_{i=0}^{+\infty}x^i\\=\sum_{i=0}^{+\infty}x^i-x\sum_{i=0}^{+\infty}x^i\\=\sum_{i=0}^{+\infty}x^i-\sum_{i=1}^{+\infty}x^i\\=x^0=1$

这个式子在后面会有很大的用处：

$$\frac 1{1-x}=1+x+x^2+x^3+…$$
变形有：
$$\frac 1{1-\lambda x}=1+\lambda+\lambda^2x^2+\lambda^3x^3+…$$
$$\frac 1{1-\alpha x}+\frac 1{1-\beta x}=(\alpha+\beta)x+(\alpha^2+\beta^2)x^2+(\alpha^3+\beta^3)x^3+…$$

---

引入例题

我们先用生成函数证明一个定理：$\sum_{i=0}^kC_n^iC_m^{k-i}=C_{n+m}^k$

这明显可以用定义证明：在$n$个白球和$m$个黑球中取$k$个的方案数为(在$n$个白球中取$i$个并在$m$个黑球中取$k-i$个的方案数)对于$i\in[0,k]$的和

但这里用生成函数进行证明：

对于序列$\{C_n^i\}$的生成函数为$G_1(x)=C_n^0+C_n^1x+C_n^2x^2+…$

我们知道$(1+x)^n$的式子展开后对应的$x^i$项的系数为$C_n^i$（这一部分看Pascal三角后理解会更透彻一些）

即这个序列的生成函数等价于$G_1(x)=(1+x)^n$，同样对于序列$\{C_m^i\}$的生成函数为$G_2(x)=(1+x)^m$

原定理中等号左边可视作卷积形式，将等号左边视作一个整体，其生成函数为$G_1(x)\cdot G_2(x)=(1+x)^n\cdot(1+x)^m=(1+x)^{n+m}$

也就是说等号右边的生成函数为$G_3=(1+x)^{n+m}$

将这个生成函数展开：$G_3=(1+x)^{n+m}=C_{n+m}^0+C_{n+m}^1x+C_{n+m}^2x^2+…$

将其中各项系数提出来即为原序列：$\{C_{n+m}^k\}$

即可得证

---

思路

整理下思路，我们的目的是证明这个式子等号两边相等，则我们可以将$x=1,2,3,…$带进去得到序列$\{a_i\}$
求出式子左边形成数列$\{a_i\}$和右边形成数列$\{b_i\}$的生成函数$G_1$与$G_2$，看是否相等，由于数列与生成函数之间唯一对应，所以只要生成函数相等则原数列相等，即在$x$取任意值下等号左右两边都相等

建议看完这一段后理一下思路，再看看上下的例题结合理解一下

---

例题1

我们来解决另外一个问题：给定递推式$a_n=2a_{n-1}+1$求通项公式$a_n$（$a_1=1$）

这个问题在我之前写的一篇博客上有解法：点我

接下来考虑生成函数的解法

设这个序列的生成函数

$G(x)=a_1x+a_2x^2+…=\sum_{i=0}^{+\infty}a_ix^i$

转化（其中注意$a_0=0$）

$G(x)=a_1x+a_2x^2+…$
$=(2a_0+1)x+(2a_1+1)x^2+…$
$=(2a_0x+2a_1x^2+…)+(x+x^2+…)$
$=2x(a_1x+a_2x^2+…)+(x+x^2+…)$

记得之前证的式子：$1+x+x^2+x^3+…=\frac 1{1-x}$

即
$G(x)=2xG(x)+\frac 1{1-x}$
$G(x)=\frac 1{(1-x)(1-2x)}=\frac 1{1-2x}-\frac 1{1-x}$

即$G(x)=\frac 1{1-2x}-\frac 1{1-x}=(1+2x+2^2x^2+…)-(1+x+x^2+…)$

算出每一项系数：$a_n=2^n-1$，而这个系数所组成的序列即为所求原序列
得通项公式：$a_n=2^n-1$

---

例题2

再用这种方法求一下斐波那契数列通项公式（$a_n=a_{n-1}+a_{n-2},a_1=a_2=1$）

设生成函数$G(x)$

根据递推式得$G(x)-x^2-x=x[G(x)-x]+x^2G(x)$（如果看不懂的话想想$a_1=a_2=1$）

继续变形：
$G(x)=\frac x{1-x-x^2}$
$=\frac 1{\sqrt 5}\bigl [\frac 1{1-\frac {1+\sqrt 5}2x}-\frac 1{1-\frac {1-\sqrt 5}2x}\bigr ]$

为方便表示，设$\alpha =1-\frac {1+\sqrt 5}2,\beta = 1-\frac {1-\sqrt 5}2x$：
$=\frac 1{\sqrt 5}[(\alpha-\beta)x+(\alpha^2-\beta^2)x^2+…]$

算出每一项系数即最终通项：$a_n=\frac {\alpha^n-\beta^n}{\sqrt 5}$

注意到这个式子是在无穷定义下得到的通项公式，所以其中的$n$应该具有$+\infty$的性质

所以这个式子中由于$|\beta|=\bigl |1-\frac {1-\sqrt 5}2\bigr |\leq 1$，即$\beta^n=0$

最终得到通项公式：$a_n=\frac {\alpha^n}{\sqrt 5}=\frac {\sqrt 5}5 \bigl ( \frac {1+\sqrt 5}2 \bigr )^n$

---

总结

由上我们可以总结一下用生成函数由递推公式推得通项公式的步骤：

• 设数列生成函数为$G(x)=\sum_{i=0}^{+\infty}a_ix^i$
• 由递推式推得$G(x)$的数值表达式
• 将$G(x)$的数值表达式转换出形似$\frac 1{1-\lambda x}$的式子并用一开始证明的式子进行转换
• 去除在无穷定义下无用项
• 将$x^i$的系数提出即为数列的项$a_i$

---

留坑

生成函数还不能熟练应用，许多性质的前置技能未激活，需要循序渐进，在接下来的时间里慢慢补上