POJ3070 Fibonacci(矩阵快速幂)
2016-01-19 20:18
465 查看
用矩阵快速幂求fibonacci序列的第n项。
/* *********************************************** Author :devil Created Time :2016/1/19 20:11:43 ************************************************ */ #include <iostream> #include <algorithm> #include <cstring> #include <cmath> #include <queue> #include <map> #include <set> #include <vector> #include <cstdio> using namespace std; const int MOD=1e4; struct matrix { int m[2][2]; } ans,base; matrix multi(matrix a, matrix b) { matrix tmp; for(int i=0; i<2; ++i) { for(int j=0; j<2; ++j) { tmp.m[i][j]=0; for(int k=0; k<2; ++k) tmp.m[i][j]=(tmp.m[i][j]+a.m[i][k]*b.m[k][j])%MOD; } } return tmp; } int fast_mod(int n)//求矩阵base的n次幂 { base.m[0][0]=base.m[0][1]=base.m[1][0]=1; base.m[1][1]=0; ans.m[0][0]=ans.m[1][1]=1;//ans初始化为单位矩阵 ans.m[0][1]=ans.m[1][0]=0; while(n) { if(n&1)//实现ans*=t;其中要先把ans赋值给tmp,然后用ans=tmp*t ans=multi(ans,base); base=multi(base,base); n>>=1; } return ans.m[0][1]; } int main() { //freopen("in.txt","r",stdin); int n; while(scanf("%d",&n)&&n!=-1) { printf("%d\n",fast_mod(n)); } return 0; }
相关文章推荐
- 松下蓄电池与UPS使用和维护
- C++的运算符重载
- MFC操作SQlite,打开数据库路径存在中文,解决方案
- android 4.0.3 设置时间无法更新到RTC
- foreach
- CSS的::selection使用方法
- Reactor模式
- 异步打开一个文件进行读取
- confluence安装
- POJ 3416 Sunscreen(贪心)
- 每个程序员需掌握的20个代码命名小贴士
- Single Number leetcode
- C++之智能指针(二)
- java 读取XML文件作为配置文件
- 视频采集处理之V4L2编程
- Android Studio导入GitHub第一个项目PullToRefresh
- leetcode 之valid anagram
- windows程序托盘图标以及消息提示--示例程序
- Java——BufferedReader和BufferedWriter
- C++二维指针内存开辟和释放