数论文章----关于逆元的求法（欧拉定理，阶乘逆元，费马小定理，模质数p的情况）

乘法逆元

对于缩系中的元素，每个数a均有唯一的与之对应的乘法逆元x，使得ax≡1(mod n)
一个数有逆元的充分必要条件是gcd(a,n)=1，此时逆元唯一存在
逆元的含义：模n意义下，1个数a如果有逆元x，那么除以a相当于乘以x。

 

下面给出求逆元的几种方法：

1.扩展欧几里得

给定模数m，求a的逆相当于求解ax=1(mod m)
这个方程可以转化为ax-my=1
然后套用求二元一次方程的方法，用扩展欧几里得算法求得一组x0,y0和gcd
检查gcd是否为1
gcd不为1则说明逆元不存在
若为1，则调整x0到0~m-1的范围中即可

PS：这种算法效率较高，常数较小，时间复杂度为O(ln n)

 

1. typedef long long ll;

2. void extgcd(ll a,ll b,ll& d,ll& x,ll& y){

3. if(!b){ d=a; x=1; y=0;}

4. else{ extgcd(b,a%b,d,y,x); y-=x*(a/b); }

5. }

6. ll inverse(ll a,ll n){

7. ll d,x,y;

8. extgcd(a,n,d,x,y);

9. return d==1?(x+n)%n:-1;

10. }

 

 

2.费马小定理

在模为素数p的情况下，有费马小定理
a^(p-1)=1（mod p）
那么a^(p-2)=a^-1(mod p)
也就是说a的逆元为a^(p-2)

而在模不为素数p的情况下，有欧拉定理
a^phi(m)=1（mo