MITE4 论如何科学求幂

快就完事了，acm最速传说！

从初中开始，我们就知道a^b这个叫做a的b次幂（但是我直到高中才接受这个说法

对于求幂，大家肯定就想到pow（a，b）这个math.h里有的函数,但是这个函数是浮点型的，会有误差。

详见[C/C++]C语言中math.h和cmath的pow()精度问题

但是问题不大，自己写一个整数的pow有多难呢是吧

ll pow(ll a,int b,ll p )
{
	ll ret=1;
	for(int i=0;i<b;i++)
		ret=(ret*a)%p;
	return ret;
}

看吧，很简单吧，就几行，甚至还带上了取模。但是问题来了，虽然这个能完成幂操作，但是速度是不是太慢了

每个数字都要o（n），如果b大一点岂不是要炸了对吧。这样是上不了秋名山的！

于是有我们就有了高改ae86幂，也叫快速幂，只要用了这个，就可以惯性过弯了

ll qpow(ll a,ll b,ll p)
{
	ll ret=1;
	while(b)
	{
		if(b&1)ret=(ret*a)%p;
		a=(a*a)%p;
		b>>=1;
	}
}

原理也很简单比如 

2^18=(2^2)^9=4^9=4^8*4=(4^2)^4*4=16^4*4=(16^2)^2*4=256^2*4=65536*4=262144

复杂度为o（logn）

好了，这个方法其实以及足够快了。但是要注意一些细节，就是如果p很大，a*ret和a*a可能会溢出。

于是又有了快速乘法····（其实一点都不快），它可以避免溢出（只要p不要大的离谱，比如说1e18，但对那么大的数取模有个屁意义，不过1e18其实也是可以避免溢出的···所以无视括号里这段废话吧），有时也可以判断溢出。

ll qmul(ll a,ll b,ll p)
{
	ll ret = 0;
	while(b)
	{
		if(b&1)ret=(ret+a)%p;
		a=(a<<1)%p;
		b>>=1;
	}
	return ret;
}

然后，ll范围里的幂运算我们基本都能解决了，但是，如果有ll范围外的幂运算呢

比如 a<1e6 b<1e1000 p<1e5这种范围的题呢？可能有人会想高精度除法什么的，那个也不是不行，但又更好的方法

那就是欧拉降幂

讲到欧拉降幂肯定要说到欧拉函数和欧拉定理了，但是那个我迟点开一个博客专门说吧，这里就直接给出公式和代码就算了。

a^b mod p = a^(b mod phi(p)+phi(p)) mod p   phi()是欧拉函数

ll phi(ll x)
{
	ll rex = a,a = x;
	for(ll i=2;i*i<=x;i++)
		if(a % i == 0)
		{
			res = res/i*(i-1);
			while(a%i==0)
				a/=i;
		}
		if(a>1)
			res = res/a*(a-1);
		return res;
}

int main()
{
	ll a,p;
	char b[1000];
	cin>>a>>b>>p;
	p = phi(p);
	ll ans=0;
	for(int i=0;i<len;i++)
		ans=(ans<<3+ans<<1+b[i]-'0')%p;
	ans+=p;
	cout<<qpow(a,ans,p)<<endl;
}

复杂度大概sqrt(p)+log p吧 我也没算过，代码的话都没编译过，都是直接打的。。。。有错的话直接留言吧（极其不负责的我）