快速幂与快速乘(c语言实现)

快速幂思想

快速幂是什么？
简单来说，快速幂就是快速算出a的b次方模c的一种算法。

那么，它的本质是什么？
本质上是一个数学公式

a^b % c = （a % c） * （a % c） * … *（a % c） % c（总共b个小括号） 
若b为偶数 = （a^2 % c) * (a^2 % c) * … * (a^2 % c) % c （总共b/2个小括号）
若b为奇数 = （a^2 % c) * (a^2 % c) * … * (a^2 % c)  * a % c % c（总共（b-1）/ 2个小括号） 
……
然后一直指数减半，底数平方取模。最后得出答案。

在快速幂运算当中，每一次运算我们都能保证答案不会爆掉c，而且可以极大的减少运算时间，快速幂也被叫做蒙哥马利算法，蒙哥马利算法使得数据规模极大减小。使原本看似不可能完成的程序变得可行。

c语言实现

#include<stdio.h>
long long a;//底数
long long b;//指数
long long c;//取模数
long long Montgomery()
{
    long long int res = 1;
    while(b)
    {
        if (b%2==1)//如果为奇数
            res = (res*a) % c;
        b/=2;//指数对半
        a = (a*a) % c;//底数平方
    }
    return res;
}
int main()
{
    scanf("%lld%lld%lld", &a, &b, &c);
    printf("%lld", Montgomery());
    return 0;
}

那快速幂中a * a也会爆怎么办？

那就用快速乘，快速乘其实就是快速幂的一种变种，与快速幂同理

a*53 = a*110101（二进制）= a*（100000+10000+100+1）=a*（100000*1+10000*1+1000*0+100*1+10*0+1*1）。

我们把乘数拆分成2进制，乘数每向前移一位，a就*=2，相当于加自己本身，这时候就与快速幂同理了

当然，因为你要进位那些，所以用了快速乘会慢一些，但该用的地方还是要用

直接上代码

#include<stdio.h>
long long a;//底数
long long b;//指数
long long c;//取模数
long long ksc(long long x,long long y){
	long long ans=0;
	while(y){
		if(y%2==1){
			ans+=x;
			ans%=c;
		}
		x+=x;   //x*=2
		x%=c;
		y/=2; //乘数向前移 
	}
	return ans;
}
long long ksm()
{
    long long int res = 1;
    while(b)
    {
        if (b%2==1){//如果为奇数
        	res = ksc(res,a);
		}         
        b/=2;//指数对半
        a = ksc(a,a);//底数平方
    }
    return res;
}
int main()
{
    scanf("%lld%lld%lld", &a, &b, &c);
    printf("%lld",ksm());
    return 0;
}