幂取模

方法一：

利用模等式
（a×b）mod p = （a mod p × b mod p）mod p
2^k mod p = (2^k-1 mod p) × 2 mod p
2^k-1 mod p = (2^k-2 mod p) × 2 mod p
2^k-2 mod p = (2^k-3 mod p) × 2 mod p
2^k-3 mod p = (2^k-4 mod p) × 2 mod p
…
2² mod p = (2¹ mod p) × 2 mod p

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

方法二：

逐次平方法（计算次数较多时）
利用二进制展开指数，同理：
（a×b）mod p = （a mod p × b mod p）mod p
7³²⁷ (mod 853) = 7²⁵⁶ × 7⁶⁴ × 7⁴ × 7² × 7¹ (mod 853)
= 298 × 123 × 695 × 49 × 7 (mod 853)

int P_Mod(int a, int b, int p)
{
	int result = 1;
	while (b) {
		if (b & 1){
			result = result * a % p;
		}
		a = a * a % p;
		b >>= 1;
	}
	return result;
}