又见斐波那契数列
时间限制:
1000 ms | 内存限制:
65535 KB
难度:4
-
描述
-
斐波那契数列大家应该很熟悉了吧。下面给大家引入一种新的斐波那契数列:M斐波那契数列。 M斐波那契数列F[n]是一种整数数列,它的定义如下:
F[0] = a
F[1] = b
F[n] = F[n-1] * F[n-2] ( n > 1 )
现在给出a, b, n,聪明的你能求出F[n]的值吗?-
输入
-
输入包含多组测试数据;
每组数据占一行,包含3个整数a, b, n( 0 <= a, b, n <= 10^9 )
输出
- 对每组测试数据请输出一个整数F[n],由于F[n]可能很大,你只需输出F[n]对1000000007取模后的值即可,每组数据输出一行。 样例输入
-
0 1 0 6 10 2
样例输出
-
0 60
-
输入包含多组测试数据;
根据列举猜想出公式,然后用矩阵快速幂求出指数,然后再用快速幂求出答案
摸除技巧:
(a+b)%mod=((a%mod)+(b%mod))%mod
(a-b)%mod=((a%mod)-(b%mod)+mod)%mod
ab%mod=(a%mod)(b%mod)%mod
a^b%mod=a^(b%(mod-1))%mod // mod 是一个质数,该式子称费马小定理
代码如下:
#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;
typedef long long ll;
const int mod=1000000007;
struct Matrix{
ll matrix[2][2];
};
int n=2;
void init(Matrix &res)
{
memset(res.matrix,0,sizeof(res.matrix));
for(int i=0;i<n;i++)
res.matrix[i][i]=1;
}
Matrix multiplicative(Matrix a,Matrix b)
{
Matrix res;
memset(res.matrix,0,sizeof(res.matrix));
for(int i = 0 ; i < n ; i++){
for(int j = 0 ; j < n ; j++){
for(int k = 0 ; k < n ; k++){
res.matrix[i][j]=(res.matrix[i][j]%(mod-1)+ a.matrix[i][k]*b.matrix[k][j]%(mod-1))%(mod-1);
}
}
}
return res;
}
Matrix powm(Matrix mx,ll m)
{
Matrix res,base=mx;
init(res);
while(m)
{
if(m&1)
res=multiplicative(res,base);
base=multiplicative(base,base);
m>>=1;
}
return res;
}
ll pow(ll a,ll b)
{
ll res=1;
while(b)
{
if(b&1)
res=res*a%mod;
a=a*a%mod;
b>>=1;
}
return res;
}
int main()
{
ll a,b,m;
while(~scanf("%lld%lld%lld",&a,&b,&m))
{
if(m==0)
printf("%lld\n",a%mod);
else if(m==1)
printf("%lld\n",b%mod);
else
{
ll fn1,fn;
Matrix ans,res;
ans.matrix[0][0]=ans.matrix[0][1]=ans.matrix[1][0]=1;
ans.matrix[1][1]=0;
res=powm(ans,m-1);
fn1=res.matrix[0][1];
fn=res.matrix[0][0];
ll anc=(pow(a,fn1)%mod*pow(b,fn)%mod)%mod;
printf("%lld\n",anc);
}
}
return 0;
}