HDU - 3117 Fibonacci Numbers 矩阵快速幂 + 取大数前4位
2015-06-02 00:19
387 查看
题目大意:要求输出第n个fibonacci数,如果该数超过1e9,就输出该数的前4位和后四位
解题思路:通过打表可得,第40个fibonacci数是大于1e9的,所以40之前的可以直接计算
40之后的比较麻烦,参考了别人的题解
http://blog.sina.com.cn/s/blog_9bf748f301019q3t.html
解题思路:通过打表可得,第40个fibonacci数是大于1e9的,所以40之前的可以直接计算
40之后的比较麻烦,参考了别人的题解
http://blog.sina.com.cn/s/blog_9bf748f301019q3t.html
[code]#include<cstdio> #include<cmath> using namespace std; typedef long long ll; const int N = 2; struct Matrix{ int mat ; }A, B, tmp; int n; Matrix matMul(Matrix x, Matrix y, int mod) { for(int i = 0; i < N; i++) for(int j = 0; j < N; j++) { tmp.mat[i][j] = 0; for(int k = 0; k < N; k++) tmp.mat[i][j] += x.mat[i][k] * y.mat[k][j] % mod; tmp.mat[i][j] %= mod; } return tmp; } void solve(int mod) { while(n) { if(n & 1) B = matMul(B,A,mod); A = matMul(A,A,mod); n >>= 1; } } void init() { B.mat[0][0] = B.mat[1][1] = 1; A.mat[0][0] = A.mat[0][1] = A.mat[1][0] = 1; B.mat[0][1] = B.mat[1][0] = A.mat[1][1] = 0; } int main() { while(scanf("%d", &n) != EOF) { init(); if(n == 0) printf("0\n"); else if(n == 1) printf("1\n"); else { init(); if(n <= 39) { n -= 1; solve(1e9); printf("%d\n", B.mat[0][0]); } else { double t; double s = (sqrt(5.0) + 1.0) / 2; t = -0.5 * log(5.0) / log(10.0) + ((double)n) * log(s)/ log(10.0); t -= floor(t); t = pow(10.0,t); while(t < 1000) t *= 10; n--; solve(10000); printf("%04d...%04d\n", int(t), B.mat[0][0]); } } } return 0; }
相关文章推荐
- CentOS6.6默认安装完虚机化平台和虚拟化工具后没有virt-install命令
- 音视频播放器源码
- 第三章 布局
- 在阿里云centos6.5系统上用docker搭建php环境的失败过程
- Redis配置文件参数说明
- 初稿-学习Photoshop创建一个下载应用界面实例
- poj3695&bzoj1853 容斥定理
- jps命令使用
- https
- 6.1学习内容 ndk编程释放内存,C语言中char*添加字符得方法,java二维数组的读取与生成
- iOS NSFileHandle 基本用法介绍
- 总结SSH框架整合基础与需要的jar包
- 像写SQL一样编写Java数据应用-TinySqlDsl
- iOS网络编程开发笔记1—JSON和XML数据解析
- 企业库一全套增删改
- Java-Fork/Join的简单例子
- 甘特图与网络图
- JavaScript Cookies
- [转载]nil,Nil,NULL 和NSNull的小结
- Python Show-Me-the-Code 第 0021题 密码加密