【hdu 4549】 M斐波那契数列

Problem Description

M斐波那契数列F[n]是一种整数数列，它的定义如下：

F[0] = a
F[1] = b
F[n] = F[n-1] * F[n-2] ( n > 1 )

现在给出a, b, n，你能求出F[n]的值吗？

 

Input

输入包含多组测试数据；
每组数据占一行，包含3个整数a, b, n（ 0 <= a, b, n <= 10^9 ）

 

Output

对每组测试数据请输出一个整数F[n]，由于F[n]可能很大，你只需输出F[n]对1000000007取模后的值即可，每组数据输出一行。

 

Sample Input

0 1 0

6 10 2

 

Sample Output

0

60

 

Source

2013金山西山居创意游戏程序挑战赛——初赛（2）

 

这道题是一个矩阵快速幂+费马小定理降幂，首先设Fibonacci数列第n项为fib[n],可以推出F[n]=a^fib[n-1]*b^fib[n]，然后用矩阵快速幂求fib[n-1]和fib[n]，再根据费马小定理得到当a,p互质时a^b%p=a^(b%p)%p，这样就可以直接快速幂解决了，下面是程序：

#include<stdio.h>
#include<iostream>
#define ll long long
using namespace std;
struct node{
	ll x[2][2];
	node(){
		x[0][0]=0;
		x[0][1]=x[1][0]=x[1][1]=1;
	}
};
node mul(node &a,node &b,ll &mod){
	node s;
	int i,j,k;
	for(i=0;i<2;i++){
		for(j=0;j<2;j++){
			s.x[i][j]=0;
			for(k=0;k<2;k++){
				s.x[i][j]+=a.x[i][k]*b.x[k][j];
				s.x[i][j]%=mod;
			}
		}
	}
	return s;
}
node fib(ll n,ll mod){
	node s,r;
	s.x[0][0]=s.x[1][1]=1;
	s.x[1][0]=s.x[0][1]=0;
	while(n){
		if(n&1){
			s=mul(s,r,mod);
		}
		r=mul(r,r,mod);
		n>>=1;
	}
	return s;
}
ll pow(ll a,ll b,ll p){
	ll s=1,r=a;
	while(b){
		if(b&1){
			s*=r;
			s%=p;
		}
		r*=r;
		r%=p;
		b>>=1;
	}
	return s;
}
int main(){
	ll a,b,n;
	const ll mod=1000000007;
	while(~scanf("%lld%lld%lld",&a,&b,&n)){
		node s=fib(n,mod-1);
		printf("%lld\n",(pow(a,s.x[0][0],mod)*pow(b,s.x[0][1],mod))%mod);
	}
	return 0;
}