【算法讲2.5：拓欧求逆元】算法讲2的补充

前置知识

• 拓展欧几里得
• 逆元：在取模意义下 $ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}MOD)$，$a$ 即为在模 $MOD$ 意义下 $x$ 的逆元

简单推导

• 由于我们有拓展欧几里得算法，能方便求出该方程 $ax+by=c$ 的通解。
  我们想求 $x$ 的逆元，本质就是求出 $ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}MOD)$，这时候的 $a$ 为多少。
  如果方程在取模 $MOD$ 意义下，方程一边加上 $b$ 倍的 $MOD$，方程等号不会变，即变为：
  $a\times x+b\times MOD\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}MOD)$
• 由于在取模意义下， 若存在逆元，则 $c=\gcd (x,MOD)=1\text{gcd(x,MOD)=1}\gcd (x,MOD)=1$。
  注意0没有逆元，这一点也是满足上式的。
  式子变成 $ax+b\times MOD\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}MOD)$
• 我们使用拓欧算法求出的系数 $a$ 就直接是 $x$ 的逆元了。注意该系数可能为负数，需要变成正数即可。

核心算法

ll ex_gcd(ll a,ll b,ll& x,ll& y)
{
    if(b==0)
    {
        x=1;y=0;
        return a;
    }
    ll ans=ex_gcd(b,a%b,x,y);
    ll tmp=x;
    x=y;
    y=tmp-a/b*y;
    return ans;
}
ll inv(ll a,ll mod)//存在逆元条件：gcd(a,mod)=1
{
    ll x,y;
    ll g=ex_gcd(a,mod,x,y);
    if(g!=1)return -1;
    return (x % mod + mod) % mod;
}