POJ 3070 斐波那契数列(矩阵快速幂(版题))
来源:互联网 发布:iptv管理系统php 编辑:程序博客网 时间:2024/06/06 06:31
链接:http://poj.org/problem?id=3070
- Fibonacci
Time Limit: 1000MS Memory Limit: 65536K
Total Submissions: 16351 Accepted: 11470
Description
In the Fibonacci integer sequence, F0 = 0, F1 = 1, and Fn = Fn − 1 + Fn − 2 for n ≥ 2. For example, the first ten terms of the Fibonacci sequence are:
0, 1, 1, 2, 3, 5, 8, 13, 21, 34, …
An alternative formula for the Fibonacci sequence is
.
Given an integer n, your goal is to compute the last 4 digits of Fn.
Input
The input test file will contain multiple test cases. Each test case consists of a single line containing n (where 0 ≤ n ≤ 1,000,000,000). The end-of-file is denoted by a single line containing the number −1.
Output
For each test case, print the last four digits of Fn. If the last four digits of Fn are all zeros, print ‘0’; otherwise, omit any leading zeros (i.e., print Fn mod 10000).
Sample Input
0
9
999999999
1000000000
-1
Sample Output
0
34
626
6875
Hint
As a reminder, matrix multiplication is associative, and the product of two 2 × 2 matrices is given by
.
Also, note that raising any 2 × 2 matrix to the 0th power gives the identity matrix:
.
.
#include <iostream>#include <string.h>using namespace std;const int maxn =15;const int mod =10000;struct matrix{ long long int a[2][2];};matrix multply(matrix x,matrix y){ int i,j,k; matrix result; memset(result.a,0,sizeof(result.a));//** for (k=0; k<2; k++) for (i=0; i<2; i++) { if(x.a[i][k]) for (j=0; j<2; j++) { if (y.a[k][j]) { result.a[i][j]+=x.a[i][k]*y.a[k][j]; if (result.a[i][j]>=mod) result.a[i][j]%=mod; } } } return result;}void matrix_pow(matrix a,int nn){ matrix result; memset(result.a,0,sizeof(result.a)); int i,j,k; for (i=0; i<2; i++) result.a[i][i]=1; while (nn) { if (nn&1) result=multply(result,a); nn>>=1; a=multply(a,a);// int i,j;// for (i=0; i<2; i++)// for (j=0; j<2; j++)// if(!j)// cout<<a.a[i][j]<<" ";// else// cout <<a.a[i][j]<<endl; } cout<<result.a[0][1]<<endl;}int main (){ long long int nn; while (cin>>nn) { if (nn==-1) break; matrix a; a.a[0][0]=a.a[0][1]=a.a[1][0]=1; a.a[1][1]=0;// int i,j;// for (i=0; i<2; i++)// for (j=0; j<2; j++)// if(!j)// cout<<a.a[i][j]<<" ";// else// cout <<a.a[i][j]<<endl; matrix_pow(a,nn); } return 0;}
- POJ 3070 斐波那契数列(矩阵快速幂(版题))
- POJ - 3070 - Fibonacci (矩阵快速幂 + 斐波那契数列)
- poj 3070-Fibonacci (矩阵快速幂 求 斐波那契数列)
- 斐波那契数列(矩阵快速幂)
- *矩阵快速幂(斐波那契数列模板)
- 斐波那契数列 (矩阵快速幂)
- poj-3070 Fibonacci(矩阵快速幂 + 斐波那契数列)
- POJ 3070 (矩阵快速幂 ——斐波那契数列)
- 矩阵快速幂优化递推式 例:斐波那契数列 poj 3070
- 【矩阵快速幂】HDU 4549 : M斐波那契数列(矩阵嵌套)
- 快速幂、矩阵快速幂+斐波那契数列(Fibonacci Sequence)
- hdu 4549 M斐波那契数列(矩阵快速幂,快速幂降幂)
- POJ 3070 矩阵快速幂求斐波那契数列
- POJ 3070-Fibonacci(矩阵快速幂求斐波那契数列)
- hdoj 4549 M斐波那契数列(费马小定理降幂,矩阵快速幂)
- HDU 4549 M斐波那契数列(矩阵快速幂)(费马小定理)
- HDOJ 题目4549 M斐波那契数列(数学,矩阵快速幂,费马小定理)
- HDU - 4549 M斐波那契数列(矩阵快速幂+费马小定理)
- CodeForces Round#433
- 一道面试题
- [JZOJ5363]生命之树
- inline函数必须定义在头文件吗
- Codeblocks environment error编译器找不到问题解决方法
- POJ 3070 斐波那契数列(矩阵快速幂(版题))
- 【算法】图的最短路径(Floyd算法)
- org.hibernate.HibernateException: JDBC Driver class not found: com.mysql.jdb
- 心路里程碑
- 微信小程序UI------实现携程首页顶部的界面(弹性布局)
- Qt之模型/视图
- OkHttp
- Leetcode538. 采用中序遍历将BST转化为节点值更大的树
- css 中的属性