关于逆元

逆元详细解释：

[对于缩系(https://baike.baidu.com/item/缩系/4763557?fr=aladdin)中的数字]逆元一般用扩展欧几里得算法来求得，

扩展欧几里得：
这个方程可以转化为ax-my=1 
然后套用求二元一次方程的方法，用扩展欧几里得算法求得一组x0,y0和gcd 
检查gcd是否为1 
gcd不为1则说明逆元不存在 

给定模数m，求a的逆相当于求解ax=1(mod m)

若为1，则调整x0到0~m-1的范围中即可

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

[cpp]  view plain  copy

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. }  

费马小定理：

如果为素数，那么还可以根据费马小定理得到逆元为。

 

一般情况(条件b|a已知)：

 ans=a/bmodm=amod(mb)/b

PS:费马小定理和扩展欧几里得算法求逆元是有局限性的，它们都会要求与互素，如果a与m不互素，那就没有逆元，这个时候需要a mod (bm)/b来搞（此时就不是逆元的概念了）。但是当a与m互素的时候，bm可能会很大，不适合套这个一般公式，所以大部分时候还是用逆元来搞

相关链接：

http://blog.csdn.net/guhaiteng/article/details/52123385

http://blog.csdn.net/acdreamers/article/details/8220787

打表法求逆元：

其实有些题需要用到1-p模p的所有逆元，这里p为奇质数。那么如果用快速幂求时间复杂度为O(p log(p))， 
如果对于一个1000000级别的素数，这样做的时间复杂度是很高了。实际上有的算法复杂度可以压缩到O(p)，有一个递推式如下

inv[i]=(M-M/i)*inv[M%i]%M (其中M为模数，要求为奇质数)

相关链接：

http://blog.csdn.net/Frods/article/details/53868810