hdu 4549(快速幂）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4549

思路：大牛思路：简单推导下有,  a, b, a^1*b^1, a^1*b^2 .... 可以知道a,b的幂满足Fib, 然后构造矩阵快速幂...就好了.还需要个性质: A^X = A^( X mod Eular(M) ) ( mod M ) .

然后我就郁闷啦，自己的代码怎么都过不了，orz...

一下是wa代码，哪个大神帮忙看看：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
#define MOD 1000000007
typedef long long ll;
struct Matrix{
    ll map[2][2];
};
Matrix Mata,Unit;


void Initaite(){
    Unit.map[0][0]=Unit.map[1][1]=1;
    Unit.map[0][1]=Unit.map[1][0]=0;

    Mata.map[0][0]=0;
    Mata.map[0][1]=1;
    Mata.map[1][0]=1;
    Mata.map[1][1]=1;
}

Matrix Mul(const Matrix &p,const Matrix &q){
    Matrix r;
    memset(r.map,0,sizeof(r.map));
    for(int i=0;i<2;i++){
        for(int j=0;j<2;j++){
            for(int k=0;k<2;k++){
                r.map[i][j]=(r.map[i][j]+p.map[i][k]*q.map[k][j])%(MOD-1);
            }
        }
    }
    return r;
}

Matrix Pow(int n){
    Matrix q=Unit,p=Mata;
    while(n){
        if(n&1){
            q=Mul(p,q);
        }
        n>>=1;
        p=Mul(p,p);
    }
    return q;
}

ll PPow(ll x,ll n){
    ll p=x,q=1;
    while(n){
        if(n&1){
            q=p*q%MOD;
        }
        n>>=1;
        p=p*p%MOD;
    }
    return q;
}


int main(){
    Initaite();
    int a,b,n;
    while(~scanf("%d%d%d",&a,&b,&n)){
        if(n==0)printf("%d\n",a);
        else if(n==1)printf("%d\n",b);
        else {
            Matrix Mat1=Pow(n-1);
            ll x1=Mat1.map[1][1],x2=Mat1.map[0][1];
        //    cout<<x1<<x2<<endl;
            ll ans=PPow(a,x1)*PPow(b,x2)%MOD;
           printf("%d\n",(int)ans);
        }
    }
    return 0;
}