矩阵快速幂求斐波那契数列

最新推荐文章于 2024-05-31 11:57:04 发布

尾流

最新推荐文章于 2024-05-31 11:57:04 发布

阅读量94

点赞数

文章标签：矩阵线性代数

本文链接：https://blog.csdn.net/weixin_62868784/article/details/128710079

版权

原视频

【竞赛向】斐波那契数列的矩阵快速幂求法_哔哩哔哩_bilibili最后用python写的，没用numpy之类的，因为算法竞赛时不能用第三方库，这个视频也许对蓝桥杯有帮助。预备知识：线性代数。绘图软件：ipad上的GoodNotes, 视频播放量 4943、弹幕量 4、点赞数 283、投硬币枚数 75、收藏人数 330、转发人数 40, 视频作者阿-岳同学, 作者简介学习区和编程区up主、时而认真，时而不正经，相关视频：TimSort 一个几乎没人知道的排序算法 | 时间复杂度最快达到了o(n)，小兔老师数正方体原视频(非恶意发送)，[证毕QED]如何用3个π表示自然数，二分查找为什么总是写错？，质数通项公式，《 n 倍角公式》，KMP匹配太麻烦了？试试巨简单的RK算法，切比雪夫定理，告诉我们一些任何人永远积不出来的函数！，上海市小升初考试题，A，B，C，D为正方形四边中点，求问号处面积，1.01^100与2谁大谁小？https://www.bilibili.com/video/BV1GA411Q79k/?spm_id_from=333.1007.top_right_bar_window_history.content.click&vd_source=d3ea9fa6b291ac6cd33aae43110b598e

def matmul(m1,m2):
    res = [[0] * len(m2[0]) for _ in range(len(m1))]
    for y in range(len(m1)):
        for x in range(len(m2[0])):
            n = 0
            for a, b in zip(m1[y], (line[x] for line in m2)):
                n += a * b
            res[y][x] = n
    return res
def fab(n):
    if n < 0:
        return 0
    m = [
        [0,1],
        [1,1],
    ]
    resM = [[1,0],[0,1]]
    while n != 0:
        if n & 1 == 1:
            resM = matmul(m ,resM)
        m = matmul(m, m)
        n >>= 1
    c, d = resM[1]
    return d