題面
大意就是求斐波那契數列第n項，做法為矩陣快速冪。

代碼

#include<iostream> #include<cstdio> #include<cstring> #define LL long longusing namespace std; const int mod = 10000; struct Mat{LL a[4][4];Mat(){memset(a,0,sizeof(a));}Mat operator*(const Mat &h){Mat c;for(register int i=1;i<=2;i++)for(register int j=1;j<=2;j++)for(register int k=1;k<=2;k++){c.a[i][j]+=a[i][k]%mod*h.a[k][j]%mod;c.a[i][j]%=mod;}return c;} }ans,f;int n;inline void fast_pow(Mat b,int k){for(;k;k>>=1){if(k&1) f=f*b;b=b*b;} }int main(){while(~scanf("%d",&n)){if(n==-1) break;f.a[1][1]=0;f.a[1][2]=1;ans.a[2][1]=ans.a[1][2]=ans.a[2][2]=1;ans.a[1][1]=0;fast_pow(ans,n);printf("%lld\n",f.a[1][1]);}return 0; }

轉載于:https://www.cnblogs.com/sdfzsyq/p/9677007.html

總結

以上是生活随笔為你收集整理的poj 3070的全部內容，希望文章能夠幫你解決所遇到的問題。

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