【组合数学：五】线性齐次递推关系（下）| 非齐次递推关系

线性齐次递推关系（下）

• 例子：递推关系
  $$h_n=4h_{n-1}-4h_{n-2}$$
  有特征方程
  $$x^2-4x+4=(x-2{)}^2=0$$
  于是 $2$ 是二重特征根。在这种情况下，变成
  $$h_n=c_1{2}^n+c_2{2}^n=c{2}^n$$
  其中 $c=c_1+c_2$ 是一个新常数。但这不总是可以做到的。假如初始值是 $h_0=1,h_1=3$，那么容易得到 $c$ 是无解的。
  因此，$h_n=c{2}^n$ 不是给定递推关系的通解。
  某个特征根是重根，我们就需要寻找与这个根相关的其他解法。
  下面证明 $n{2}^n$ 也是一个解。
  $$h_n-4h_{n-1}+4h_{n-2}=n{2}^n-4(n-1)2^{n-1}+4(n-2)2^{n-2}=0$$
  我们现在断言：
  $$h_n=c_1{2}^n+c_{2}n{2}^n$$
  更一般的，如果 $q$ (可能是复数) 是常系数线性齐次递推关系的特征方程的 $s$ 重根，那么可以证明
  $$h_n=c_1{q}^n+c_{2}n{q}^n+c_3{n}^2{q}^n+\cdots +c_s{n}^{s-1}{q}^n$$
  对常数 $c_1,\cdots ,c_s$ 的每一种选择也是一个解
• 定理：设 $q_1,\cdots ,q_t$ 为常系数线性齐次递推关系
  $$h_n=a_1{h}_{n-1}+a_2{h}_{n-2}+\cdots +a_k{h}_{n-k},a_k\ne 0(n\ge k)$$
  的特征方程的互不相同的根。如果 $q_i$ 是特征方程的 $s_i$ 重根，那么这个递推关系中对应于 $q_i$ 的部分是：
  $$H_n^{(i)}=(c_1+c_{2}n+\cdots +c_{s_i}{n}^{s_i-1})q_i^n$$
  这个递推关系的通解是：
  $$h_n=H_n^{(1)}+H_n^{(2)}+\cdots +H_n^{(t)}$$
  但在实际应用上，我们很难求得多项式的所有根而受到限制。
  但是我们还是可以利用生成函数来求解（至少在理论上）任意 $k$ 阶常系数线性齐次递推关系。
  相应的生成函数是形如 $p(x)/q(x)$ 的函数，其中 $p(x)$ 是次数小于 $k$ 的多项式，而 $q(x)$ 是常数项等于 $1$ 的 $k$ 阶多项式。为了求数列项的生成函数，首先，我们用部分分式法把 $p(x)/q(x)$ 表示成如下形式的代数分式的和
  $$\frac{c}{(1-rx{)}^t}$$
  其中 $t$ 是正整数，$r$ 是实数，$c$ 是常数。于是我们可以求出该分式的幂级数，把项加起来，就得到生成函数的幂级数。
• 例子：设 $h_0,h_1,\cdots ,h_n,\cdots$ 是满足下面递推关系的数列：
  $$h_n+h_{n-1}-16h_{n-2}+20h_{n-3}=0(n\ge 3)$$
  其中 $h_0=0,h_1=1,h_2=-1$ 。求 $h_n$ 的通项公式。
  设 $g(x)={\sum }_{i=0}^{\infty }{h}_i{x}^i$ 是上述数列的生成函数。把下面四个方程
  $$\begin{array}{rl}g(x)& =h_0+h_{1}x+h_2{x}^2+\cdots +h_n{x}^n+\cdots \\ xg(x)& =h_{0}x+h_1{x}^2+\cdots +h_{n-1}{x}^n+\cdots \\ -16x^{2}g(x)& =-16h_0{x}^2-\cdots -16h_{n-2}{x}^n-\cdots \\ 20x^{3}g(x)& =20h_0{x}^3+\cdots +20h_{n-3}{x}^n+\cdots \end{array}$$
  加起来，我们得到
  $$(1+x-16x^2+20x^3)g(x)=h_0+(h_1+h_0)x+(h_2+h_1-16h_0)x^2$$
  带入 $h_0,h_1,h_2$，我们得到
  $$g(x)=\frac{x}{1+x-16x^2+20x^3}$$
  我们看到分子 $=(1-2x{)}^2(1+5x)$。因此我们能得到：
  $$g(x)=\frac{x}{1+x-16x^2+20x^3}=\frac{c_1}{1-2x}+\frac{c_2}{(1-2x{)}^2}+\frac{c_3}{1+5x}=\frac{-2/49}{1-2x}+\frac{7/49}{(1-2x{)}^2}+\frac{-5/49}{1+5x}$$
  然后，展开幂级数项，我们直接能得到：
  $$g(x)=\sum _{k=0}^{\infty }(-\frac{2}{49}{2}^k+\frac{7}{49}(k+1)2^k-\frac{5}{49}(-5{)}^k)x^k$$
  所以我们就得到了：
  $$h_n=-\frac{2}{49}{2}^n+\frac{7}{49}(n+1)2^n-\frac{5}{49}(-5{)}^n$$

