poj3070 Fibonacci(矩阵乘法初学)
来源:互联网 发布:网络学历教育报名 编辑:程序博客网 时间:2024/05/22 13:32
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:
.
Source
[Submit] [Go Back] [Status] [Discuss]
代码:
#include<iostream>#include<cstdio>#include<cstring>using namespace std;#define mod 10000struct mat{ int A[3][3];}s;mat matmul(mat x,mat y) //x*y {mat t; t.A[0][0]=(x.A[0][0]*y.A[0][0]+x.A[0][1]*y.A[1][0])%mod; t.A[0][1]=(x.A[0][0]*y.A[0][1]+x.A[0][1]*y.A[1][1])%mod; t.A[1][0]=(x.A[1][0]*y.A[0][0]+x.A[1][1]*y.A[1][0])%mod; t.A[1][1]=(x.A[1][0]*y.A[0][1]+x.A[1][1]*y.A[1][1])%mod; return t;}mat matpow(mat ss,int k) //ss^k{if(k==1) return s;ss=matpow(ss,k/2);if(k%2==0){return matmul(ss,ss);}elsereturn matmul(matmul(ss,ss),s);}int main(){int i,j,n;while(~scanf("%d",&n)&&n!=-1){ if(n==0){printf("0\n");continue;}s.A[0][0]=s.A[0][1]=s.A[1][0]=1;s.A[1][1]=0; s=matpow(s,n);printf("%d\n",s.A[0][1]);}return 0;}
- poj3070 Fibonacci(矩阵乘法初学)
- 【poj3070】 Fibonacci 【矩阵乘法】
- 【poj3070】Fibonacci(矩阵)
- 【poj3070】Fibonacci (矩阵乘法+快速幂)
- poj3070(矩阵快速幂,矩阵乘法)
- POJ3070 Fibonacci (矩阵连乘)
- POJ3070 - Fibonacci (矩阵快速幂)
- poj3070 Fibonacci (矩阵快速幂)
- [POJ3070]Fibonacci(矩阵快速幂)
- poj3070 Fibonacci(矩阵快速幂)
- poj3070 Fibonacci(矩阵快速幂)
- poj3070 Fibonacci(矩阵倍增模板)
- poj3070 矩阵乘法求斐波拉契
- poj3070 Fibonacci 矩阵快速幂
- poj3070 Fibonacci 矩阵快速幂
- POJ3070 - Fibonacci - 矩阵快速幂
- POJ3070 Fibonacci(矩阵快速幂)
- poj3070 Fibonacci(矩阵快速幂)
- 11b 一个简单实用好学易于理解的模态对话框 附截图
- 唯有自由人才能创造出金字塔
- 白话数据结构之【线段树】
- windows 资源保护找到了损坏文件但无法修复其中某些文件
- 《数据结构编程实验》 1.4.1Hangover
- poj3070 Fibonacci(矩阵乘法初学)
- Palindromic Permutation - UVa 11027 组合数
- php substr()根据起点和终点访问字符串
- 考试排名
- oracle的购买价格研究
- Android Studio 错误: 非法字符: '\ufeff' 解决方案|错误: 需要class, interface或enum
- Food (hdu 4292 网络流sap模板题)
- CreateJS基础 学习笔记(上)
- 从头到尾彻底理解KMP