1.扩展欧几里得求乘法逆元
原理就是解线性同余方程;
void exgcd (ll a, ll b, ll &x, ll &y) {
if (b == 0) {
x = 1;
y = 0;
return ;
}
exgcd (b, a % b, x, y);
ll tmp = x;
x = y;
y = tmp - a / b * y;
}
ll x, y;
exgcd (a, b, x, y);//求a在模b意义下的乘法逆元
x = (x % p + p) % p;
printf ("%d\n", x);
2.利用费马小定理
前提是模数必须是质数;因为a^(p-1) mod p=1,所以a模p的乘法逆元就是a^(p-2)
运用快速幂
#include<cstdio>