M斐波那契数列

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

6 10 2

Sample Output

60

#include<stdio.h>
#include<string.h>
struct Edge
{
    long long num[2][2];
}x;
int a,b,n;
Edge play(Edge p,Edge y)
{
    Edge idx;
    memset(&idx,0,sizeof(idx));
    for(int i=0;i<2;i++)
        for(int j=0;j<2;j++)
            for(int k=0;k<2;k++)
                idx.num[i][j]=(idx.num[i][j]+p.num[i][k]*y.num[k][j])%1000000006;
    return idx;
}
Edge power(Edge k,long long y)
{
if(y==1)
return k;
Edge idx=power(k,y/2);
idx=play(idx,idx);
if(y%2)
idx=play(idx,k);
return idx;
}
long long ksm(long long k,long long y)
{
    long long ans=1;
    while(y)
    {
        if(y%2)
        ans=ans*k%1000000007;
        k=k*k%1000000007;
        y/=2;
    }
    return ans;
}
int main()
{
    scanf("%d%d%d",&a,&b,&n);
x.num[0][0]=0;
    x.num[0][1]=1;
x.num[1][0]=1;
x.num[1][1]=1;
    Edge idx=power(x,n);
    int ans=(ksm(a,idx.num[0][0])*ksm(b,idx.num[1][0]))%1000000007;
    printf("%d",ans);
}