非齐次递推关系

• 通常情况下，非齐次递推关系更难求解，而且通常需要依赖于其非齐次部分。
• 例子：(汉诺塔问题)
  设 $h_n$ 是转移 $n$ 个圆盘所需的移动次数。可以验证 $h_0=0,h_1=1,h_2=3$
  为了把 $n$ 个圆盘转移到另一根柱子上，我们必须首先把一根珠子上顶部 $n-1$ 个盘子转移到穿有最大圆盘的那根珠子上。因此
  $$h_n=2h_{n-1}+1(n\ge 1)$$
  通过反复递推带入，我们可以得到 $h_n=2^n-1$
• 例子：(汉诺塔的生成函数解法)
  首先我们有
  $$h_n=2h_{n-1}+1(n\ge 1),h_0=0$$
  设
  $$g(x)=\sum _{n=0}^{\infty }{h}_n{x}^n$$
  于是我们有：
  $$\begin{gathered} g(x)=h_0+h_{1}x+h_2{x}^2+\cdots +h_n{x}^n+\cdots \\ -2xg(x)=2h_{0}x+\cdots +2h_{n-1}{x}^n+\cdots \end{gathered}$$
  把上面两个方程相加，得到
  $$(1-2x)g(x)=x+x^2+\cdots +x^n+\cdots =\frac{1}{1-x}$$
  因此
  $$g(x)=\frac{x}{(1-x)(1-2x)}=\frac{1}{1-2x}-\frac{1}{1-x}=\sum _{n=0}^{\infty }(2^n-1)x^n$$

---

• 例子：求解：
  $$h_n=3h_{n-1}-4n(n\ge 4),h_0=2$$
  首先我们考虑相应的齐次递推关系 $h_n=3h_{n-1}$，它的通解是 $h_n=c{3}^n(n\ge 1)$
  下面我们要求原来的非齐次递推关系的一个特殊的解。对于适当的 $r,s$，我们尝试着寻找下面形式的一个解
  $$h_n=rn+s$$
  首先我们必须有
  $$\begin{gathered} rn+s=3(r(n-1)+s)-4n \\ rn+s=(3r-4)n+(-3r+3s) \end{gathered}$$
  我们得到
  $$\begin{gathered} r=3r-4 \\ s=-3r+3s \end{gathered}$$
  因此得到 $r=2,s=3$，且 $h_n=2n+3$ 满足。现在将齐次关系的通解与费其次关系的特殊解合成，得到
  $$h_n=c{3}^n+2n+3$$
  根据初始条件 $h_0=2$ 得到 $c=-1$，因此
  $$h_n=-3^n+2n+3$$
  是原来问题的解。

---

• 上述解法是求解费其次微分方程方法的离散模式。总结如下：
  求齐次关系的通解
  求非齐次关系的一个特殊解
  将通解和特殊解合成，确定通解中的常数值

---

• 如果 $b_n$ 是 $n$ 的 $k$ 次多项式，可以做如下尝试
  $h_n=r$ 如果 $b_n=d$
  $h_n=rn+s$ 如果 $b_n=dn+e$
  $h_n=r{n}^2+sn+t$ 如果 $b_n=d{n}^2+en+f$
  $h_n=p{d}^n$ 如果 $b_n=d^n$

