快速幂及矩阵快速幂模板

快速幂是基于减少运算次数来提高运算速度的方法

普通快速幂

//求x的n次幂 
int QuickPow(int x,int N)
{
	int res=x;
	int ans=1;
	while(N)
	{
		if(N&1)
		{
			ans=ans*res;
		}
		res=res*res;
		N=N>>1;
	}
	return ans;
}

矩阵快速幂：

为了便于理解快速幂是如何和矩阵结合起来了，下面就斐波那契数列进行举例说明

斐波那契数列的表达式：f[i]=f[i-1]+f[i-2];可表示为f[i]=1*f[i-1]+1*f[i-2];

                                       我们把f[i-1]表示成f[i-1]=1*f[i-1]+0*f[i-2];

                                        用矩阵来表示上面两组关系

 

下面模板

//结构矩阵 
struct Matrix
{
	int m[maxn][maxn];
}ans,res;
//矩阵乘法 
Matrix Mul(Matrix a,Matrix b,int n)
{
	Maxtrix temp;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			temp.m[i][j]=0;
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(j=1;j<=n;j++)
		{
			for(k=1;k<=n;k++)
			{
				temp.m[i][j]=a.m[i][j]*b.m[k][j];
			}
		}
	}
	return temp;
}
//矩阵快速幂，求矩阵res的n次幂
void quickpower(int N,int n)
{
	//整数快速幂默认的ans是1，矩阵的话ans应为单位矩阵
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			if(i==j)
			ans.m[i][j]=1;
			else
			ans.m[i][j]=0;
		}
	 }
	 while(N)
	 {
	 	if(N&1)
	 	ans=Mul(res,res);
	 	res=Mul(res,res);
	 	N=N>>1;
	  } 
 }