快速幂2020/10/20

反反复复搞这个快速幂

快速幂和普通幂不同在于：

1、求幂方法

普通a^b=aa……a
此方法a^b=((aa)(aa)(aa)(a*a))…
eg:3^5= (3^2) ^2 *3^1
O(n)=>O(logn)

2、二进制运算

判断奇偶性不用if(a/2),而用if(a&1);
乘2除2变成a=a>>1,a=a<<1(二进制位运算，二进制向左/右移一格，即乘2或除2);

3、在只要求求幂的余数时

由于a为底数，n为指数，p为余数，
则(a^ n)%p=((a^ b)%p* (a^(n-b))%p)%p
例：2^40除以13的余数，1212要比512512快吧。
这题其实我这些点都在其他地方学到了，真正折磨我的不是这些，而是以下几点：
1、幂运算放在一个独立函数里时最后一个节点过不了；
2、改到主函数里时有几个变量没改好导致结果完全错了（我怎么会这么粗心，唉）。
AC代码

#include<iostream>
using namespace std;
int main()
{
	long long b,p,k,s;
	cin>>b>>p>>k;
	long long result=1,B=b,P=p;
	b%=k;
	while(p>0)
	{
		if(p&1)
			result=(result*b)%k;
		b=(b*b)%k;
		p=p>>1;
	}
	result=result%k;
	cout<<B<<"^"<<P<<" mod "<<k<<"="<<result;
	return 0;
}

这里放一下后来写的通用快速幂函数（用的时候要当心数据超出长整型就变0了）

long long pow(long long base,long long index)
{
	long long ans=1;
	while(index)
	{
		if(index&1)
			ans=ans*base;
		base=base*base;
		index=index>>1;
	}
	return ans;
}