2016算法第一次练习赛——C 斐波那契进阶
来源:互联网 发布:怎么看电脑端口号 编辑:程序博客网 时间:2024/06/06 11:42
C 斐波那契进阶
时间限制:1000ms 内存限制:65536kb
题目描述
AlvinZH听说大家都掌握了斐波那契数列,先简单了解一下,斐波那契整数序列中,有F(0)=0,F(1)=1,F(n)=F(n-1)+F(n-2) (n≥2)。这样我们可以快速地求出斐波那契数列的前几项:0,1,1,2,3,5,8,13,21,34,...
你知道吗?斐波那契序列还有另外一种定义:
AlvinZH不满足于此,他想知道当n很大很大时,F(n)有多大,你来告诉他吧!(结果对10007取模)
输入
输入将由多组测试数据组成,以EOF结尾。
每组数据只含一行,为题中描述的整数n(0≤n≤10^9)。
输出
对于每组数据,输出一行,为F(n)对10007取模的值。
输入样例
011000000000
输出样例
017300
HINT
2×2矩阵的零次幂为单位矩阵(如图1),计算两个2×2矩阵的乘积如图2。
题目分析
f[i-1]=f[i-1]
f[i]=f[i-1]+f[i-2]
可能你会觉得第一句有一点废话,但是是有不一样的意义的,这体现了递推的过程,就是说我们每次扫两个数,根据前面的两个数推出后面的两个数,根据这个我们可以建一个2*2的矩阵A
1110
然后把我们每次存的两个数放在另一个矩阵B里面:
那么把这两个矩阵相乘就可以得到另一个矩阵:
这样就可以得到f[i]了,至于为什么乘了之后会变成这两个数,我们根据矩阵乘法的乘法规律可以很容易推出来。
这样子的话,我们每次用A*B替换B,最后得到的矩阵的第一个数就是f[n]了。
那么,矩阵乘法的优越性究竟体现在哪里呢。其实,矩阵乘法只是体现了我们从之前求的数到现在要求的数的递推过程,就是说矩阵乘法可以完成多个元素的递推。不过这个我们用普通的递推就可以实现的啊~~认真想想我们就能发现,我们在矩阵乘法的过程中把上见面的A矩阵自己相乘了很多遍。就是说,我们可以求A矩阵的幂最后乘上B矩阵,既然要求幂,矩阵乘法满足结合律,那么我们就可以用快速幂啦~~矩阵乘法的优越性就体现在这里:在递推过程变成不断乘以一个矩阵,然后用快速幂快速求得从第一个到第n个的递推式,这样子就可以在短时间内完成递推了。
快速幂写法:
while(n){ if(n&1) b=a*b; a=a*a; n>>=1;}
示例代码
#include<cstdio>#include<cstdlib>#include<cstring>#define Mod 10007struct node{ int v[3][3]; int m,l;};node get_mul(node a,node b){ node c; c.m=a.m;c.l=b.l; for(int i=1;i<=c.m;i++) for(int j=1;j<=c.l;j++) { c.v[i][j]=0; for(int k=1;k<=a.l;k++) c.v[i][j]=(c.v[i][j]+a.v[i][k]*b.v[k][j])%Mod; } return c;}int main(){ int n; while(scanf("%d",&n)!=EOF) { if(n==0) {printf("0\n");continue;} node a,b,c; a.m=a.l=2,a.v[1][1]=1,a.v[1][2]=1,a.v[2][1]=1,a.v[2][2]=0; b.m=b.l=2,b.v[1][1]=1,b.v[1][2]=0,b.v[2][1]=0,b.v[2][2]=1; c.m=2,c.l=1,c.v[1][1]=1,c.v[2][1]=0; n--; while(n) { if(n&1) b=get_mul(a,b); a=get_mul(a,a); n>>=1; } b=get_mul(b,c); printf("%d\n",b.v[1][1]); }}
- 2016算法第一次练习赛——C 斐波那契进阶
- 练习赛一 C 斐波那契进阶
- 算法练习---斐波那契数列
- 算法练习——”递推之斐波那契“
- 算法题目——python基础巩固练习一 斐波那契数列
- 算法练习——求斐波那契第n项
- C 练习实例11 斐波那契数列
- ACM第一次练习—1002C
- 2016算法第一次练习赛——D AlvinZH的儿时回忆——跳房子
- 2016算法第一次练习赛——A 群鸦的盛宴
- 2016算法第一次练习赛——B 朴素的中位数
- C语言之基本算法06—斐波那契数列求和
- 剑指offer-算法题练习:part7 斐波那契数列-斐波那契数列问题
- 斐波那契查找算法完整C代码
- 【C】斐波那契数列(递归算法)
- 斐波那契数列算法c语言实现
- 【C++练习】斐波那契数列
- Scala练习-斐波那契查找
- 【状压DP】【cofun1372】售货员难题
- Thymeleaf模板引擎的一些知识点
- Unity SteamVR插件集成
- 2.2.12多线程的死锁(jps命令、jstack命令)
- 【20171011】python_语言设计(6)程序设计方法与面向对象程序设计
- 2016算法第一次练习赛——C 斐波那契进阶
- MySQL索引类型总结和使用技巧以及注意事项
- C++STL之vector详解(转)
- Android 代码设置来电铃声
- Servlet登录传递注意事项
- 日记(周中)
- c语言 计算闰年
- html —— 查看元素属性 与 js 全局属性
- Android Killer反编译卡死的一个解决方案