【刷题】数学知识——快速幂：快速幂求逆元

逆元定义

若$a$能被$b$整除，且$b$和$m$互质，则一定存在一个整数$x$，使得
$a/b\equiv a\ast x(modm)$，则$x$为$b$在模$m$下的乘法逆元，记作$b^{-1}(modm)$

逆元的三条性质

• $b\ast b^{-1}\equiv 1(modm)$

• $b$存在乘法逆元的充要条件是$b$和$m$互质

• 当$m$为质数时，$b^{m-2}$是$b$的乘法逆元

---

性质1证明

因为$a$是$b$的整数倍，所以$a=kb$
$a/b\equiv a\ast b^{-1}(modm)$
$a/b\ast b\equiv a\ast b^{-1}\ast b(modm)$
$kb\equiv kb\ast b^{-1}\ast b(modm)$
又因为b和m互质，可消去b
$k\equiv k\ast b^{-1}\ast b(modm)$
$k\ast (1-b^{-1}\ast b)\equiv 0(modm)$
因为$a$是任意数，所以$k$也是任意数，对于任意$k$上式都要满足，所以
$b^{-1}\ast b\equiv 1(modm)$

性质2证明

当b存在乘法逆元时，根据性质1有$b^{-1}b\equiv 1(modm)$
移项得$1-b{b}^{-1}\equiv 0(modm)$，即$1-b{b}^{-1}$可被$m$整除
等价于$mx=1-b{b}^{-1}$
令$a=m$，$y=b^{-1}$，则有
$ax=1-by$
移项得$ax+by=1$
根据裴蜀定理推论：$a$和$b$互质的充要条件是方程$ax+by=1$有整数解
所以a(也就是m)和b互质。

当a和b互质时，同样根据裴蜀定理推论，$ax=1-by$有整数解。
令$a=m$，又因为$x$时整数，则$1-by$能被$m$整除
即$1-by\equiv 0(modm)$，$by\equiv 1(modm)$
因此这个整数$y$就是$b$的逆元。

性质3证明

由费马小定理，我们知道当p是质数时，有$a^{p-1}\equiv 1(modp)$。
也就是$a\ast a^{p-2}\equiv 1(modp)$
根据定义知$a^{p-2}$是$a$的逆元

---

快速幂求逆元

例如$a=2,p=3$，则$2^{3-2}=2$，$2\ast 2\equiv 1(mod3)$，2在mod 3下的逆元是2。
当$a$和$p$互质时，存在乘法逆元，又因为$p$是质数，所以只要$a$是$p$的倍数就输出impossible。

#include <iostream>
using namespace std;

int n, a, p;

int quick_pow(int a, int b, int p) {
    int ans = 1;
    while(b) {
        if (b & 1) ans = (long long) ans * a % p;
        b >>= 1;
        a = (long long) a * a % p;
    }
    return ans;
}

int main() {
    scanf("%d", &n);
    while(n --) {
        scanf("%d%d", &a, &p);
        if (!(a % p)) printf("impossible\n");
        else printf("%d\n", quick_pow(a, p - 2, p));
    }
    return 0;
}