乘法逆元(费小定理)

乘法逆元(费小定理)

乘法逆元
若$a,b$互质，且满足同余方程$$ax\equiv 1(\mathrm{mod}b)$$
则称$x$为$a$模$b$的乘法逆元，记作$a^{-1}$
费马小定理
若$a$和$p$互质且$p$为质数，则有 $$a^{p-1}\equiv 1(\mathrm{mod}p)$$
证明:
取与$p$互质的序列$A=1,2,\mathrm{3...},p-1$
即证: $$\prod _{i=1}^{p-1}a\times A_i\equiv \prod _{i=1}^{p-1}{A}_i(\mathrm{mod}p)$$
可证得:对于每一个$a\times A_i(\mathrm{mod}p)$的余数是唯一的

假设:$a{A}_i(\mathrm{mod}p)$ 和 $a{A}_j(\mathrm{mod}p)$余数都是$r$
则有
$a{A}_i=xp+r$
$a{A}_j=yp+r$
消去$r$有$a(A_i-A_j)=p(x-y)$
左右一定不等，因为左边不是p的倍数，矛盾

故有$a\times A_i(\mathrm{mod}p)$的最终结果一定是等于序列$A$的，即可得证
根据费马小定理
若有$a,p$，其中$a,p$互质且$p$为质数，则有
$$a^{p-1}\equiv 1(\mathrm{mod}p)$$
$$a\times a^{p-2}\equiv 1(\mathrm{mod}p)$$
由乘法逆元的定义可知$a^{p-2}$就是$a$模$p$意义下的乘法逆元

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
ll ksm(ll a,ll b,ll p){
    ll res = 1;
    while(b){
        if(b & 1)res = res * a % p;
        a = a * a % p;
        b >>= 1;
    }
    return res;
}
int main(){
    ll a,p;cin >> a >> p;
    if(a % p){
        cout << ksm(a, p - 2, p) << endl;
    }
    return 0;
}