矩陣快速冪及斐波那契數列模板

本篇博客先給出矩陣快速冪以及利用矩陣快速冪求斐波那契數列的模板，講解待更新……?

const int N=10;
int tmp[N][N];
void multi(int a[][N],int b[][N],int n)
{memset(tmp,0,sizeof tmp);for(int i=0;i<n;i++)for(int j=0;j<n;j++)for(int k=0;k<n;k++)tmp[i][j]+=a[i][k]*b[k][j];for(int i=0;i<n;i++)for(int j=0;j<n;j++)a[i][j]=tmp[i][j];
}
int res[N][N];
void Pow(int a[][N],int n)
{memset(res,0,sizeof res);//n是冪，N是矩陣大小for(int i=0;i<N;i++) res[i][i]=1;while(n){if(n&1)multi(res,a,N);//res=res*a;復制直接在multi里面實現了；multi(a,a,N);//a=a*an>>=1;}
}

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2;
const int MOD=1000000009;
struct mat
{ll a[N][N];
};
mat mat_mul(mat x,mat y)
{mat res;memset(res.a,0,sizeof(res.a));for(int i=0;i<2;i++)for(int j=0;j<2;j++)for(int k=0;k<2;k++)res.a[i][j]=(res.a[i][j]+x.a[i][k]*y.a[k][j])%MOD;return res;
}
void mat_pow(ll n)
{mat c,res;c.a[0][0]=c.a[0][1]=c.a[1][0]=1;c.a[1][1]=0;memset(res.a,0,sizeof(res.a));for(int i=0;i<2;i++) res.a[i][i]=1;while(n){if(n&1) res=mat_mul(res,c);c=mat_mul(c,c);n=n>>1;}printf("%I64d\n",res.a[0][1]);
}
int main()
{ll n;scanf("%lld",&n);mat_pow(n);return 0;
}

?