---

• 例子：求解
  $$h_n=2h_{n-1}+3^n(n\ge 1),h_0=2$$
  解一：因为齐次关系 $h_n=2h_{n-1}$，通解为 $h_n=c{2}^n(n\ge 1)$
  为找其特殊解，我们做下面的尝试：$h_n=p{3}^n$
  那么必须满足 $p{3}^n=2p{3}^{n-1}+3^n$
  解得 $p=3$
  因此 $h_n=c{2}^n+3^{n+1}$
  带入 $h_0=2$，解得 $c=-1$
  于是得到 $h_n=-2^n+3^{n+1}(n\ge 0)$
  解二：利用生成函数。设 $g(x)={\sum }_{i=0}^{\infty }{h}_i{x}^i$
  利用原来的递推关系和 $h_0=2$，我们得到
  $$\begin{array}{rl}g(x)-2xg(x)& =h_0+(h_1-2h_0)x+\cdots +(h_n-2h_{n-1})x^n+\cdots \\ & =2+3x+3^2{x}^2+\cdots +3^n{x}^n+\cdots \\ & =2-1+(1+3x+3x^2{x}^2+\cdots +3^n{x}^n+\cdots )\\ & =1+\frac{1}{1-3x}\end{array}$$
  因此
  $$g(x)=\frac{1}{1-2x}+\frac{1}{(1-3x)(1-2x)}=\frac{1}{1-2x}+\frac{3}{1-3x}-\frac{2}{1-2x}=\sum _{n=0}^{\infty }(3^{n+1}-2^n)x^n$$

---

求解卡特兰数的例子

• 令 $h_n$ 表示把凸 $n+1$ 多边形区域分成三角形区域的方案数。则 $h_n$ 满足递推关系
  $$h_n=\sum _{k=1}^{n=1}{h}_k{h}_{n-k}$$
  我们现在希望求出它的通解。令 $h_1=1$，设
  $$\begin{array}{rl}g(x)& =\sum _{i=1}^{\infty }{h}_i{x}^i\\ (g(x){)}^2& =h_1^2{x}^2+\cdots +(h_1{h}_{n-1}+h_2{h}_{n-2}+\cdots +h_{n-1}{h}_1)x^n+\cdots \\ (g(x){)}^2& =h_1^2{x}^2+h_3{x}^3+h_4{x}^4+\cdots +h_n{x}^n+\cdots \\ & =g(x)-h_{1}x=g(x)-x\end{array}$$
  因此，$g(x)$ 满足方程
  $$(g(x){)}^2-g(x)+x=0$$
  根据二次公式，$g(x)$ 有两个根 $g_1(x),g_2(x)$，其中
  $$g_1(x)=\frac{1+\sqrt{1-4x}}{2},g_2(x)=\frac{1-\sqrt{1-4x}}{2}$$
  由 $g(x)$ 的定义可知 $g(0)=0$，因为 $g_1(0)=1,g_2(0)=0$，所以得到 $g(x)=g_2(x)$
  根据牛顿二项式定理得到
  $$(1+z{)}^{1/2}=1+\sum _{n=1}^{\infty }\frac{(-1{)}^{n-1}}{n\times 2^{2n-1}}{C}_{2n-2}^{n-1}{z}^n(|z|<1)$$
  用 $-4x$ 替代 $z$ ，则得到
  $$\begin{array}{rl}(1-4x{)}^{1/2}& =1+\sum _{n=1}^{\infty }\frac{(-1{)}^{n-1}}{n\times 2^{2n-1}}{C}_{2n-2}^{n-1}(-1{)}^n{4}^n{x}^n\\ & =1+\sum _{n=1}^{\infty }(-1{)}^{2n-1}\frac{2}{n}{C}_{2n-2}^{n-1}{x}^n\end{array}$$
  因此
  $$g(x)=\sum _{n=1}^{\infty }\frac{1}{n}{C}_{2n-2}^{n-1}{x}^n$$
  所以
  $$h_n=\frac{1}{n}{C}_{2n-2}^{n-1}$$