洛谷P1962 斐波那契数列 （矩阵快速幂）

学了矩阵，练一下手。。。

 1 #include<bits/stdc++.h>
 2 typedef long long ll;
 3 const ll mod=1e9+7;
 4 using namespace std;
 5
 6 struct Matrix{
 7     ll g[3][3];
 8     Matrix() {
 9         memset(g,0,sizeof(g));
10     }//矩阵初始化为0
11     Matrix operator *(const Matrix &b) const{//重载乘号
12         Matrix res;
13         for(int i=1;i<=2;i++)
14             for(int j=1;j<=2;j++)
15                 for(int k=1;k<=2;k++)
16                     res.g[i][k]=(res.g[i][k]+g[i][j]*b.g[j][k]%mod)%mod;
17         return res;
18     }
19 }a,ans;
20
21 void init(){
22     a.g[1][1]=1,a.g[1][2]=1,a.g[2][1]=1;
23     ans.g[1][1]=ans.g[1][2]=1;
24 }
25
26 void qpow(ll x){//矩阵快速幂
27     while(x){
28         if(x&1) ans=ans*a;
29         a=a*a;
30         x>>=1;
31     }
32 }
33
34 int main(){
35     ll n;
36     cin>>n;
37     if(n<=2){
38         cout<<1;
39         return 0;
40     }
41     init();
42     qpow(n-2);
43     cout<<ans.g[1][1]%mod;
44     return 0;
45 }
46
47 //没开long long始终都过不了的数据
48 //65748392011234567 in
49 //188363182 out