【组合数学：四】生成函数 | 指数生成函数 | 线性齐次递推关系（上）

斐波那契数

• 之前博客整理过了 【小组专题三：斐波那契专题】

生成函数

• 一方面，我们可以把生成函数看成是代数对象，其形式上的处理使得人们可以通过代数手段计算一个问题的可能性的数目；另一方面，生成函数是无限可微分函数的泰勒级数。
  通常我们默认级数是收敛的且只在其形式上操作幂级数。
• 设：
  $$h_0,h_1,h_2,\cdots ,h_n,\cdots$$
  是无穷数列。它的生成函数定义为无穷级数
  $$g(x)=\sum _{i=0}^{\infty }{h}_i{x}^i$$
  每一个有限数列都有一个生成函数（后面的所有项都为 $0$），它是一个多项式。
• 例子：每一项都等于 $1$ 的无穷数列
  $$1,1,1,\cdots ,1,\cdots$$
  的生成函数是
  $$g(x)=\sum _{i=0}^{\infty }{x}^i=\frac{1}{1-x}$$
• 例子：设 $m$ 是正整数。二项式系数
  $$C_{m,0},C_{m,1},C_{m,2},\cdots ,C_{m,m}$$
  的生成函数是
  $$g_m(x)=\sum _{i=0}^m{C}_m^i{x}^i=(1+x{)}^m$$
• 例子：设 $k$ 是整数，并设数列 $h_n$ 等于
  $$e_1+e_2+\cdots +e_k=n$$
  的非负整数解的个数。
  它的生成函数是
  $$\begin{array}{rl}g(x)& =(\sum _{e_1=0}^{\infty }{x}^{e_1})(\sum _{e_2=0}^{\infty }{x}^{e_2})\cdots (\sum _{e_k=0}^{\infty }{x}^{e_1})\\ & =\frac{1}{(1-x)}\times \frac{1}{(1-x)}\times \cdots \times \frac{1}{(1-x)}\\ & =\frac{1}{(1-x{)}^k}\end{array}$$
  根据之前我们可知
  $$g(x)=\sum _{n=0}^{\infty }{C}_{n+k-1}^n{x}^n$$
• 例子：求装有苹果、香蕉、橘子和梨子的果篮的数量 $h_n$，其中在每个果篮中苹果的个数是偶数，香蕉数是 $5$ 的倍数，橘子数最多是 $4$ 个，梨子要么一个要么没有
  解：直接写出生成函数
  $$\begin{array}{rl}g(x)& =(1+x^2+x^4+\cdots )(1+x^5+x^{10}+\cdots )(1+x+x^2+x^3+x^4)(1+x)\\ & =\frac{1}{1-x^2}\frac{1}{1-x^5}\frac{1-x^5}{1-x}(1+x)\\ & =\frac{1}{(1-x{)}^2}\\ & =\sum _{n=0}^{\infty }{C}_{n+1}^n{x}^n\end{array}$$
  于是我们可以看到 $h_n=n+1$

指数生成函数

• 之前的生成函数特别适合涉及二项式系数的数列；然而对于某些计数排序的项数列，更有效的方法是考虑关于下面单项式集合的生成函数
  { 1 , x , x 2 2 ! , ⋯   , x n n ! , ⋯   } \{1,x,\frac{x^2}{2!},\cdots,\frac{x^n}{n!},\cdots\} {1,x,2!x2​,⋯,n!xn​,⋯}
  这些单项式出现在泰勒级数
  $$e^x=\sum _{n=0}^{\infty }\frac{x^n}{n!}$$
  中。数列 $h_0,h_1,\cdots ,h_n,\cdots$ 的指数生成函数定义为
  $$g^{(e)}(x)=\sum _{n=0}{h}_n\frac{x^n}{n!}$$
  后文简记为 $g(x)$
• 例子：设 $n$ 是正整数。确定下面数列的指数生成函数
  $$P(n,0),P(n,1),\cdots ,P(n,n)$$
  解：直接带入
  $$g(x)=\sum _{i=0}^{n}P(n,i)\frac{x^i}{i!}=\sum _{i=0}^n{C}_n^i{x}^i=(1+x{)}^n$$
• 例子：数列 { 1 , 1 , 1 , ⋯   , 1 , ⋯   } \{1,1,1,\cdots,1,\cdots\} {1,1,1,⋯,1,⋯} 的指数生成函数是 $g(x)=e^x$
  更一般的，数列 { a 0 , a 1 , a 2 , ⋯   , a n , ⋯   } \{a^0,a^1,a^2,\cdots,a^n,\cdots\} {a0,a1,a2,⋯,an,⋯} 的指数生成函数是 $g(x)=e^{ax}$
