4. 快速幂

1. 求快速幂

给定 n 组 ai,bi,pi，对于每组数据，求出 abiimodpi 的值。

输入格式
第一行包含整数 n。

接下来 n 行，每行包含三个整数 ai,bi,pi。

输出格式
对于每组数据，输出一个结果，表示 abiimodpi 的值。

每个结果占一行。

数据范围
1 ≤ n ≤ 100000,
1 ≤ ai,bi,pi ≤ 2×109

注意的一点就是灵巧用二进制来进行表示
关于二进制的用法：

• b&1就是判断b的二进制表示中第0位上的数是否为1，若为1，b&1=true，反之b&1=false
• b&1也可以用来判断奇数和偶数，b&1=true时为奇数，反之b&1=false时为偶数
• 取整数n的二进制数的第k位数 n>>k&1
  

#include<iostream>
using namespace std;

typedef long long LL;

int qmi(int a, int k, int p) {
	int res = 1;	//存储答案

	while (k) {	//将l转化为二进制表示
		//如果k的二进制在当前位（即0位）为1即乘以当前的a，a耶在循环中进行变换
		if (k & 1) {	//判断k的0位是否为1
			res = (LL)res * a % p;
		}

		k = k >> 1;	//将k的二进制右移一位

		a = (LL)a * a % p;	//将a同步变换，即求a的不断的平方
	}
	return res;
}

int main() {
	int n, a, k, p;
	cin >> n;
	
	while (n--) {
		cin >> a >> k >> p;
		cout << qmi(a, k, p) << endl;
	}
	return 0;
}

2. 快速幂求逆元

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

给定 n 组 ai,pi，其中 pi 是质数，求 ai 模 pi 的乘法逆元，若逆元不存在则输出 impossible。

注意：请返回在 0∼p−1 之间的逆元。

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

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

数据范围
1 ≤ n ≤ 10^5,
1 ≤ ai,pi ≤ 2∗10^9

费马小定理

！费马小定理： 若p是一个质数并且a不是p的倍数，则有a^(p-1) == 1 (mod p)

#include<iostream>
using namespace std;

typedef long long LL;

int qmi(int a, int k, int p) {
	int res = 1;	//存储答案

	while (k) {	//将l转化为二进制表示
		//如果k的二进制在当前位（即0位）为1即乘以当前的a，a耶在循环中进行变换
		if (k & 1) {	//判断k的0位是否为1
			res = (LL)res * a % p;
		}

		k = k >> 1;	//将k的二进制右移一位

		a = (LL)a * a % p;	//将a同步变换，即求a的不断的平方
	}
	return res;
}

int main() {
	int n, a, k, p;
	cin >> n;
	while (n--) {
		cin >> a >> p;
		int res = qmi(a, p - 2, p);
		if (a % p)	cout << res << endl;	//不能用res的结果作为判断条件，因为是看a是否为p的倍数来进行判断是否输出，若用res则是看a的p-2次方是否为p的倍数
		else cout << "impossible" << endl;
	}
	return 0;
}