乘法逆元

在取模运算的除法取模那里，我们说取模除法并不能直接对除数取模。在离散数学中有群的概念，同余乘法构成一个群，我们在求除法取模时要把除以一个数当成乘以一个数的倒数，这个倒数在群论中被叫做逆元，因此我们要求的就是除数$x$的乘法逆元$inv(x)$。

乘法逆元

已知$a/b(modp)$，$p$为质数

若$ax\equiv 1(modp)$ ，且$a$与$p$互质，那么我们就能定义: $x$为$a$模$p$意义下的逆元，记为 $a^{-}1$或$inv(a)$。

下面简称乘法逆元为逆元

扩展欧几里得求逆元

此方法的优点是$p$不是质数但满足$a$和$p$互质时也可以使用

由$ax\equiv 1(modp)$，我们可以转化为求解$ax+py=1$，再根据扩展欧几里得算法即可求解$x$

由于最后的求得的$x$可能是负数，因此我们再进行处理：$(x\%p+p)\%p$

typedef long long ll;

ll exgcd(ll a,ll b,ll &x,ll &y){
    if(!b){
        x=1;y=0;
        return a;
    }
    ll gcd=exgcd(b,a%b,y,x);
    y-=(a/b)*x;
    return gcd;
}

ll inv(ll a,ll p){
	ll x,y;
	exgcd(a,p,x,y);
    return (x%p+p)%p;
}

快速幂取模求逆元

费马小定理

若$p$为素数，$a$为正整数，且 $a,p$互质。 则有 $a^{p-1}\equiv 1(modp)$

由于证明过程较为繁琐，就不学证明了，直接运用：

由$ax\equiv 1(modp)$，再根据费马小定理得$ax\equiv a^{p-1}(modp)$

所以$x\equiv a^{p-2}(modp)$

这样我们直接使用快速幂取模就可以求出逆元:

ll quick_mod(ll x,ll n,ll p){
	ll ans=1;
    while(n){
        if(n&1) ans=ans*x%p;
        x=x*x%p;    
        n>>=1;
    }
	return ans;
}

ll inv(ll a){
    return quick_mod(a,p-2,p);
}

线性求逆元

1.求$1-n$在模$p$下的所有乘法逆元（$n<p$且$p$为质数）

如果一个个求，显然时间复杂度有点高，因此我们试着找下递推的规律

令$x=p/i$ , $y=p\%i$, 则有：$x\ast i+y=p$

然后 $x\ast i+y\equiv 0(modp)$

再有 $y\equiv -x\ast i(modp)$

同余式两边同除以 $y\ast i$ 得 $y/(y\ast i)\equiv -x\ast i/(y\ast i)(modp)$

将除法转化为乘法 $inv[i]\equiv -x\ast inv[y](modp)$

最后将$x=p/i,y=p\%i=p-p/i$带入就可以得到:

递推公式$inv[i]\equiv (p-p/i)\ast inv[p\%i](modp)$

const int N=?;
int inv[N]
void solve(int n,int p){
	inv[1]=1;
	for(int i=2;i<=n;i++)
    	inv[i]=(p-p/i)*inv[p%i]%p;
}

2.求阶乘$1!-n!$的逆元

设$fac$数组保存的是$i$的阶乘，类似上面的推导，我们有：

$fac[i]\ast inv[i]\equiv 1(modp)$

由$fac[i+1]\ast inv[i+1]\equiv 1(modp)$ 得 $fac[i]\ast (i+1)\ast inv[i+1]\equiv 1(modp)$

由同余的除法可得 : $inv[i]\equiv inv[i+1]\ast (i+1)$

故我们先由快速幂取模求出$inv[n!]$再往下递推求出所有的逆元

ll fac[N];
int inv[maxn];
ll p;

ll quick_mod(ll x,ll n,ll p){
	ll ans=1;
    while(n){
        if(n&1) ans=ans*x%p;
        x=x*x%p;    
        n>>=1;
    }
	return ans;
}

void solve(int n){ //求1!-n!的所有逆元
	fact[0]=1;
	for (int i=1;i<=n; i++) {
		fact[i]=fact[i-1]*i%p;
	}
	inv[n]=quick_mod(fact[n],p-2,p);
	for (int i=n-1;i>=0;i--) {
		inv[i]=inv[i+1]*(i+1)%mod;
	}
}