矩阵二分幂求Fib（非结构体版）

来源：互联网发布：移动短信平台软件编辑：程序博客网时间：2024/05/22 02:27

矩阵二分幂求fib模板（0<=n<=10^9)，如果取模的数的平方没超过long long,那么n可以满足条件为n<=10^18,速度很快。

#include<stdio.h>#include<string.h>#include<algorithm>#define ll long long#define mod 1000000007using namespace std;ll solve(ll n){ll Ans[2][2]={1,0,0,1};ll Base[2][2]={0,1,1,1};ll temp00=0,temp01=0,temp10=0,temp11=0;while(n){if(n&1){temp00=Ans[0][0],temp01=Ans[0][1];temp10=Ans[1][0],temp11=Ans[1][1];Ans[0][0]=(temp00*Base[0][0]+temp01*Base[1][0])%mod;Ans[0][1]=(temp00*Base[0][1]+temp01*Base[1][1])%mod;Ans[1][0]=(temp10*Base[0][0]+temp11*Base[1][0])%mod;Ans[1][1]=(temp10*Base[0][1]+temp11*Base[1][1])%mod;}temp00=Base[0][0],temp01=Base[0][1];temp10=Base[1][0],temp11=Base[1][1];Base[0][0]=(temp00*temp00+temp01*temp10)%mod;Base[0][1]=(temp00*temp01+temp01*temp11)%mod;Base[1][0]=(temp10*temp00+temp11*temp10)%mod;Base[1][1]=(temp10*temp01+temp11*temp11)%mod;n>>=1;}return Ans[0][1];}int main(){    int t;    scanf("%d",&t);    while(t--){    ll n;    scanf("%lld",&n);        ll ans = solve(n);        printf("%lld\n",ans);}return 0;}

0 0