问题 B: 【例题2】Fibonacci第n项

本文介绍了一种使用矩阵快速幂算法高效求解大规模Fibonacci数列问题的方法,通过矩阵递推公式,避免了传统递归或迭代算法的时间复杂度过高的问题。文章提供了详细的代码实现,包括矩阵乘法和矩阵快速幂的具体步骤。

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

问题 B: 【例题2Fibonaccin

时间限制: 1 Sec  内存限制: 128 MB
提交: 18  解决: 10
[提交][状态][讨论版][
命题人:quanxing][Edit] [TestData]

题目链接:http://acm.ocrosoft.com/problem.php?cid=1704&pid=1

题目描述

大家都知道Fibonacci数列吧,f[1]=1,f[2]=1,f[3]=2,f[4]=3.......也就是f[n]=f[n-1]+f[n-2]。现在,问题很简单,输入nm,求第n项取模m

输入

输入nm

1<=n<=2 000 000 000

1<=m<=1 000 000 010

输出

输出第n项取模m

样例输入

5 1000

样例输出

5

思路:矩阵快速幂,加推导斐波那契数列的矩阵递推式

代码:

#include<bits/stdc++.h>

using namespace std;

#define ll long long

struct MUL

{

    ll m[3][3];

}ans, res;

MUL mul(MUL a, MUL b, int n, int mod)//矩阵乘法,定义结构体函数,传入结构体,返回结构体,这样比较方便

{

    MUL tmp;

    for (int i = 1; i <= n; i++)

    {

         for (int j = 1; j <= n; j++)

         {

             tmp.m[i][j] = 0;

         }

    }

    for (int i = 1; i <= n; i++)

    {

         for (int j = 1; j <= n; j++)

         {

             for (int k = 1; k <= n; k++)

             {

                  tmp.m[i][j] += ((a.m[i][k] % mod) *(b.m[k][j] % mod)) % mod;

             }

         }

    }

    return tmp;

}

ll quickpower(int N, int n, int mod)//矩阵快速幂

{

    MUL C, res;//C是常数矩阵

    memset(res.m, 0, sizeof(res.m));

    C.m[1][1] = 1;

    C.m[1][2] = 1;

    C.m[2][1] = 1;

    C.m[2][2] = 0;

    for (int i = 1; i <= n; i++)res.m[i][i] = 1;//单位矩阵

    while (N)

    {

         if (N & 1)

             res = mul(res, C, n, mod);

         C = mul(C, C, n, mod);

         N = N >> 1;

    }

    return res.m[1][1];

}

int main()

{

    ll n, mod;

    cin >> n >> mod;

    cout << quickpower(n - 1, 2, mod) % mod << endl;

}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值