快速幂

快速幂

快速幂的主要作用是求a^b对mod取模，尤其适合于指数b和mod很大的情况，它的原理也十分简单，我们举一个简单的例子：2¹⁰%100

如果用循环乘法来实现，一共要执行10次乘法，然后再对100取余数；

而快速幂的思想则是减少执行的次数，把2¹⁰变成（2*2）⁵，这样做的好处是用一次底数的平方，减少了一半的执行量，指数b越大，减少的操作越多。

我们继续，4⁵还能不能继续这样变换呢？当然可以，但是因为5是奇数，我们不能将它直接除以2，但是可以提出一个4¹，这样我们又可以对剩下的4⁴进行操作了，一直继续下去：

4¹*4⁴=4¹*16²=4¹*256¹

对于取模的问题，我们都知道，(a*b)%mod=((a%mod)*(b%mod))%mod，因此在每次乘以底数a的时候都取一次模即可。

然后我们用流程图来表示一下整个过程：

然后只要在b>0的时候循环这个流程即可。

代码：

long long quick_pow(long long a,long long b,long long mod){
	long long ans=1;
	while(b>0){
		if(b&1){//b为奇数
			ans*=a%mod;
		}
		b>>=1;//b除以2
		a=(a*a)%mod;
	}
	return ans;
}

其中b&1相当于b%2==1；b>>=1相当于b/=2。

这里是运用了位运算，在C/C++中相比于逻辑运算执行速度更快。

快速乘法
根据快速幂的原理我们可以得出快速乘法的流程。在快速幂中我们是把幂运算降为乘法运算，那么快速乘法只需要将乘法运算降为加法运算即可。
代码：

long long quick_mult(long long a,long long b,long long mod){
	long long ans=0;
	while(b>0){
		if(b&1){//b为奇数
			ans+=a%mod;
			b--; //可以不写，因为C++整数除法向下取整
		}
		b>>=1;//b除以2
		a=(a+a)%mod;
	}
	return ans;
}