若干递推数列的通项公式求解方法

最新推荐文章于 2022-02-26 11:29:21 发布

hezlik

最新推荐文章于 2022-02-26 11:29:21 发布

阅读量1.1k

点赞数 3

分类专栏：算法入门

本文链接：https://blog.csdn.net/hzk_cpp/article/details/102537500

版权

算法入门专栏收录该内容

43 篇文章 7 订阅

订阅专栏

一.等差与等比数列.

等差数列：以下数列被称为等差数列：
$f_n=a+f_{n-1}$

这个数列的通项公式为：
$f_{n}=f_{0}+na$

设 $s_{n}=\sum_{i=0}^{n}f_{i}$ ，则有：
$s_{n}=(n+1)f_{0}+a\frac{n(n+1)}{2}$

这个公式是通过将 $f_{i}$ 和 $f_{n-i+1}$ 配起来得到的.

等比数列：以下数列被称为等比数列：
$f_{n}=af_{n-1}$

这个数列的通项公式为：
$f_{n}=a^{n}f_{0}$

设 $s_{n}=\sum_{i=0}^{n}f_{i}$ ，则有：
$s_{n}=f_0\frac{1-a^{n+1}}{1-a}$

这个公式可以用错位相减法得到，具体过程如下：
$s_{n}=\sum_{i=0}^{n}f_i\\ =\sum_{i=0}^{n}\frac{f_{i+1}-f_{i}}{a-1}\\ =\frac{1}{a-1}\sum_{i=0}^{n}(f_{i+1}-f_{i})\\ =\frac{1}{a-1}(f_{1}-f_{0}+f_{2}-f_{1}+\cdots+f_{n+1}-f_{n})\\ =\frac{f_{n+1}-f_{0}}{a-1}\\ =f_{0}\frac{1-a^{n+1}}{1-a}$

将数列转化为等比数列是求数列通项公式的一个重要方法.

二.等比数列求解数列通项式的应用.

考虑一类一阶的递推数列：
$f_{n}=af_{n-1}+b\,\,(a\neq 1)$

先假设存在一个 $t$ 满足：
$f_{n}=af_{n-1}+(a-1)t$

那么有：
$f_{n}+t=a(f_{n-1}+t)$

这时就可以让 $b = (a - 1) t$ ，有：
$f_{n+1}+\frac{b}{a-1}=a\left(f_{n}+\frac{b}{a-1}\right)$

此时我们知道了数列 $g_n=f_{n}+\frac{b}{a-1}$ 是一个等比数列了.

根据等比数列通项公式有：
$f_{n}+\frac{b}{a-1}=a^{n}\left(f_{0}+\frac{b}{a-1}\right)\\ f_{n}=a^nf_0+b\frac{1-a^{n}}{1-a}$

三.特征方程求解数列通项式.

其实特征方程是线性代数里的一个东西，不过这里我们先介绍求解二阶递推式的特征方程.

仿照上面的方式，我们来求一求一类二阶递推数列的通项公式：
$f_{n+1}=af_{n}+bf_{n-1}$

假设存在 $t$ 满足：
$f_{n+1}+tf_{n}=c(f_{n}+tf_{n-1})$

那么有：
$f_{n+1}=(c-t)f_{n}+ctf_{n-1}$

