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

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; }

?

轉載于:https://www.cnblogs.com/aerer/p/9930934.html

總結

以上是生活随笔為你收集整理的矩阵快速幂及斐波那契数列模板的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。