Computational Complexity of Fibonacci Sequence / 斐波那契数列的时空复杂度

Fibonacci Sequence 维基百科

\(F(n) = F(n-1)+F(n-2)\),其中 \(F(0)=0, F(1)=1\),即该数列由 0 和 1 开始,之后的数字由相邻的前两项相加而得出。

递归

def fibonacci(n):
    assert n >= 0, 'invalid n'
    if n < 2: return n
    return fibonacci(n - 1) + fibonacci(n -2)

递归方法的时间复杂度为高度为 \(n-1\) 的不完全二叉树的节点数,所以近似为 \(O(2^n)\)

数学求解方法:

\(T(n) = T(n-1) + T(n-2) \qquad (n>1)\)

设 \(f(n)\) 为参数为 n 时的时间复杂度 \(f(n) = f(n-1) + f(n-2)\)

转化为求二阶常系数齐次差分方程,设通解为 \(f_{n} = C_{1}f_{n-1} + C_{2}f_{n-2}\)

设有特解 \(f_{n} = \lambda\),\(\lambda\) 为非零待定常数,将 \(\lambda\) 代入方程,易得特征方程 \(\lambda ^{2} = \lambda + 1\),则 \(λ = \frac{1\pm \sqrt{5}}{2}\)

再根据 \(f(0) = 0, f(1) = 1.\) 求出 \(C_{1}\) 和 \(C_{2}\)

得出通项公式

\[
f(n) = \frac{1}{\sqrt{5}}[(\frac{1+ \sqrt{5}}{2})^{n} - (\frac{1- \sqrt{5}}{2})^{n}]
\]

当 \(n->\infty\) 时, \(\left |(\frac{1- \sqrt{5}}{2})^{n} \right |<1\),所以趋于零。

时间复杂度为 \(O((\frac{1+ \sqrt{5}}{2})^{n})\) ,约等于 \(O(1.618^{n})\) ,即指数级复杂度 \(O(2^n)\),递归算法空间复杂度取决于递归的深度,显然为 \(O(n)\)。

迭代

def fibonacci(n):
    assert n >= 0, 'invalid n'
    if n == 0: return 0
    a, b = 0, 1
    for _ in range(n - 1):
        a, b = b, a+b
    return b

时间复杂度 \(O(n)\),空间复杂度 \(O(1)\)

矩阵

\(F(n)\) 和 \(F(n - 1)\) 写成一个 2 x 1 的矩阵,然后对其进行变形。

\(\begin{bmatrix}F_{n}\\F_{n-1}\end{bmatrix}=\begin{bmatrix}F_{n-1}+F_{n-2}\\F_{n-1}\end{bmatrix}=\begin{bmatrix}1\times F_{n-1}+1\times F_{n-2}\\1\times F_{n-1}+0\times F_{n-2}\end{bmatrix}=\begin{bmatrix}1 & 1\\ 1 & 0\end{bmatrix}\times\begin{bmatrix}F_{n-1}\\F_{n-2}\end{bmatrix}=...=\begin{bmatrix}1&1\\1&0\end{bmatrix}^{n-1}\times\begin{bmatrix}F_{1}\\F_{0}\end{bmatrix}=\begin{bmatrix}1&1\\1&0\end{bmatrix}^{n-1}\times\begin{bmatrix}1\\0\end{bmatrix}\)

因此要求 \(F_{n}\),只要对这个二阶方阵求 \(n - 1\) 次方,最后取结果方阵第一行第一列的数字就可以了。

等式中的矩阵 \(\begin{bmatrix}1&1\\1&0\end{bmatrix}\) 被称为斐波那契数列的 Q- 矩阵。

通过 Q- 矩阵,我们可以利用如下公式进行计算:

\(F_{n} = Q^{n-1}_{1,1}\)

如此一来,计算斐波那契数列的问题就转化为了求 \(Q^{n-1}\) 的问题。

\[
A^{n} = \begin{cases} A\cdot (A^{2})^{\frac{n-1}{2}}, & \text{if $n$ is odd} \\\ (A^{2})^{\frac{n}{2}}, & \text{if $n$ is even} \\ \end{cases}
\]

可见时间复杂度满足 \(T(n) = T(n / 2) + O(1)\)
Master 定理 可得

时间复杂度 \(O( \log n)\) ,空间复杂度显然为 \(O(1)\)

from numpy import matrix

def MatrixPower(mat, n):
    assert n > 0, 'invalid n'
    res = None
    temp = mat
    while True:
      if n & 1:
        if res is None: res = temp
        else: res = res * temp
      n >>= 1
      if n == 0: break
      temp = temp * temp
    return res

def fibonacci(n):
    assert n >= 0, 'invalid n'
    if n < 2: return n  # F(0) = 0, F(1) = 1
    mat = matrix([[1, 1], [1, 0]], dtype=object)
    mat = MatrixPower(mat, n - 1)
    return mat[0, 0]

通项公式法

def fibonacci(n):
    root_Five = 5**0.5
    result = (((1 + root_Five) / 2)**n - ((1 - root_Five) / 2)**n) / root_Five
    return int(result)

