Fibonacci数列

Zewbie

于 2021-08-24 12:06:24 发布

阅读量188

点赞数

分类专栏：笔记文章标签：数学算法矩阵

本文链接：https://blog.csdn.net/qq_44729222/article/details/119886947

版权

笔记专栏收录该内容

6 篇文章 1 订阅

订阅专栏

$F i b o n a c c i$ 数列： $F_0=0,F_1=1,F_n=F_{n-1}+F_{n-2}(n\ge 2)$ .

本文重点关注 $F i b o n a c c i$ 数列的通项公式与第 $n$ 项的计算.

文章目录

一、通项公式
二、计算

一、通项公式

构造矩阵递推式：
$\begin{cases} F_{n+1}=F_n+F_{n-1},\\ F_{n}=F_n+0\cdot F_{n-1}. \end{cases}$

$\begin{aligned} \Rightarrow \begin{pmatrix} F_{n+1} \\ F_{n} \end{pmatrix} &=\begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}\begin{pmatrix} F_{n} \\ F_{n-1} \end{pmatrix}\\ &=\begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}^{2}\begin{pmatrix} F_{n-1} \\ F_{n-2} \end{pmatrix}\\ &=\cdots \\ &=\begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}^{n}\begin{pmatrix} F_{1} \\ F_{0} \end{pmatrix}. \end{aligned}$

令 $A=\begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}$ ， $A$ 为实对称矩阵，因此必定存在可逆矩阵 $P$ ，使 $P^{-1}AP=\Lambda$ .

取 $P=\begin{pmatrix} -2 & 2 \\ \sqrt{5}+1 & \sqrt{5}-1 \end{pmatrix}$ ，则

$P^{-1}=\begin{pmatrix} \frac{1}{10+2\sqrt{5}} & \\ & \frac{1}{10-2\sqrt{5}} \end{pmatrix}\begin{pmatrix} -2 & \sqrt{5}+1 \\ 2 & \sqrt{5}-1 \end{pmatrix}$ ，

满足 $P^{-1}AP=\Lambda=\begin{pmatrix} \dfrac{1-\sqrt{5}}{2} & \\ & \dfrac{1+\sqrt{5}}{2} \end{pmatrix}$ .

于是
$\begin{aligned} \Rightarrow \begin{pmatrix} F_{n+1} \\ F_{n} \end{pmatrix} &=A^{n}\begin{pmatrix} F_{1} \\ F_{0} \end{pmatrix}\\ &=P\Lambda ^{n}P^{-1}\begin{pmatrix} 1 \\ 0 \end{pmatrix}\\ &=\begin{pmatrix} \dfrac{-1}{\sqrt{5}}\left( \dfrac{1-\sqrt{5}}{2}\right) ^{n+1} +\dfrac{1}{\sqrt{5}}\left( \dfrac{1+\sqrt{5}}{2}\right) ^{n+1} \\ \dfrac{-1}{\sqrt{5}}\left( \dfrac{1-\sqrt{5}}{2}\right) ^{n} +\dfrac{1}{\sqrt{5}}\left( \dfrac{1+\sqrt{5}}{2}\right) ^{n} \end{pmatrix}. \end{aligned}$
故 $F i b o n a c c i$ 数列的通项公式为
$F_n=\dfrac{1}{\sqrt{5}}\left( \dfrac{1+\sqrt{5}}{2}\right) ^{n}-\dfrac{1}{\sqrt{5}}\left( \dfrac{1-\sqrt{5}}{2}\right) ^{n}.(n\ge 0)$
由于
$\left| F_{n}-\dfrac{1}{\sqrt{5}}\left( \dfrac{1+\sqrt{5}}{2}\right) ^{n}\right| =\dfrac{1}{\sqrt{5}}\left( \dfrac{\sqrt{5}-1}{2}\right) ^{n}\le \dfrac{1}{\sqrt{5}} \lt\dfrac 12.(n\ge 0)$
因此
$F_{n}=\left\lfloor \dfrac{1}{\sqrt{5}}\left( \dfrac{\sqrt{5}+1}{2}\right) ^{n}+\dfrac{1}{2}\right\rfloor.(n\ge 0)$

二、计算

方法一：利用递推公式

$\begin{cases}F_{n}=F_{n-1}+F_{n-2},n\geq 2,\\ F_{0}=0,F_{1}=1.\end{cases}$

时间复杂度 $O (n)$ .

def get_fib1(n: int):
    """递推式求Fibonacci第n项"""
    if n < 0:
        return None
    if n == 0 or n == 1:
        return n
    a, b = 0, 1
    for i in range(2, n + 1):
        a += b
        a, b = b, a
    return b

方法二：矩阵快速幂

令 $A=\begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}$ ，则 $\begin{pmatrix} F_{n+1} \\ F_{n} \end{pmatrix}=A^{n}\begin{pmatrix} 1 \\ 0 \end{pmatrix}$ .

显然 $F_n=A^n_{2,1}$ .

时间复杂度 $O(\log_2n)$ .

$F_{47}$ 超过 $4$ 字节整数范围导致结果溢出.

def get_fib2(n: int):
    """矩阵快速幂计算F_n
        n=47时溢出
    """
    if n < 0:
        return None
    import numpy as np
    tmp = [[1, 0], [0, 1]]
    a = [[1, 1], [1, 0]]
    while n > 0:
        if n & 1:
            tmp = np.dot(tmp, a)
        a = np.dot(a, a)
        n >>= 1
    return tmp[1][0]

方法三：通项公式计算

通项公式：
$F_{n}=\left\lfloor \dfrac{1}{\sqrt{5}}\left( \dfrac{\sqrt{5}+1}{2}\right) ^{n}+\dfrac{1}{2}\right\rfloor.$
时间复杂度取决于计算 ((1 + sqrt5) / 2) ** n 的速度.

$n = 71$ 时会由于浮点误差过大导致结果不对.

def get_fib3(n: int):
    """利用通项公式求F_n
        n=71时浮点误差过大
    """
    if n < 0:
        return None
    import math
    sqrt5 = math.sqrt(5)
    return int(1 / sqrt5 * ((1 + sqrt5) / 2) ** n + 0.5)

使用高精度可以在一定程度上解决浮点误差问题.

若运算精度取 $28$ ， $n = 123$ 时结果错误.

若运算精度取 $100$ ， $n = 464$ 时结果错误.

若运算精度取 $1000$ ， $n = 4768$ 时结果错误.

def get_fib4(n: int):
    """利用通项公式求F_n"""
    import decimal
    from decimal import Decimal
    decimal.getcontext().prec = 1000  # 运算精度，默认值 28
    sqrt5 = Decimal(5).sqrt()
    _1 = Decimal(1)
    _2 = Decimal(2)
    return int(_1 / sqrt5 * ((_1 + sqrt5) / _2) ** n + _1 / _2)