快速幂算法

3^10=333333333*3

//尽量想办法把指数变小来，这里的指数为10

3^10=(33)(33)(33)(33)(3*3)

3^10=(3*3)5

3¹⁰⁼⁹5

//此时指数由10缩减一半变成了5，而底数变成了原来的平方，求3^{10原本需要执行10次循环操作，求9}5却只需要执行5次循环操作，但是3^10却等于95,我们用一次（底数做平方操作）的操作减少了原本一半的循环量，特别是在幂特别大的时候效果非常好，例如2^10000=45000,底数只是做了一个小小的平方操作，而指数就从10000变成了5000，减少了5000次的循环操作。

//现在我们的问题是如何把指数5变成原来的一半，5是一个奇数，5的一半是2.5，但是我们知道，指数不能为小数，因此我们不能这么简单粗暴的直接执行5/2，然而，这里还有另一种方法能表示9^5

9^5=（94）*（9^1）

//此时我们抽出了一个底数的一次方，这里即为9^1，这个91我们先单独移出来,剩下的9^4又能够在执行“缩指数”操作了，把指数缩小一半，底数执行平方操作

9^5=（812）*(9^1)

//把指数缩小一半，底数执行平方操作

9^5=（65611）*(9^1)

//此时，我们发现指数又变成了一个奇数1，按照上面对指数为奇数的操作方法，应该抽出了一个底数的一次方，这里即为6561^{1，这个6561}1我们先单独移出来，但是此时指数却变成了0，也就意味着我们无法再进行“缩指数”操作了。

9^5=（65610）(9^1)*(65611)=1(9^1)*(65611)=(9^1)*(65611)=9*6561=59049

我们能够发现，最后的结果是9*6561，而9是怎么产生的？是不是当指数为奇数5时，此时底数为9。那6561又是怎么产生的呢？是不是当指数为奇数1时，此时的底数为6561。所以我们能发现一个规律：最后求出的幂结果实际上就是在变化过程中所有当指数为奇数时底数的乘积。

#include<bits/stdc++.h>
using namespace std;
long long fastpower(long long base,long long power)
{
	long long result =1;
	while(power>0)
	{
		if(power%2==0)
		{
			power =power/2;
			base =base*base %1000;
		}
		else{
			power =power-1;
			result=result *base%1000;
			power=power/2;
			base=base*base%1000;
		}
	}
	return result;
}
int main(){
	long long base,power;
	long long  result;
		cin>>base>>power;
	long long t=fastpower(base,power);


	cout<<t;
	return 0;
}
输入   2               100
输出 326