poj 3070
来源:互联网 发布:高价回收域名的事例 编辑:程序博客网 时间:2024/05/01 22:10
http://poj.org/problem?id=3070
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
099999999991000000000-1
Sample Output
0346266875
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 <stdio.h>#include <iostream>#include <string.h>#include <math.h>using namespace std;typedef long long LL;const int N=2;const LL MOD=10000;struct Matrix{ LL m[N][N];};Matrix I={ 1,0, 0,1,};Matrix multi(Matrix a,Matrix b){ Matrix c; for(int i=0; i<N; i++) for(int j=0; j<N; j++) { c.m[i][j]=0; for(int k=0; k<N; k++) { c.m[i][j]+=a.m[i][k]*b.m[k][j]%MOD; } c.m[i][j]=c.m[i][j]%MOD; } return c;}Matrix quick_mod(Matrix a,LL k){ Matrix ans=I; while(k!=0) { if(k&1) { ans=multi(ans,a); } k>>=1; a=multi(a,a); } return ans;}int main(){ LL n; while(~scanf("%lld",&n)) { if(n==-1) break; Matrix A= {1,1, 1,0 }; Matrix x1=quick_mod(A,n); printf("%lld\n",x1.m[0][1]); } return 0;}
- POJ 3070
- poj 3070
- POJ 3070
- poj 3070
- poj 3070
- POJ 3070
- poj 3070
- poj 3070
- (3070)POJ
- poj 3070
- poj 3070
- POJ 3070
- poj 3070
- poj 3070
- POJ 3070
- 【POJ】3070
- poj 3070 矩阵乘法
- POJ 3070 Fibonacci
- 沃尔法士大夫
- 除VISIO外,开源免费的流程图设计软件
- HDU 1864 最大报销额
- bzoj1024: [SCOI2009]生日快乐 DFS
- 只是做点笔记
- poj 3070
- apache2.2+django1.6.5+python34+mod_wsgi3.5环境搭建
- 在fragment中调用SharedPreferences
- 所得征收税费计算方式发掘
- JSP简单练习-javaBean的简单应用
- 增加VMware虚拟硬盘文件(.vmdk) 大小
- 青春,死亡的盛夏
- nodejs中npm常用命令
- [C++]_[获取Utf8字符串的字符个数和子字符串]