C/C++ 快速幂

计算a^b mod p。
时间复杂度：O(logb)
注：十进制数b在二进制下的位数为 int(log₂b)+1。

算法原理：
十进制数b可表示为 b = c_k-12^k-1+ … + c₁2¹+c₀2⁰ （c_i=0或1，k为b的位数）
所以 $a^b=(a^{2^{k-1}}{)}^{c_{k-1}}\cdots (a^{2^1}{)}^{c_1}\cdot (a^{2^0}{)}^{c_0}$
且有 $a^{2^{n+1}}=(a^{2^n}{)}^2$

# include <bits/stdc++.h>
using namespace std;


int main() {
	int a, b, p; scanf("%d %d %d", &a, &b, &p);
	// alg
	int ans = 1 % p;  // 处理 p=1的情况
	while (b) {
		if (b & 1) ans = (long long)ans*a % p;  // ci为 1
		a = (long long)a*a % p;  // 转换类型防止溢出
		b >>= 1;  // i +1
	}
	// output
	printf("%d\n", ans);
	
	return 0;
}