斐波那契矩阵快速

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<string>
#include<set>
#include<map>
#include<queue>
#include<vector>
using namespace std;
#define LL long long
#define maxn 100005
const int N=2;
const LL mod=1000000007;
struct node
{
    LL v[N][N];
};
node matrixmul(node A,node B)
{
    node ans;
    for(int i=0; i<N; i++)
        for(int j=0; j<N; j++)
        {
            ans.v[i][j]=0;
            for(int k=0; k<N; k++)
                ans.v[i][j]+=(A.v[i][k]*B.v[k][j])%mod;
            ans.v[i][j]=ans.v[i][j]%mod;
        }
    return ans;
}
node power(node a,LL n)
{
    node ans= {1,1,1,0};
    while(n)
    {
        if(n&1)
            ans=matrixmul(ans,a);
        a=matrixmul(a,a);
        n>>=1;
    }
    return ans;
}
int main()
{
    int t;
    LL n;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld",&n);
        n++;
        node C={1,1,1,0};
        if(n<2)puts("1");
        else {
        node ans=power(C,n-2);
        printf("%d\n",ans.v[0][0]);
       }
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值