求解乘法逆元

费马小定理

适用条件：$mod$是素数
整数 $a$，质数 $p$
如果$a$${}^p$ $-a$ 是 $p$ 的倍数，$a^p\equiv a(mod$ $p)$
如果 $a$ 不是 $p$ 的倍数，$a$${}^{p-1}$ $\equiv 1(mod$ $p)$
所以 $a$${}^{p-2}$ * a $\equiv 1(mod$ $p)$
所以$a$${}^{p-2}$是 $a$ 在模$p$下的逆元

ll M(ll a) {return a%mod;}//取模
ll Pow(ll a,ll b = mod-2){ll c; while(b){if(b&1) c=M(c*a); b>>=1; a=M(a*a);} return c;}//快速幂

扩展欧几里得

适用条件：存在逆元，个数不多，$mod$很大
$ax\equiv 1(modp)$即$ax-yp=1$.
把$y$写成$+$的形式就是$ax+py=1$
$x$ 是 $a$ 的模 $p$ 乘法逆元，$y$ 是 $b$ 的模 $a$ 乘法逆元。

void exgcd(ll a,ll b,ll &d,ll &x,ll &y) {//扩展欧几里得算法 
    if(!b) { d = a; x = 1; y = 0; return;}
    exgcd(b, a%b, d, y, x); 
    y -= x*(a/b);
}

ll inv(ll a, ll p) {//求a在p下的逆元，不存在逆元返回-1 
    ll d, x, y;
    exgcd(a, p, d, x, y);
    return d==1 ? (x+p)%p : -1;
}

递归法

ll inv(int i) {return i==1 ? 1 : (mod - mod/i) * inv(mod%i) % mod;}