显然,该方法的时空复杂度均为 \(O(1)\) , 但使用公式计算的方法由于有大量的浮点运算,在 n 增大时浮点误差不断增大会导致返回结果不正确甚至数据溢出。

方法 时间复杂度 空间复杂度
递归 \(O(2^{n})\) \(O(n)\)
迭代 \(O(n)\) \(O(1)\)
矩阵 \(O(log^{n})\) \(O(1)\)
公式 \(O(1)\) \(O(1)\)

原文地址:https://www.cnblogs.com/yexuesong/p/12321448.html

时间: 2024-10-15 01:23:46

Computational Complexity of Fibonacci Sequence / 斐波那契数列的时空复杂度的相关文章

Fibonacci series(斐波纳契数列)的几种常见实现方式

费波那契数列的定义: 费波那契数列(意大利语:Successione di Fibonacci),又译费波拿契数.斐波那契数列.斐波那契数列.黄金分割数列. 在数学上,费波那契数列是以递归的方法来定义: (n≧2) 用文字来说,就是费波那契数列由0和1开始,之后的费波那契系数就由之前的两数相加. 首几个费波那契系数是:0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233-- 特别指出:0不是第一项,而是第零项. 下面是费波那契数列的几种常见编程实现:

POJ 3070 Fibonacci【斐波那契数列/矩阵快速幂】

Fibonacci Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 17171   Accepted: 11999 Description In the Fibonacci integer sequence, F0 = 0, F1 = 1, and Fn = Fn − 1 + Fn − 2 for n ≥ 2. For example, the first ten terms of the Fibonacci sequen

Fibonacci(斐波那契数列)的第N位数的第N位数

无穷数列1,1,2,3,5,8,13,21,34,55...称为Fibonacci数列,它可以递归地定义为F(n)=1 ...........(n=1或n=2)F(n)=F(n-1)+F(n-2).....(n>2)现要你来求第n个斐波纳奇数.(第1个.第二个都为1): 输入第一行是一个整数m(m<5)表示共有m组测试数据每次测试数据只有一行,且只有一个整形数n(输出对每组输入n,输出第n个Fibonacci数 #include<stdio.h> long int fib(long

python实现斐波那契数列(Fibonacci sequence)

使用Python实现斐波那契数列(Fibonacci sequence) 斐波那契数列形如 1,1,2,3,5,8,13,等等.也就是说,下一个值是序列中前两个值之和.写一个函数,给定N,返回第N个斐波那契数字.例如,1返回1 6返回8 我选择了两种方法,一种是将list变成一个队列,另一个则是使用环形队列.不多说,直接上代码:后面我会对为什么这样实现做一个解释 第一个是使用队列的方式: 1 def fibonacciSeq(num): 2 fibonacciSeqList = [] 3 for

递归:阶乘、斐波那契数列

阶乘 public static void main(String[] args) { System.out.println(factorial(5));; } //factorial 阶乘 public static long factorial(int n) { if (n == 1) { return 1; } return n*factorial(n-1); } 阶乘画图理解 斐波那契数列递归实现: public static void main(String[] args) { Sys

python实现斐波那契数列

斐波那契数列的发明者是意大利数学家昂纳多.斐波那契(Leonardo Fibonacci).斐波那契数列又被称为黄金分割数列,或兔子数列.它指的是这样一个数列:0 1 1 2 3 5 8 13 21 34 ....在数学上,斐波那契数列以递归的方法定义:F(0)=0,F(1)=1,F(n)=F(n-1)+F(n-2(N>=2,n属于N*).简单的说明斐波那切数列的规律为:第1个数为0,第2个数为1,之后每个数值都是前两位的和. #!/usr/bin/python3 # coding=utf-8

【笔试】27、斐波那契数列

写一个函数,输入n,求斐波那契数列的第n项.斐波那契数列的定义如下 我们的第一反应是使用递归来做,但是实际上递归的效率非常之低 public static long TiaoJi(int n) { if(n <= 0) return 0; if(n == 1) return 1; return TiaoJi(n - 1) + TiaoJi(n - 2); } 其实我们可以用循环 /** *题目:1.写一个函数,输入n,求斐波那契数列的第n项. * 2.一只青蛙一次可以跳上1级台阶,也可以跳上2级

斐波拉契数列带来的思考

无意间看到一道题: 编写一个计算前100位斐波那契数的函数.根据定义,斐波那契序列的前两位数字是0和1,随后的每个数字是前两个数字的和.例如,前10位斐波那契数为:0,1,1,2,3,5,8,13,21,34. 我一看这不是熟悉的斐波那契数列嘛,简单!一个递归搞定,就没重视了.然后晚上自己尝试写着玩,却发现原来自己智力被爆. 首先写了个这样的程序: public class FibonnacciA{ public static void main(String[] args){ System.o

SQL Server -&gt;&gt; 斐波那契数列(Fibonacci sequence)

斐波那契数列(Fibonacci sequence)的T-SQL实现 ;WITH T AS ( SELECT 1 AS NUM, CAST(1 AS BIGINT) AS curr, CAST(NULL AS BIGINT) AS prv UNION ALL SELECT curr.NUM + 1 AS NUM, CAST(CASE WHEN prv IS NULL THEN curr ELSE curr + prv END AS BIGINT) AS curr, CAST(curr AS BI