矩阵快速幂&矩阵乘法的一些神奇操作

最新推荐文章于 2025-03-17 17:58:34 发布

置顶 Zhang_sir00

最新推荐文章于 2025-03-17 17:58:34 发布

阅读量447

点赞数 1

本文链接：https://blog.youkuaiyun.com/Zhang_sir00/article/details/104193435

版权

本文介绍了一种使用矩阵快速幂算法高效求解斐波那契数列第n项的方法，适用于n大于1e10的大数计算场景。通过定义特定的矩阵运算和快速幂操作，实现了对斐波那契数列的高效求解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

牛客第一场J题让我明白了自己有多少东西都还不会,今天来学习一下神(jian)奇(dan)的矩阵快速幂并如何利用矩阵快速幂求斐波那契数列第n项(n>1e10)

矩阵快速幂(斐波那契数列)

代码:

#include <iostream>
#include <cstring>
#include <string>
#define ll long long
#define mod 1000000007
using namespace std;
struct node {
    ll a[2][2];
    void init_1()                    //1 1
    {                                //0 0
        a[0][0] = a[0][1] = 1;
        a[1][0] = a[1][1] = 0;
    }

    void init_2()                     //1 1
    {                                 //1 0
        a[0][0] = a[0][1] = a[1][0] = 1;
        a[1][1] = 0;
    }

    void init_3()                     //1 0
    {                                 //0 1
        a[0][0] = a[1][1] = 1;
        a[0][1] = a[1][0] = 0;
    }

    friend ll getfabo(node a,node b,ll n)
    {
        if (n == 1)
            return 1;
        a.init_1();
        b.init_2();
        b = qmul(b,n-2);
        a = mul(a, b);
        return a.a[0][0];
    }

    friend node mul(node a,node b)//朴实无华的矩阵乘法
    {
        node ans;
        for (int i = 0; i < 2; i++)
        {
            for (int j = 0; j < 2; j++)
            {
                ans.a[i][j] = 0;
                for (int k = 0; k < 2; k++)
                {
                    ans.a[i][j] = (ans.a[i][j] + (a.a[i][k] * b.a[k][j] % mod)) % mod;
                }
            }
        }
        return ans;
    }

    friend node qmul(node a, ll n)//(神奇而简单的矩阵快速幂)和快速幂原理一样
    {
        node b;
        b.init_3();
        while (n)
        {
            if (n & 1)
                b = mul(a, b);
            n >>= 1;
            a = mul(a, a);
        }
        return b;
    }

}fabo1,fabo2;

int main()
{
    cout << getfabo(fabo1, fabo2, 1000000000000)<<endl;
    return 0;
}