【刷题】数学知识——快速幂

用于快速求$a^{b}modp$的结果，时间复杂度是$O(logk)$
其中$1\le a,b,p\le 10^9$

核心思想：

• 预处理$a^{2^b}modp$的结果，$b=0,1,2...$

• 对于任意求的$a^{k}modp$的结果，将其拆成$a^{2^{x_1}+2^{x_2}+...}=a^{2^{x_1}}\ast a^{2^{x_2}}\ast ...$
  例如预处理了$a^1,a^2,a^4,a^8$，则$a^{14}=a^8\ast a^4\ast a^2$，$14$的二进制是$1110$，在1的位置乘上$a^{x_i}$。

• 预处理只需要对$a^{2^b}$平方，就可以得到$a^{2^{b+1}}$，迭代$b-1$次即可。
  如$a^1$平方得到$a^2$，$a^2$平方得到$a^4$，依次类推

---

#include <iostream>
using namespace std;

int n; 

int quick_pow(int a, int b, int p) {
	int ans = 1;
	while (b) {
		if (b & 1) ans = (long long)ans * a % p;  // 转化成二进制，1的位置乘上a对应次幂 
		b >>= 1;
		a = (long long)a * a % p; // 对a^(2^b)进行平方
	}
	return ans; 
}

int main() {
	int a, b, p;
	scanf("%d", &n);
	while (n -- ) {
		scanf("%d%d%d", &a, &b, &p);
		printf("%d\n", quick_pow(a, b, p));
	}
	return 0;
}