【概述】
1.定义:若

,a、b 互质,则称 x 为 a 的逆元,记为

2.同余模公式

3.应用:当题目要求对结果求 m 的模,且当过程需要计算

时,需要对

取模,即

,有时 b 过于大,会出现爆精度的情况,所以需要变除法为乘法。
即:设 c 是 b 的逆元
则:

故:

即:

4.求解逆元方法
1)费马小定理
2)扩展欧几里德算法
3)线性求逆元
【费马小定理求逆元】
1.费马小定理
若 a 为一整数,p 是一质数,且

,那么
2.适用
题目要求模 p 为素数的情况。
3.方法
由

可得

,即得

是 a 的逆元。
4.实现
#define MOD 1000000007
long long quick_pow(long long a,long long b)
{
int res=1,base=a%MOD;
while(b)
{
if(b&1)
res=(base*res)%MOD;
base=(base*base)%MOD;
b>>=1;
}
return res;
}
long long inv(long long a)
{
return quick_pow(a,MOD-2);
}
【扩展欧几里德算法求逆元】
1.适用
题目要求模 m 为素数的情况。
2.方法
由逆元定义可知 ,求 a 模 m 的逆元 x ,即为求解同余方程

将方程转化为

,然后套用解同余方程的方法,用扩展欧几里得算法求得一组

和

然后检查

是否为 1 ,若不为 1 说明逆元不存在,若为 1 ,则调整

到
![[0,m-1]](data/attach/1904/9056vstvjaj0zaqr0aqjphqp144qi9aq.jpg)
的范围中即可。
3.实现
long long Extended_GCD(long long a,long long b,long long &x,long long &y)
{
if(b==0)
{
x=1;
y=0;
return a;
}
long long gcd=Extended_GCD(b,a%b,y,x);
y=y-(a/b)*x;
return gcd;
}
long long inv(long long a)
{
long long x,y;
Extended_GCD(a,MOD,x,y);
x=(x%MOD+MOD)%MOD
return x
}
【线性求逆元】
1.适用
在模质数 p 下,求 1~n 所有逆元,用上面两种方法复杂度差不多都是 O(nlogn),通过递推关系可以线性求出所有逆元。
2.方法
已知,1的逆元一定是1,故有:

①
设

②
将 ② 代入 ① 有:

即:

两边同除

,得:

代入 t、 k, 有:
![inv[i] = (p-p/i) * inv[p:, mod:,i] :, mod:, p](data/attach/1904/43iigdy6fgtc9vkst3d84yn8ezrcf74z.jpg)
综上,易得:
3.实现
int inv[N];
long long Inv(int n)
{
inv[1]=1;
for(int i=2;i<=n;i++)
inv[i]=(p-p/i)*inv[p%i]%p;
}
【说明】
通常情况下,遇到

的问题时,一般通过费马小定理来解决,但是只有当

时,b 的逆元才存在。
对于

的情况,有一通用公式:

其推导如下

对于公式:

,其适用于所有的情况,无需区分互不互素,而费马小定理和扩展欧几里德算法求逆元具有局限性的,它们都会要求 b 与 k 互素,如果两者不互质,那就没有逆元。
当两者互质的时候,b 与 k 可能会很大,不适合套用一般公式,因此大部分时还是使用逆元处理。