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

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;復(fù)制直接在multi里面實(shí)現(xiàn)了；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; }

?

轉(zhuǎn)載于:https://www.cnblogs.com/aerer/p/9930934.html

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。