除法逆元

转自：https://blog.csdn.net/baodream/article/details/77822634

总结一下几种求逆元的方法

//费马小定理求逆元
ll quick_mod(ll a,ll b,ll c) //快速幂计算(a^b)%c
{
    ll ans = 1;
    while(b)
    {
        if(b&1)   //相当于b%2==1
            ans = (ans*a)%c;
        a = (a*a)%c;
        b>>=1;    //相当于b/=2
    }
    return ans;
}
 
ll inv(ll b,ll c)   //计算b的逆元
{
    return quick_mod(b,c-2,c);
}
 
ll div(ll a,ll b,ll c)  //计算(a/b)%c
{
    return ((a%c)*(inv(b,c)%c))%c;
}
 
//扩展GCD求逆元
void ex_gcd(ll a, ll b, ll &d, ll &x, ll &y){
    if (!b) {d = a, x = 1, y = 0;}
    else{
        ex_gcd(b, a % b, d, y, x);
        y -= x * (a / b);
    }
}
 
ll inv(ll a, ll p){//如果不存在，返回-1
    ll d, x, y;
    ex_gcd(a, p, d, x, y);
    return d == 1 ? (x % p + p) % p : -1;
}
 
//递归求逆元
//当p是个质数的时候有inv(a) = (p - p / a) * inv(p % a) % p,且1的逆元就是1
ll inv(ll t, ll p) {
    //求t关于p的逆元，注意:t要小于p，最好传参前先把t%p一下 ,即inv(a%p, p)求a对p的逆元
    return t == 1 ? 1 : (p - p / t) * inv(p % t, p) % p;
}
 
//打表求逆元
ll inv[maxn];
void Prepare_inv(ll n,ll M){
    inv[1]=1;
    for(ll i=2;i<=n;i++){
        inv[i]=(ll)(M-M/i)*inv[M%i]%M;
    }
}