乘法逆元

乘法逆元

​ 在组合计数中，由于答案过大，需要取模。比如计算$C_n^m$.

​ 为了求出$\frac{a}{b}\equiv x(modm)$中的$x$，就需要求出$b$的逆元$b^{-1}$，因为$b{b}^{-1}\equiv 1$，所以$\frac{a}{b}\equiv a{b}^{-1}$。

1-费马小定理

​ （的推论）$a^p\equiv a(modp)$，当$p$为质数且$(a,p)=1$。

​ 这样$a^{p-2}=a^{-1}$，用快速幂求就好。

ll Reverse(ll op){
    ll poi = mod - 2, base = op, res = 1;
    while(poi){
        if(poi&1)
            res = (res*base) % mod; //在mod>INT_MAX时不适用。
        base = (base*base) % mod;
        poi >>= 1;
    }
    return res;
}

2-扩展欧几里得算法

​ $ax+by=(a,b)$一定有解。

​ 又设 $b{x}^{\prime }+(a\%b)y^{\prime }=(a\%b,b)=(a,b)$，

​ 则$ax+by=bx‘+(a\%b)y’=b{x}^{\prime }+(a-b\ast \lfloor a/b\rfloor )y’$

​ $ax+by=a{y}^{\prime }+b(x^{\prime }-\lfloor a/b\rfloor y^{\prime })$

​ 因此$x=y^{\prime }$，$y=x^{\prime }-\lfloor a/b\rfloor y^{\prime }$。

​ 只要求解$b{x}^{\prime }+(a\%b)y^{\prime }=(a,b)$，将$x^{\prime },y^{\prime }$带入即可求出$x,y$。

​ 这样，我们可以像欧几里得算法一样不断减小系数的规模。只要递归地求解更小规模的方程，最后回溯时将解依次带回就好。递归终点是当$b_i=0$时，这时$(a_i,b_i)=a_i=(a,b)$，$a_i{x}_i+b_i{y}_i=(a_i,b_i)$一组解是$x_i=1,y_i=0$。

​ 当$(a,b)=1$时，$ax+by=1$可用来求解$b$在模$a$意义下的逆元$y$，或者$a$在模$b$意义下的逆元$x$.

void exEuclid(ll a, ll b, ll &d, ll &x, ll &y){
    if(!b)
    	d=a, x=1, y=0;		//d就是gcd
    else{
        exEuclid(b, a % b, d, y, x);
        y -= x * (a / b);	//x=y', y=x' - y'*(a/b)
    }
}

​ 复杂度与辗转相除法相同$O(logab)$，因为每一次都使较大的数减小一半以上。

3-线性递推法

​ 需要$mod$是质数。

​ 设$mod=k\ast i+r\equiv 0$，

​ 两边同乘$i^{-1}{r}^{-1}$得

​ $k\ast r^{-1}+i^{-1}\equiv 0$即

​ $i^{-1}\equiv mod-k{r}^{-1}$。

​ 其中$k=\lfloor mod/i\rfloor$，$r=mod\%i$。

rev[1]=1;
for(int i=2; i<mod; i++)	
    rev[i]=(mod - mod/i)*rev[mod%i]%mod;

例题：[SCOI2010]生成字符串

许多需要求组合数的题目中都需要求逆元。本题求解的思路与求卡特兰数的推导过程非常相似。