求逆元

一、O(n)求逆元

$inv[i] = ( MOD - \frac{MOD}{i} ) * inv[MOD\%i] \% MOD$

证明：

设$t=\frac{MOD}{i} , k=MOD % i$
则有 $t * i + k\equiv0 \pmod {MOD}$
有 $-t * i\equiv k \pmod { MOD}$
两边同时除以$i*k$得到
$-t * inv[k]\equiv inv[i]\pmod { MOD}$
即$inv[i]\equiv -\frac{MOD}{i}* inv[MOD\%i]\pmod { MOD}$
即$inv[i]\equiv ( MOD - \frac{MOD}{i} * inv[MOD\%i]\pmod { MOD}$
证毕

适用于MOD是质数的情况，能够O(n)时间求出1~n对模MOD的逆

代码实现：

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

二、费马小定理

在模为素数p的情况下，有费马小定理
$a^{p-1}\equiv1\pmod p$
∴$a^{p-2}\equiv a^{-1}\pmod p$
也就是说a的逆元为$a^{p-2}$
而在模不为素数p的情况下，有欧拉定理
$a^{\phi m}\equiv1\pmod m$
∴$a^{-1}\equiv a^{\phi m-1}\pmod m$
因此逆元x便可以套用快速幂求得了$x=a^{\phi m-1}$
但是似乎还有个问题？如何判断a是否有逆元呢?
检验逆元的性质，看求出的幂值x与a相乘是否为1即可
PS:这种算法复杂度为$O（log_2N）$在几次测试中，常数似乎较上种方法大所以当p比较大的时候需要用快速幂求解。

代码实现：

typedef  long long ll;  
ll pow_mod(ll x, ll n)
{  
    ll res=1;  
    while(n>0)
    {  
        if(n&1)res=res*x%mod;  
        x=x*x%mod;  
        n>>=1;  
    }  
    return res;  
}  
pow_mod(x,mod-2);