快速幂求模+矩阵的快速幂

int quickMod(int a,int b)
{
	int ans = 1;
	while (b)
	{
		if (b&1)
			ans = ans * a % MOD;
		a = a*a % MOD;
		b >>= 1;
	}
	return ans;
}

矩阵的快速幂模板

#include<stdio.h>  
#include<string.h>  
#include<algorithm>  
using namespace std;  
struct Mat//定义矩阵  
{  
    int a[2][2];  
    void init()  
    {  
        memset(a,0,sizeof(a));  
        a[0][0]=a[1][1]=1;  
     }   
};  
Mat mul(Mat a,Mat b)//矩阵乘法   
{  
    Mat ans;  
    ans.init();  
    for(int i=0;i<2;i++)  
    {  
        for(int j=0;j<2;j++)  
        {  
            ans.a[i][j]=0;  
            for(int k=0;k<2;k++)  
            {  
                ans.a[i][j]+=a.a[i][k]*b.a[k][j];  
            }  
            ans.a[i][j]%=10000;  
        }  
    }  
    return ans;  
}  
Mat power(Mat a,int num)//矩阵快速幂模板   
{  
    Mat ans;  
    ans.init();  
    while(num)  
    {  
        if(num%2==1)  
        {  
            ans=mul(ans,a);  
        }  
        num/=2;  
        a=mul(a,a);  
    }  
    return ans;  
}