long long quickpow(long long m,long long n)
{
long long ans=1;
while(n)
{
if(n&1)
ans=(ans%mod)*(m%mod)%mod;
n=n>>1;
m=(m*m)%mod;
}
return ans;
}
mod是要求的取mod值。至于矩阵的高速幂就把乘法和计算的对象换为矩阵的即可了。
原文：http://www.cnblogs.com/yangykaifa/p/6880426.html