首先对逆元不太了解的 强推这篇 可爱风 的博客 (反正我jio得太TM可爱了)https://www.cnblogs.com/linyujun/p/5194184.html
注释: inv 为逆元
方法一:费马小定理
a^(p-1) ≡1 (mod p)
两边同除以a
a^(p-2) ≡1/a (mod p)
什么(,,• ₃ •,,),这可是数论,还敢写1/a
应该写a^(p-2) ≡ inv(a) (mod p)
所以inv(a) = a^(p-2) (mod p)
这个用快速幂求一下,复杂度O(logn)(ง •̀_•́)ง
LL pow_mod(LL a, LL b, LL p){//a的b次方求余p
LL ret = 1;
while(b){
if(b & 1) ret = (ret * a) % p;
a = (a * a) % p;
b >>= 1;
}
return ret;
}
LL Fermat(LL a, LL p){//费马求a关于p的逆元
return pow_mod(a, p-2, p);
}
方法二:扩展欧几里德算法
a*x + b*y = 1
如果ab互质,有解
这个解的x就是a关于b的逆元
y就是b关于a的逆元
为什么呢?