乘法逆元

前言: 乘法逆元,一般用于将模意义下的除法转换为乘法

求解普通逆元($ax\equiv 1(modm)$)

拓展欧几里得

首先我们看看逆元的定义:若$ax\equiv 1(modm)$,且a与m互质,则$x$为$a$在$modm$意义下的乘法逆元,记作$a^{-1}$.

根据这个同余式,我们可以得到$ax+km=1$,且a和m是互质的,那么就满足贝祖等式$ax+by=gcd(a,b)$,可以用拓欧求解.

ll exgcd(ll a,ll b,ll &x,ll &y)//ax+by
{
	if(b==0)
	{
		x=1,y=0;
		return a;
	}
	int ret=exgcd(b,a%b,x,y),t=x;
	x=y,y=t-(a/b)*x;
	return ret;
}

这个板子虽然稍长,但返回值是gcd,用处更广泛

另外,由于gcd求得负数和1最大公约数是-1,exgcd不能用来求解a为负数的情况

快速幂

观察乘法逆元的形式:$ax\equiv 1(modm)$.

考虑到欧拉定理:$a^{\phi (m)}\equiv 1(modm)$.

由于洛谷模版题中给出m为质数,则有:
$$\begin{gathered} ax\equiv a^{m-1}(modm) \\ a\equiv a^{m-2}(modm) \end{gathered}$$
$那么,a在modm意义下的乘法逆元即为a^{m-2}modm$,满足快速幂形式,可使用快速幂求解.

ll fpm(ll x, ll power, ll mod) 
{
    x%=mod;
    ll ans=1;
    for (;power;power>>=1,(x*=x)%=mod)
    	if(power&1)(ans*=x)%=mod;
    return ans;
}

ll inv(ll x)
{
	return fpm(x,mod-2,mod);
}

代码中的(ans*=x)%=mod即为(ans=ans*x)=(ans=ans*x)%mod.