• 定理：设 $S$ 是多重集合 { n 1 ⋅ a 1 , n 2 ⋅ a 2 , ⋯   , n k ⋅ a k } \{n_1\cdot a_1,n_2\cdot a_2,\cdots,n_k\cdot a_k\} {n1​⋅a1​,n2​⋅a2​,⋯,nk​⋅ak​} ，其中 $n_i$ 是非负整数。设 $h_n$ 是 $S$ 的 $n$ 排列数。那么数列 $h_0,h_1,\cdots ,h_n,\cdots$ 的指数生成函数 $g(x)$ 为：
  $$g(x)=\prod f_{n_i}(x)$$
  其中，$\forall i\in [1,k]$，有：
  $$f_{n_i}(x)=1+x+\frac{x^2}{2!}+\cdots +\frac{x^{n_i}}{n_i!}$$
  证明略。
  注意到当 $n_i=\infty$ 时，上面的定理仍然成立。
• 例子：用红、白、蓝三种颜色给 $1\times n$ 的棋盘着色，要求红色格子是偶数个的方案数。
  解：直接写出相应的指数生成函数：
  $$\begin{array}{rl}g(x)& =(1+\frac{x^2}{2!}+\frac{x^4}{4!}+\cdots )(1+\frac{x^1}{1!}+\frac{x^2}{2!}+\cdots )(1+\frac{x^1}{1!}+\frac{x^2}{2!}+\cdots )\\ & =\frac{e^x+e^{-x}}{2}{e}^x{e}^x\\ & =\frac{e^{3x}+e^x}{2}\\ & =\frac{1}{2}(\sum 3^i\frac{x^i}{i!})(\sum \frac{x^i}{i!})\\ & =\frac{1}{2}\sum (3^i+1)\frac{x^i}{i!}\end{array}$$
  因此，我们能得到 $h_n=\frac{1}{2}(3^n+1)$
• 例子：一个 $n$ 位数，每个数字都是奇数且数字 $1$ 和数字 $3$ 出现偶数次的个数 $h_n$
  解：直接写出相应的指数生成函数：
  $$\begin{array}{rl}g(x)& =(1+\frac{x^2}{2!}+\frac{x^4}{4!}+\cdots {)}^2(1+x+\frac{x^2}{2!}+\cdots {)}^3\\ & =(\frac{e^x+e^{-x}}{2}{)}^2\cdot e^{3x}\\ & =\frac{1}{4}(e^{5x}+2e^{3x}+e^x)\\ & =\frac{1}{4}(\sum (5^i+2\cdot 3^i+1)\frac{x^i}{i!})\end{array}$$
  因此，我们得到 $h_n=\frac{1}{4}(5^n+2\cdot 3^n+1)$
• 例子：用红、白、蓝三色给 $1\times n$ 的棋盘着色，其中红色格子有偶数个，蓝色格子至少一个的方案数
  解：直接写指数生成函数：
  $$\begin{array}{rl}g(x)& =(1+\frac{x^2}{2!}+\frac{x^4}{4!}+\cdots )(1+x+\frac{x^2}{2!}+\cdots )(x+\frac{x^2}{2!}+\cdots )\\ & =\frac{e^x+e^{-x}}{2}{e}^x\times (e^x-1)\\ & =\frac{e^{3x}-e^{2x}+e^x-1}{2}\\ & =-\frac{1}{2}+\sum \frac{3^i-2^i+1}{2}\frac{x^i}{i!}\end{array}$$
  注意到，$h_0=-\frac{1}{2}+\frac{1}{2}=0$
  $h_n=\frac{3^n-2^n+1}{2},n>0$

求解线性齐次递推关系 （上）

• 设 $h_0,h_1,\cdots ,h_n,\cdots$ 是一个数列。称这个数列满足 $k$ 阶线性递推关系是指存在量 $a_1,a_2,\cdots ,a_k(a_k\ne 0)$ 和量 $b_n$ ，使得
  $$h_n=a_1{h}_{n-1}+a_2{h}_{n-2}+\cdots +a_k{h}_{n-k}+b_n(n\ge k)$$
  其中，$a_i，b_i$ 或是常数，或是依赖于 $n$ 的量。
• 我们称线性递推关系是 齐次的 ($homogeneous$) ,如果 $b_n$ 是常数 $0$。
  称它是 常系数 ($constantcoefficient$) 的，如果量 $a_1,a_2,\cdots ,a_k$ 是常数。
  我们后文讨论的是常系数线性齐次递推关系，即形如
  $$h_n-a_1{h}_{n-1}-a_2{h}_{n-2}-\cdots -a_k{h}_{n-k}=0(n\ge k)$$
  的多项式。