转化为：
$\left\{\begin{matrix} c-t=a\\ ct=b \end{matrix}\right.$

消去 $t$ ，则有：
$c^{2}-ac-b=0\\$

这就是所谓的特征方程.

若这个特征方程的两根不同，根据等比数列通项公式就有：
$f_{n+1}+t_1f_{n}=c_{1}^{n}(f_{1}+t_1f_{0})\\ f_{n+1}+t_2f_{n}=c_{2}^{n}(f_{1}+t_2f_{0})$

消去 $f_{n+1}$ 得：
$f_{n}=c_{1}^{n}\frac{f_1+t_{1}f_0}{t_1-t_2}-c_{2}^{n}\frac{f_{1}+t_2f_{0}}{t_{1}-t_{2}}$

否则特征方程的两个解相同，那么 $c=\frac{a}{2}$ ，此时有：
$t=-\frac{a}{2}=-c$

接下来我们开始操作：
$\frac{f_{n+1}}{c^{n+1}}-\frac{f_{n}}{c^{n}}=\frac{f_{n+1}-cf_{n}}{c^{n+1}}=\frac{f_{n+1}+tf_{n}}{c^{n+1}}=\frac{f_{n}+tf_{n-1}}{c^{n}}=\cdots=\frac{f_{1}+tf_{0}}{c}=\frac{f_{1}-cf_{0}}{c}$

故 $g_{n}=\frac{f_{n}}{c^{n}}$ 必然是一个等差数列.

那么根据等差数列通项公式有：
$f_{n}=c^{n}\left(f_{0}+n\frac{f_1+tf_0}{c}\right)$

于是接下来我们可以总结，对于一个递推数列满足 $f_{n}=af_{n-1}+bf_{n-2}$ ，若其特征方程 $x^{2}=ax+b$ 的解为 $c_1,c_2$ ，则有：
1. $c_1\neq c_2$ 时，存在 $t_1,t_2$ 满足 $f_{n}=c_1^{n}t_1-c_2^{n}t_2$ .
2. $c_1=c_2$ 时，存在 $t_1,t_2$ 满足 $f_{n}=c^{n}(t_1+nt_2)$ .

其中 $t_1,t_2$ 均可通过初项确定.

特征方程求解数列通项式的一个经典例子是Fibonacci数列，即 $f_0=0,f_1=1,f_n=f_{n-1}+f_{n-2}$ .

列出其特征方程：
$x^{2}-x-1=0$

解得：
$x=\frac{1\pm \sqrt{5}}{2}$

根据上面的结论，直接设 $t_1,t_2$ 使得：
$f_{n}=t_1\left(\frac{1+\sqrt{5}}{2}\right)^{n}-t_2\left(\frac{1-\sqrt{5}}{2}\right)^{n}$

用 $f_0$ 和 $f_1$ 列出方程得：
$\left\{\begin{matrix} 0=t_1-t_2\\ 1=\frac{1+\sqrt{5}}{2}t_1-\frac{1-\sqrt{5}}{2}t_2 \end{matrix}\right.$

解得 $t_1=t_2=\frac{\sqrt{5}}{5}$ ，Fibonacci数列的通项公式即为：
$f_{n}=\frac{1}{\sqrt{5}}\left(\frac{1+\sqrt{5}}{2}\right)^{n}-\frac{1}{\sqrt{5}}\left(\frac{1-\sqrt{5}}{2}\right)^{n}$

四.生成函数方法.

生成函数是一类非常通用的组合方法，具体参见生成函数相关.

生成函数求解数列递推式的常见思路是利用递推式列出生成函数关系式，然后解方程并通过各种方法展开.

例如经典的Catalan数列，设Catalan数列为 $C_{n}$ ，那么有：
$C_n=\sum_{i=0}^{n-1}C_{i}C_{n-1-i}$

初始 $C_{0}=1$ .

设其生成函数为 $C (x)$ ，那么我们有：
$C(x)=xC^2(x)+1\\ xC^2(x)-C(x)+1\\ C(x)=\frac{1\pm \sqrt{1-4x}}{2x}$

由于 $C (x)$ 的常数项为 $1$ ，所以有 $C (0) = 1$ ，利用极限理论对两根验证：
$\lim_{x\rightarrow 0}\frac{1+\sqrt{1-4x}}{2x}=\lim_{x\rightarrow 0}\frac{2}{2x}=+\infty\\ \lim_{x\rightarrow 0}\frac{1-\sqrt{1-4x}}{2x}=\lim_{x\rightarrow 0}\frac{1}{\sqrt{1-4x}}=1$

我们发现取负的时候才是正确的，于是考虑暴力展开：
$C(x)=\frac{1-\sqrt{1-4x}}{2x}\\ =\frac{1}{2x}-\frac{1}{2x}\sum_{i=0}^{+\infty}\binom{\frac{1}{2}}{i}(-1)^{i}(4x)^{i}\\ =\frac{1}{2x}\sum_{i=1}^{+\infty}(-1)^{i-1}2^{2i}\binom{\frac{1}{2}}{i}x^{i}$

考虑求解 $\binom{\frac{1}{2}}{n}$ ：
$\binom{\frac{1}{2}}{n}\\ =\frac{\left(\frac{1}{2}\right)^{\underline{n}}}{n!}\\ =\frac{1}{n!}\prod_{i=1}^{n}\left(\frac{3}{2}-i\right)\\ =\frac{(-1)^{n-1}}{n!}2^{-n}\prod_{i=1}^{n-1}(2i-1)\\ =\frac{(-1)^{n-1}}{n!}2^{-n}\frac{\prod_{i=1}^{2n-2}i}{\prod_{i=1}^{n-1}2i}\\ =(-1)^{n-1}2^{1-2n}\frac{(2n-2)!}{n!(n-1)!}$

代入得到：
$C(x)=\frac{1}{2x}\sum_{i=1}^{+\infty}(-1)^{i-1}2^{2i}(-1) ^{i-1}2^{1-2i}\frac{(2n-2)!}{i!(i-1)!}x^{i}\\ =\frac{1}{x}\sum_{i=1}^{+\infty}\frac{(2i-2)!}{i!(i-1)!}x^{i}\\ =\sum_{i=0}^{+\infty}\frac{(2i)!}{(i+1)!i!}x^{i}\\ =\sum_{i=0}^{+\infty}\frac{1}{i+1}\binom{2i}{i}x^{i}$