数学知识——快速幂

1. 快速幂的用途
在O(logn)的时间复杂度下求出a^kmod p d的结果（1<= a, p, k<=10⁹）

2. 算法原理
反复平方法：

由于在循环中计算，发现每次预处理的后一个值都是前一个值的平方，如果这次循环中k的最后一位是1就乘上对应的a再对p取余，如果最后一位是0就不乘。每次把k的二进制位向前推进一位，并求出下一个a的值。

3. 例题： 给定n组a_i,b_i,p_i，对于每组数据，求出a_ib_i mod p_i的值。

输入格式
第一行包含整数n。
接下来n行，每行包含三个整数a_i,b_i,p_i。

输出格式
对于每组数据，输出一个结果，表示a_ib_i mod p_i的值。
每个结果占一行。

数据范围
1≤n≤100000,
1≤a_i,b_i,pi≤2∗10⁹
输入样例：
2
3 2 5
4 3 9
输出样例：
4
1

#include<iostream>
#include<algorithm>

using namespace std;

typedef long long LL;

int qmi(int a, int k, int p)
{
	int res = 1;
	while(k)
	{
		if(k & 1) res = (LL) res * a % p;	// 如果k的二进制位中的最后一位是1，则把结果乘以相应的幂再取模
		k >>= 1;	// 每次删去k的最后一位
		a = (LL) a * a % p;	 // 预处理出下一个a的值
	}
	return res;
}

int main()
{
	int n; 
	cin >> n;

	while(n--)
	{
		int a, k, p;
		cin >> a >> k >> p;
		printf("%d\n", qmi(a, k, p));
	}
	return 0;
}				

4. 逆元
定义：a∗b≡1(mod p)，那么a，b互为模p意义下的逆元。
推论：如果bc ≡1(mod p)，那么a/b (mod p) = bc *(a/b) (mod p) = a * c (mod p)

如果一个数x对模m存在逆元，则x与m一定互质

适应场景：如果要算x/a%p，就可以改成x*b%p，其中b是a在p下的逆元。

5. 例题2: 快速幂求逆元
给定n组a_i,p_i，其中p_i是质数,求a_i模p_i的乘法逆元，若逆元不存在则输出impossible。
注意：请返回在0∼p−1之间的逆元。

乘法逆元的定义
若整数b，m互质，并且对于任意的整数 a，如果满足b|a，则存在一个整数x，使得a/b≡a∗x(mod m)，则称x为b的模m乘法逆元，记为b⁻¹(mod m)。
b存在乘法逆元的充要条件是b与模数m互质。当模数m为质数时，b^m−2即为b的乘法逆元。

输入格式
第一行包含整数n。
接下来n行，每行包含一个数组a_i,p_i，数据保证p_i是质数。

输出格式
输出共n行，每组数据输出一个结果，每个结果占一行。
若a_i模p_i的乘法逆元存在，则输出一个整数，表示逆元，否则输出impossible。

数据范围
1≤n≤10⁵,
1≤a_i,p_i≤2∗10⁹
输入样例：
3
4 3
8 5
6 3
输出样例：
1
2
impossible

#include<iostream>
#include<algorithm>

using namespace std;

typedef long long LL;

int qmi(int a, int k, int p)
{
	int res = 1;
	while(k)
	{
		res = (LL) res * a % p;
		k >>= 1;
		a = (LL) a * a % p;
	}
	return res;
}

int main()
{
	int n;
	cin >> n;

	while(n--)
	{
		int a, p;
		cin >> a >> p;
		if(a % p == 0) printf("impossible\n");  // 如果a是p的倍数，a%p=0则不存在逆元
		else printf("%d\n", qmi(a, p - 2, p));	
	}
	return 0;
}				

算法原理：求a的逆元即求a*x ≡ 1 (mod p) 中的x，当p为质数时，由费马小定理得 b^p-1 ≡ 1 (mod p) 即 b * b^p-2 ≡ 1 (mod p)
所以b^p-2 （mod p) 即为b的逆元