【洛谷1962】 斐波那契数列
2017-07-17 19:55
295 查看
题面
题目背景
大家都知道,斐波那契数列是满足如下性质的一个数列:• f(1) = 1
• f(2) = 1
• f(n) = f(n-1) + f(n-2) (n ≥ 2 且 n 为整数)
题目描述
请你求出 f(n) mod 1000000007 的值。输入格式:
·第 1 行:一个整数 n输出格式:
第 1 行: f(n) mod 1000000007 的值输入输出样例
输入样例#1:
5输出样例#1:
5输入样例#2:
10输出样例#2:
55说明
对于 60% 的数据: n ≤ 92对于 100% 的数据: n在long long(INT64)范围内。
题解
看一看数据范围如果使用O(n)的递推显然会炸掉
那么我们有没有别的方法?
显然是有的
使用斐波那契数列的递推公式怎么样?
但是,,,里面带有根号,如果直接使用显然是会掉精度的
所以,,,应该怎么办
我们知道
f[i]=f[i-1]+f[i-2]
f[i-1]=f[i-2]+f[i-3]
所以
我们可以用矩阵来表示
因此
我们可以继续推导
可以得到
接下来使用矩阵快速幂就可以直接求解
#include<iostream> #include<cstdlib> #include<cstdio> #include<cstring> #include<cmath> #include<algorithm> using namespace std; #define MOD 1000000007 #define MAX 10 #define ll long long struct yl//矩阵 { int n;//大小 long long g[MAX][MAX]; }; yl operator *(yl a,yl b)//定义乘法 { int n=a.n; yl cool; memset(cool.g,0,sizeof(cool.g)); for(int i=1;i<=n;++i) for(int j=1;j<=n;++j) for(int k=1;k<=n;++k) cool.g[i][j]=(cool.g[i][j]+1ll*a.g[i][k]*b.g[k][j]%MOD)%MOD; cool.n=n; return cool; } void write(yl a) { int n=a.n; for(int i=1;i<=n;++i) { for(int j=1;j<=n;++j) cout<<a.g[i][j]<<' '; cout<<endl; } } yl Pow(yl a,long long b)//a的b次方 { if(b==1)return a; yl s=Pow(a,b/2); s.n=a.n; s=s*s; if(b&1)s=s*a; return s; } int main() { ll n; cin>>n; if(n==0) { cout<<0<<endl; return 0; } if(n==1||n==2) { cout<<1<<endl; return 0; } else { yl a; a.n=2; a.g[1][1]=a.g[1][2]=a.g[2][1]=1; a.g[2][2]=0; yl s=Pow(a,n-1); s.n=2; //write(s); cout<<s.g[1][1]<<endl; return 0; } }
相关文章推荐
- 【洛谷1962】 斐波那契数列
- [洛谷1962]斐波那契数列
- 洛谷P1962 斐波那契数列
- 洛谷P1962 斐波那契数列
- 洛谷P1962 斐波那契数列
- 洛谷P1962 斐波那契数列【矩阵运算】
- LuoguP1962 斐波那契数列 解题报告【矩阵快速幂】
- 洛谷P1962 斐波那契数列 || P1349 广义斐波那契数列[矩阵乘法]
- 洛谷 P2626 斐波那契数列(升级版)
- 洛谷—— P1962 斐波那契数列
- 矩阵快速幂模板及其应用(洛谷3390)(洛谷1962)
- Luogu 1962 斐波那契数列(矩阵,递推)
- AC日记——斐波那契数列 洛谷 P1962
- 【洛谷P1962 斐波那契数列】矩阵快速幂+数学推导
- [luoguP1962] 斐波那契数列(矩阵快速幂)
- 洛谷P1962 斐波那契数列
- 洛谷_P1962_斐波那契数列
- 洛谷 P1962 斐波那契数列
- AC日记——斐波那契数列(升级版) 洛谷 P2626
- 洛谷 P2626 斐波那契数列(升级版)