快速幂

快速幂

食用范围

在 $O(logk)$ 的时间内求出 $a^k\mathrm{m}\mathrm{o}\mathrm{d}p$ 的结果，其中 $a,p,k\le 10^9$

核心思想

反复平方法，预处理出

$a^{2^0}\mathrm{m}\mathrm{o}\mathrm{d}p$

$a^{2^1}\mathrm{m}\mathrm{o}\mathrm{d}p$

$a^{2^2}\mathrm{m}\mathrm{o}\mathrm{d}p....$

$a^{2^{logk}}\mathrm{m}\mathrm{o}\mathrm{d}p$

的结果（$logk$ 个）然后组合出 $a^k\to$ 把它拆成若干个前面乘积的形式 $\to$ 使用二进制表示拆 $k$

Code

返回 $a^k\mathrm{m}\mathrm{o}\mathrm{d}p$ 的值

ll qmi(ll a, ll k, ll p) {
	ll res = 1;
	while(k) {
		if(k & 1) res = res * a % p;
		k >>= 1;
		a = a * a % p;
	}
	return res;
}

快速幂求逆元

希望将除法变成乘法

如果 $a÷b\equiv a\times x(modm)$

那么如果求出 $x$ ，$a÷b$ 就可以转化为 $a\times b^{-1}$ 的形式，避免了除法的复杂运算，并且把除法变成乘法还保证了结果是一个整数。

乘法逆元的定义

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

通俗的解释：给定一个 $b$ 找到一个 $x$ 使得 $b\times x\equiv 1(modm)$

于是，根据费马小定理，有，$b^{p-1}\equiv 1(modp)\to b\times b^{p-2}\equiv 1(modp)\to b^{p-2}$ 就是 $b$ 的逆元

费马小定理：如果 $p$ 是一个质数，而整数 $a$ 不是 $p$ 的倍数，则有，则有 $a^{p-1}\equiv 1(modp)$

无解的情况：$b$ 是 $p$ 的倍数，此时，模数一定为 0，满足上述情况时，均可由费马小定理构造出一组解。

Code

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

注意：请返回在 $0\sim p-1$ 之间的逆元。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 7;
const int N = 1e6 + 7, M = N * 2;
const int inf = 0x3f3f3f3f;
const long long mod = 1e9 + 7;

ll qmi(ll a, ll k, ll p) {
	ll res = 1;
	while(k) {
		if(k & 1) res = res * a % p;
		k >>= 1;
		a = a * a % p;
	}
	return res;
}

int main() {
	int n;
	ll a, k, p;
	scanf("%d", &n);
	while(n--) {
		scanf("%lld%lld", &a, &p);

		ll res = qmi(a, p - 2, p);
		if(a % p) printf("%lld\n", res );//特判不存在的情况，注意p=2时是0次方，所以不能简单的用结果是否为0来判断 
		else printf("impossible\n");
	}
    
	return 0;
}