个人模板 快速幂取模，矩阵快速幂，快速乘

#include<bits/stdc++.h>
using namespace std;
const int MOD = 1007;
int xx(int a,int n,int MOD) //计算a的n次方对MOD取模
{
    int ret=1;
    int tmp=a%MOD;
    while(n)
    {
        if(n&1)
            ret=ret*tmp%MOD;
        tmp=tmp*tmp%MOD;
        n>>=1;
    }
    return ret;
}
int main()
{
    int m,n;
    while(scanf("%d%d",&m,&n)==2)
    {
        printf("%d\n",xx(m,n,MOD));
    }
}

#include <bits/stdc++.h>
using namespace std;
const int N = 4 ;
int n;
struct xx
{
    int a[N][N];
} ori,res;
xx mul(xx x,xx y)
{
    xx temp;
    memset(temp.a,0,sizeof(temp.a));
    for(int i=0; i<N; i++)
    {
        for(int j=0; j<N; j++)
        {
            for(int k=0; k<N; k++)
            {
                temp.a[i][j]+=x.a[i][k]*y.a[k][j];
            }
        }
    }
    return temp;
}
void init()
{
    int x;
    for(int i=0; i<N; i++)
    {
        for(int j=0; j<N; j++)
        {
            scanf("%d",&x);
            ori.a[i][j]=x;
        }
    }
    for(int i=0; i<N; i++)
    {
        for(int j=0; j<N; j++)
        {
            res.a[i][j]=i==j?1:0;
        }
    }
}
xx calc(long long k)                                          
{
    while(k)
    {
        if(k%2==1)                      ///这么写比普通的快速幂好，亲身体会，HDU 5667用普通写法 TLE，这么写就过了
        {
            res=mul(ori,res);         
            k-=1;
        }
        else
        {
            ori=mul(ori,ori);
            k/=2;
        }
    }
    return res.a[0][0];
}
int main()
{
    while(cin>>n)
    {
        init();
        calc(n);
    }
}

long long mul(long long a,long long b)  
{  
    long long res=0;  
    while(b)  
    {  
        if(b&1) res=(a+res)%Mod;  
        b>>=1;  
        a=(a+a)%Mod;  
    }  
    return res;  
}