【笔记】《离散数学》第十章 递推方程与生成函数

OP

《离散数学》指清华大学出版社《离散数学（第三版）》

本文主要整理自考试范围中书上的内容，添加了少许个人的理解，受限于个人水平，一些内容还没有完全搞清楚，敬请谅解；

同时也作为以后算法竞赛中参考的纸质资料，以解决相关问题；

如有疏漏，敬请评论区指正；

10.1 递推方程及其应用

10.1.2 常系数线性齐次递推方程的求解

定义 10.2：关于k阶常系数线性齐次递推方程的定义

设递推方程满足
$$\{\begin{array}{l}H(n)-a_{1}H(n-1)-a_{2}H(n-2)-...-a_{k}H(n-k)=f(n)\\ H(0)=b_0,H(1)=b_1,H(2)=b_2,...,H(k-1)=b_{k-1}\end{array}$$
其中 $a_1,a_2,...,a_k$ 为常数，$a_k\ne 0$ ，这个方程称为k阶常系数线性递推方程。
其中k阶为k项，常系数指系数为常数，线性为 $H(i)$ 的次数为1；
当 $f(n)=0$ 时称这个方程为齐次方程。

定义 10.3：关于特征方程与特征根

给定常系数线性齐次递推方程如下
$$\{\begin{array}{ll}H(n)-a_{1}H(n-1)-a_{2}H(n-2)-...-a_{k}H(n-k)=0& \\ H(0)=b_0,H(1)=b_1,H(2)=b_2,...,H(k-1)=b_{k-1}& \end{array}(10.2)$$
方程 $x^k-a_1{x}^{k-1}-...-a_k=0$ 称为该递推方程的特征方程，特征方程的根称为该方程的特征根。

定理 10.1：关于递推方程与特征根

设 $q$ 是非零复数，则 $q^n$ 是递推方程(10.2)的解（$H(n)=q^n$）当且仅当 $q$ 是它的特征根。

定理 10.2：关于递推方程与特征根

设 $h_1(n)$ 和 $h_2(n)$ 是递推方程(10.2)的解，$c_1,c_2$ 为任意常数，则 $c_1{h}_1(n)+c_2{h}_2(n)$ 也是这个递推方程的解。

定义 10.4：关于通解

若对于递推方程(10.2)的每个解 $h(n)$​ 都存在一组常数 $c_1^{‘},c_2^{‘},...,c_k^{‘}$

使得 $h(n)=c_1^{‘}{q}_1^n+c_2^{‘}{q}_2^n+...+c_k^{‘}{q}_k^n$​ 成立，则称 $c_1{q}_1^n+c_2{q}_2^n+...+c_k{q}_k^n$​ 为该递推方程的通解；

定理 10.3：无重特征根与通解

设 $q_1,q_2,...,q_k$ 是递推方程(10.2)不等的特征根，则 $H(n)=c_1{q}_1^n+c_2{q}_2^n+...+c_k{q}_k^n$ 为该递推方程的通解。

定理 10.4：有重特征根与通解

设 $q_1,q_2,...,q_t$ 是递推方程(10.2)不等的特征根，且 $q_i$ 的重数为 $e_i$ ，则 $H(n)={\sum }_{i=1}^t(c_{i1}+c_{i2}n+c_{i3}{n}^2+...+c_{i{e}_i}{n}^{e_i-1})q_i^n$ 为该递推方程的通解。

求解递推方程有以下步骤：

1. 构造特征方程
2. 求解特征根
3. 构造通解
4. 用初值求解待定系数

例题10.5

求解 Fibonacci 数列的递推方程。

递推方程是 $f_n-f_{n-1}-f_{n-2}=0$ ，初值是 $f_0=1,f_1=1$ ；
得到特征方程 $x^2-x-1=0$ ；
解得特征根 $\frac{1+\sqrt{5}}{2},\frac{1-\sqrt{5}}{2}$ ；
设通解为 $f_n=c_1(\frac{1+\sqrt{5}}{2}{)}^n+c_2(\frac{1-\sqrt{5}}{2}{)}^n$ ；
代入初值 $f_0=1,f_1=1$ ，有

$$\{\begin{array}{l}{c}_1+c_2=1\\ c_1(\frac{1+\sqrt{5}}{2})+c_2(\frac{1-\sqrt{5}}{2})=1\end{array}$$

解得$c_1=\frac{1}{\sqrt{5}}\frac{1+\sqrt{5}}{2},c_2=-\frac{1}{\sqrt{5}}\frac{1-\sqrt{5}}{2}$ ；
所以递推方程的通解为

$$f_n=\frac{1}{\sqrt{5}}(\frac{1+\sqrt{5}}{2}{)}^{n+1}-\frac{1}{\sqrt{5}}(\frac{1-\sqrt{5}}{2}{)}^{n+1}$$

例题10.7

求解以下递推方程

$$\{\begin{array}{l}H(n)+H(n-1)-3H(n-2)-5H(n-3)-2H(n-4)=0\\ H(0)=1,H(1)=0,H(2)=1,H(3)=2\end{array}$$

特征方程是 $x^4+x^3-3x^2-5x-2=0$ ；
特征根为 $-1，-1，-1，2$ ；
设通解为 $H(n)=(c_1+c_{2}n+c_3{n}^2)(-1{)}^n+c_4{2}^n$ ；
代入初值得

$$\{\begin{array}{l}{c}_1+c_4=1\\ -c_1-c_2-c_3+2c_4=0\\ c_1+2c_2+4c_3+4c_4=1\\ -c_1-3c_2-9c_3+8c_4=2\end{array}$$

解得 $c_1=\frac{7}{9},c_2=-\frac{1}{3},c_3=0,c_4=\frac{2}{9}$ ；
所以通解为
$$H(n)=\frac{7}{9}(-1{)}^n-\frac{1}{3}n(-1{)}^n+\frac{2}{9}{2}^n$$

10.1.3 常系数线性非齐次递推方程的求解

常系数线性非齐次递推方程的标准形是$$H(n)-a_{1}H(n-1)-...-a_{k}H(n-k)=f(n)(10.3)$$其中 $n⩾k,a_k\ne 0,f(n)\ne 0$ 。

定理 10.5：通解的结构

设 $\overline{H(n)}$ 是对应齐次方程(10.2)的通解，$H^{\ast }(n)$ 是方程（10.3）的一个特解，则
$$H(n)=\overline{H(n)}+H^{\ast }(n)$$
是递推方程（10.3）的通解。

有了此定理，求解常系数线性非齐次递推方程只需要另外求解出方程的特解：

可按以下一般顺序：

1. 构造特征方程
2. 求解特征根
3. 根据非齐次项的形式构造特解
4. 代入递推方程求解特解中的未知常数
5. 由确定的非齐次特解和待定系数的齐次部分的通解构造非齐次通解
6. 用初值求解待定系数

以下是构造特解的一些规律：

1.1 特征根不为1时，如果 $f(n)$ 为 n 的 t 次多项式，那么特解也为 n 的 t 次多项式

例题10.8
找出以下递推方程的特解：

$$a_n-2a_{n-1}=2n^2$$

设 $a_n^{\ast }=P_1{n}^2+P_{2}n+P_3$ ；
代入递推方程并整理，得

$$-P_1{n}^2+(4P_1-P_2)n+(-2P_1+2P_2-P_3)=2n^2$$

对应系数相等得到线性方程组

$$\{\begin{array}{l}-P_1=2\\ 4P_1-P_2=0\\ -2P_1+2P_2-P_3=0\end{array}$$

解得 $P_1=-2,P_2=-8,P_3=-12$ ；
所以 $a_n^{\ast }=-2n^2-8n-12$ ；

1.2 特征根为1时，如果 $f(n)$ 为 n 的 t 次多项式，特解是将 n 的 t 次多项式的每一项提高一次

例题10.9
找出以下递推方程的特解：
$$W(n)=W(n-1)+n-1$$

设 $W^{\ast }(n)=P_1{n}^2+P_{2}n$ ；
代入递推方程并整理，得

$$2P_{1}n+P_2-P_1=n-1$$

从而得到方程组

$$\{\begin{array}{l}2P_1=1\\ P_2-P_1=-1\end{array}$$

解得 $P_1=1/2,P_2=-1/2$ ；
所以 $W^{\ast }(n)=n^2/2-n/2$ ；

2.1 $f(n)$ 为指数函数 $A{\beta }^n$ ，这里A代表常数，若 $\beta$ 不是特征根，则特解为 $P{\beta }^n$ ，其中P为待定系数

例题10.11
找出以下递推方程的特解：

$$a_n=6a_{n-1}+8^{n-1}$$

设 $a_n^{\ast }=P{8}^n$ ；
代入递推方程得 $P{8}^n=6P{8}^{n-1}+8^{n-1}$ ；
消去 $8^{n-1}$ 项，解得 $P=1/2$ ；
所以特解为 $a_n^{\ast }=8^n/2$ ；

2.2 $f(n)$ 为指数函数 $A{\beta }^n$ ，这里A代表常数，若 $\beta$ 是e重特征根，则特解为 $P{n}^e{\beta }^n$

例题10.12
求递推方程 $H(n)-5H(n-1)+6H(n-2)=2^n$ 的特解；

2为特征根，因此令特解 $H^{\ast }(n)=Pn{2}^n$ ；
代入特征方程得 $Pn{2}^n-5P(n-1)2^{n-1}+6P(n-2)2^{n-2}=2^n$ ；
对应系数相等得方程 $-P/2=1$ ；
解得 P ，所以特解为 $H^{\ast }(n)=-n{2}^{n+1}$ ；

3. 特解的组合形式

例如：对于递推方程 $a_n+2a_{n-1}=n+3^n$ ；
我们可以设特解为 $a_n^{\ast }=P_1+P_{2}n+P_3{3}^n$ ；
代入原方程解得 $P_1=-1,P_2=-2,P_3=3$ ；
所以特解为 $a_n^{\ast }=-1-2n+3^{n+1}$ ；

10.2 生成函数及其应用

10.2.1 牛顿二项式定理及二项式系数

定义 10.5：关于牛顿二项式系数

设 r 为实数，n 为整数，引入形式符号
$$\left(\begin{array}{c}r\\ n\end{array}\right)=\{\begin{array}{ll}0& n<0\\ 1& n=0\\ \frac{r(r-1)...(r-n+1)}{n!}& n>0\end{array}$$
称为牛顿二项式系数

例如：

$$\left(\begin{array}{c}-2\\ 5\end{array}\right)=\frac{(-2)(-3)(-4)(-5)(-6)}{5!}=-6$$
$$\left(\begin{array}{c}1/2\\ 4\end{array}\right)=\frac{(\frac{1}{2})(\frac{1}{2}-1)(\frac{1}{2}-2)(\frac{1}{2}-3)}{4!}=-\frac{5}{128}$$
$$\left(\begin{array}{c}4\\ 3\end{array}\right)=\frac{(4)(3)(2)}{3!}=4$$

也就是说，分母为 $n!$ ，分子便有从 $r-0$ 开始的 n 项；

r 为自然数时，牛顿二项式系数就成为普通的二项式系数；

定理10.6：关于牛顿二项式定理

设 a 为实数，则对于一切实数 $x,y,|x/y|<1$ ，有
$$(x+y{)}^a=\sum _{n=0}^{\infty }\left(\begin{array}{c}a\\ n\end{array}\right)x^n{y}^{a-n}$$

特别地有推论：（下式 $|x|<1$ ）

$$\frac{1}{1-x}=1+x+x^2+...$$
$$\frac{1}{(1-x{)}^2}=\sum _{n=0}^{\infty }(n+1)x^n$$

10.2.2 生成函数的定义与性质

定义 10.6：关于生成函数

设序列 { a n } \{ a_n\} {an​} 构造形式幂级数
$$G(x)=a_0+a_{1}x+a_2{x}^2+...+a_n{x}^n+...$$
称 $G(x)$ 为序列 { a n } \{ a_n\} {an​} 的生成函数

我们除了定义之外还常用以下的函数展开式（不必太过关注定义域）：

$$\frac{1}{1-x}=\sum _{n=0}^{\infty }{x}^n$$
$$\frac{1}{1+x}=\sum _{n=0}^{\infty }(-1{)}^n{x}^n$$

由序列求生成函数

例题 10.20(1)
求序列 { a n } \{a_n\} {an​} 的生成函数，$a_n=7\ast 3^n$ ；

$$G(x)=\sum _{n=0}^{\infty }7\ast 3^n{x}^n=7\sum _{n=0}^{\infty }(3x{)}^n=\frac{7}{1-3x}$$

由生成函数求序列

例题 10.21
已知序列 { a n } \{a_n\} {an​} 的生成函数为 $G(x)=\frac{2+3x-6x^2}{1-2x}$ ，求 { a n } \{a_n\} {an​} 的通项；

$$G(x)=\frac{2}{1-2x}+3x=2\sum _{i=0}^{\infty }(2x{)}^n+3x$$

所以

$$a_n=\{\begin{array}{ll}{2}^{n+1}& n\ne 1\\ 2^2+3& n=1\end{array}$$

注：因为多出来的 $+3x$ 项一定是 $a_1$ 提供的，所以该项需要特判；

10.2.3 生成函数的应用

求多重集的r-组合数

S = { n 1 ⋅ a 1 , n 2 ⋅ a 2 , . . . , n k ⋅ a k } S=\{n_1\cdotp a_1,n_2\cdotp a_2,...,n_k\cdotp a_k\} S={n1​⋅a1​,n2​⋅a2​,...,nk​⋅ak​} 的r-组合数就是不定方程

$$\{\begin{array}{ll}{x}_1+x_2+...+x_k=r& \\ x_i⩽n_i& i=1,2,...,k\end{array}$$

的非负整数解的个数；

同时就是生成函数 $G(y)=(1+y+...+y^{n_1})(1+y+...+y^{n_2})...(1+y+...+y^{n_k})$ 的展开式中 $y^r$ 项的系数；

不定方程解的个数

基本的不定方程

$x_1+x_2+...x_k=r$ （x_i为自然数）解的个数为生成函数 $G(y)=(1+y+y^2+...{)}^k$ 中 $y^r$ 项的系数；
有 $N=\left(\begin{array}{c}k+r-1\\ r\end{array}\right)$ ；

带限制条件的不定方程

$x_1+x_2+...x_k=r,l_i⩽x_i⩽r_i$

生成函数 $G(y)=(y^{l_1}+y^{l_1+1}+...+y^{r_1})(y^{l_2}+y^{l_2+1}+...+y^{r_2})...(y^{l_k}+y^{l_k+1}+...+y^{r_k})$

带系数的不定方程

$p_1{x}_1+p_2{x}_2+...+p_k{x}_k=r$（x_i为自然数）
生成函数 $G(y)=(1+y^{p_1}+y^{2p_1}+...)(1+y^{p_2}+y^{2p_2}+...)...(1+y^{p_k}+y^{2p_k}+...)$

正整数无序拆分

基本模型：将 N 无序拆分成正整数 $a_1,a_2,...,a_n$ ；

不允许重复

$G(y)=(1+y^{a_1})(1+y^{a_2})...(1+y^{a_n})$ 中 $y^N$ 项的系数；

允许重复

$$G(y)=(1+y^{a_1}+y^{2a_1}+...)(1+y^{a_2}+y^{2a_2}+...)...(1+y^{a_n}+y^{2a_n}+...)=\frac{1}{(1-y^{a_1})(1-y^{a_2})...(1-y^{a_n})}$$

10.3 指数生成函数及其应用

指数生成函数的定义

定义 10.7：关于指数生成函数

设 { a n } \{a_n\} {an​} 为序列，称
$$G_e(x)=\sum _{n=0}^{\infty }{a}_n\frac{x^n}{n!}$$
为 { a n } \{a_n\} {an​} 的指数生成函数

例题 10.27
给定正整数 m ，$a_n=P(m,n)$ ，则 { a n } \{a_n\} {an​} 的生成函数为

$$G_e(x)=\sum _{n=0}^{\infty }P(m,n)\frac{x^n}{n!}=\sum _{n=0}^{\infty }\frac{m!}{n!(m-n)!}{x}^n=\sum _{n=0}^{\infty }\left(\begin{array}{c}m\\ n\end{array}\right)x^n=(1+x{)}^m$$

由此可见，对于给定 m ，$(1+x{)}^m$ 既是集合组合数序列 { C ( m , n ) } \{C(m,n)\} {C(m,n)} 的普通生成函数，也是集合排列数序列 { P ( m , n ) } \{P(m,n)\} {P(m,n)} 的指数生成函数；

指数生成函数的性质

设数列 { a n } , { b n } \{a_n\},\{b_n\} {an​},{bn​} 的指数生成函数分别为 $A_e(x)$ 和 $B_e(x)$ ，则

$$A_e(x)\cdot B_e(x)=\sum _{n=0}^{\infty }{c}_n\frac{x^n}{n!}$$

其中

$$c_n=\sum _{k=0}^n\left(\begin{array}{c}n\\ k\end{array}\right)a_k{b}_{n-k}$$

指数生成函数的应用

多重集排列计数

定理 10.8：指数生成函数求解多重集的排列问题

设 S = { n 1 ⋅ a 1 , n 2 ⋅ a 2 , . . . , n k ⋅ a k } S=\{n_1\cdotp a_1,n_2\cdotp a_2,...,n_k\cdotp a_k\} S={n1​⋅a1​,n2​⋅a2​,...,nk​⋅ak​} 为多重集，则 S 的 r-排列数的指数生成函数为
$$G_e(x)=f_{n_1}(x)f_{n_2}(x)...f_{n_k}(x)$$
其中
$$f_{n_i}(x)=1+x+\frac{x^2}{2!}+...+\frac{x^{n_i}}{n_i!},n=1,2,...,k$$

其中 $\frac{x^r}{r!}$ 项的系数即为所求；

例
由1, 2, 3, 4 组成的五位数中，要求1出现不超过2次，但不能不出现，2出现不超过1次，3出现可达3次，4出现偶数次. 求这样的五位数个数；

$$G_e(x)=(x+\frac{x^2}{2!})(1+x)(1+x+\frac{x^2}{2!}+\frac{x^3}{3!})(1+\frac{x^2}{2!}+\frac{x^4}{4!})$$$$=x+5\frac{x^2}{2!}+18\frac{x^3}{3!}+64\frac{x^4}{4!}+512\frac{x^5}{5!}+...$$

即答案为 512；

例
红、白、蓝涂色 $1\times n$ 的方格，要求偶数个为白色，问有多少方案？

设方案数为 $a_n$ ；

$$G_e(x)=(1+x+\frac{x^2}{2!}+...{)}^2(1+\frac{x^2}{2!}+\frac{x^4}{4!}+...)$$
$$=\frac{1}{2}(e^x+e^{-x})e^{2x}=\frac{1}{2}(e^{3x}+e^x)$$
$$=\frac{1}{2}(\sum _{n=0}^{\infty }\frac{(3x{)}^n}{n!}+\sum _{n=0}^{\infty }\frac{x^n}{n!})$$

故

$$a_n=\frac{3^n+1}{2}$$

Catalan数

Catalan数的定义

定义 10.8：关于Catalan数

一个凸 n+1边形，通过不相交于n+1 边形内部的对角线把 n+1 边形划分成三角形，划分方案个数记作h_n，称为Catalan数

Catalan数的递推方程

$$\{\begin{array}{l}{h}_n={\sum }_{k=1}^{n-1}{h}_k{h}_{n-k}\\ h_1=1\end{array}$$

解为

$$h_n=\frac{1}{n}\left(\begin{array}{c}2n-2\\ n-1\end{array}\right)$$

ED

\

版权声明：本文为CSDN博主「Tan_Yuu」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/Tan_Yuu/article/details/117930803