除法逆元

除法逆元可求（a/b）%mod

也就等于（a%mod）*（inv（b））%mod              （inv（b）就是b的逆元）

inv（b） =  （b^(mod-2)）%mod     (可用快速幂求解)

证明可查看一下博客，写得非常不错，还有更多的逆元求解方法

博客传送门：点击打开链接

-------------华丽的分割线-------------

2018.5.3

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

//费马小定理求逆元
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;
    }
}