……续上回
fss.sosei:斐波那契数列与Python的尾递归蹦床 连载【7】zhuanlan.zhihu.com之前的篇章把各种Fibonacci数列的基本算法讨论过了
那么是否可以做到更快呢,有什么加速手段
这篇来说下
首先第一个手段是改进算法的加速
16. 快速平方的矩阵解法
在Fibonacci数列高效解法大全及时间复杂度分析 连载【5】这篇里说过矩阵解法
矩阵法虽然跟二进制模幂解法时间复杂度一样,可算第100万项斐波那契数用时是二进制模幂解法的10倍。这是因为这算法的时间常数项大
里面用到了矩阵乘法,通用矩阵乘法算法的时间复杂度是阶数n的O(n^3)。也就是对一个二阶矩阵,分解步骤中有8次乘法,非常耗时,造成矩阵解法时间常数项很大
之前程序是直接用的numpy库做矩阵乘法,numpy库里就是通用矩阵乘法算法
然而,对斐波那契数的矩阵解法,只需对二阶矩阵做运算。而其中的二阶矩阵平方步骤是有更小时间复杂度算法的
二阶矩阵快速平方算法的时间复杂度为阶数2的O(2^(log2(5)))
也就是二阶矩阵一次平方运算分解步骤中只用做5次乘法。比通用矩阵算法的8次降低了不少
但是,对斐波那契数矩阵解法而言,矩阵乘法步骤是其中的时间常数项,缩短矩阵乘法用时只是降低整个算法的常数项,可并不降低整个算法的时间复杂度
下面就是我写的程序
import
import numpy.matlib
from gmpy2 import mpz
def Fibonacci_sequence_06 (n: int) -> int: #参数n是表示求第n项Fibonacci数
'返回单项的二阶矩阵快速平方解法'
assert isinstance(n, int), 'n is an error of non-integer type.'
Use_threshold_of_fast_power = 100000
if n >= 2:
fib_matrix = numpy.mat(((mpz(1), mpz(1)), (mpz(1), mpz(0))))
if n > Use_threshold_of_fast_power:
fib_matrix = fast_power_of_second_order_matrix(fib_matrix, n - 1)
else:
fib_matrix **= n - 1
return fib_matrix[0,0]
elif n == 1:
return 1
elif n == 0:
return 0
else:
return None
Fibonacci_sequence_06(1000000)
用算到第100万项Fibonacci数来测量下用时
Total time: 0.036093秒
比之前用通用矩阵乘法库的程序缩短了15%的时间
注意,那么一大段程序是有基础开销的
所以在非大数的时候用这个反而会不如通用库快
怎么办呢
加一个阈值判别。当N大于阈值之后,才用这个二阶矩阵快速平法算法,数不大时就调用通用库
搞定
------------------------
算Fibonacci数的快速算法已经写了4种,算法时间复杂度皆为O(log n),但常数项不同。
矩阵快速幂解法每步中两个大整数乘法是8次
快速平方的矩阵解法每步中两个大整数乘法是5次
二分解法每步中两个大整数乘法是3次
二进制模幂解法每步中两个大整数乘法是2次
有趣,刚好是个Fibonacci数列
欢迎点赞、留言
未完待续……
fss.sosei:斐波那契数列与Python的尾递归蹦床 连载【9】zhuanlan.zhihu.com