• 定理：如果多项式方程有 $k$ 个不同的根 $q_1,q_2,\cdots ,q_k$，则
  $$h_n=\sum _{i=1}^k{c}_i{q}_i^n$$
  无论给定什么样的初始值 $h_0,h_1,\cdots ,h_{k-1}$，都存在常数 $c_1,c_2,\cdots ,c_k$，使得上述式子满足递推关系和初始条件的唯一数列。
• 例子：初始值 $h_0=1,h_1=2,h_2=0$，求解下面的递推关系：
  $$h_n=2h_{n-1}+h_{n-2}-2h_{n-3}(n\ge 3)$$
  解：这个递推关系的特征方程是：
  $$x^3-2x^2-x+2=0$$
  它的三个根是 $1,-1,2$。根据定理，
  $$h_n=c_1{1}^n+c_2(-1{)}^n+c_3{2}^n$$
  是通解。现在，我们要求常数 $c_1,c_2,c_3$，使得他们满足 $h_0,h_1,h_2$ 是成立的
  利用消元法可以解出解是：$c_1=2,c_2=-\frac{2}{3},c_3=-\frac{1}{3}$。因此，
  $$h_n=2-\frac{2}{3}(-1{)}^n-\frac{1}{3}{2}^n$$
• 补充的内容：
  $$\frac{1}{(1-rx{)}^n}=\sum _{k=0}^{\infty }{C}_{n+k-1}^k{r}^k{x}^k$$
  例子：确定下面数列的生成函数
  $$0,1,4,\cdots ,n^2,\cdots$$
  解：带入上述例子，$r=1,n=2$，得到：
  $$\begin{array}{rl}\frac{1}{(1-x{)}^2}& =1+2x+3x^2+\cdots +n{x}^{n-1}+\cdots \\ \frac{x}{(1-x{)}^2}& =x+2x^2+3x^3+\cdots +n{x}^n+\cdots \\ \frac{1+x}{(1-x{)}^3}& =1+2^{2}x+3^2{x}^2+\cdots +n^2{x}^{n-1}+\cdots \\ \frac{x(1+x)}{(1-x{)}^3}& =x+2^2{x}^2+3^2{x}^3+\cdots +n^2{x}^n+\cdots \end{array}$$
• 例子：求解下面的递推关系
  $$h_n=5h_{n-1}-6h_{n-2}$$
  初始值 $h_0=1,h_1=2$
  我们把递推关系写成下面这样的形式：
  $$h_n-5h_{n-1}+6h_{n-2}=0(n\ge 2)$$
  设 $g(x)=h_0+h_{1}x+h_2{x}^2+\cdots +h_n{x}^n+\cdots$ 是数列 $h_0,h_1,h_2,\cdots ,h_n,\cdots$ 的生成函数。
  则有下面的一些方程，参照原来的递推关系，分别用 $-5x,6x^2$ 作为乘数乘以 $g(x)$ 得到：
  $$\begin{array}{rl}g(x)& =h_0+h_{1}x+h_2{x}^2+\cdots +h_n{x}^n+\cdots \\ -5xg(x)& =-5h_{0}x-5g_{1}x-\cdots -5h_{n-1}{x}^n-\cdots \\ 6x^{2}g(x)& =6h_0{x}^2+\cdots +6h_{n-2}{x}^n+\cdots \end{array}$$
  把上面三个方程加起来，我们得到
  $$(1-5x+6x^2)g(x)=h_0+(h_1-5h_0)x+(h_2-5h_1+6h_0)x^2+\cdots +(h_n-5h_{n-1}+6h_{n-2})x^n+\cdots$$
  因为 $h_n-5h_{n-1}+6h_{n-2}=0(n\ge 2)$
  又因为 $h_0=1,h_1=-2$，所以有：
  $$(1-5x+6x^2)g(x)=h_0+(h_1-5h_0)x=1-7x$$
  因此：
  $$g(x)=\frac{1-7x}{1-5x+6x^2}=\frac{c_1}{1-2x}+\frac{c_2}{1-3x}=\frac{5}{1-2x}-\frac{4}{1-3x}$$
  因此我们又可以得到：
  $$g(x)=\sum _{i=0}^{\infty }(5\cdot 2^i-4\cdot 3^i)x^i$$
  所以得到 $h_n=5\cdot 2^n-4\cdot 3^n$
  不过前提是特征方程的根是互不相同的，定理才成立。