a*x ≡ 1 (mod p) 满足a乘以x对p取模等于1 ,此时 ,称 x为a对p的逆元。
只有a与p互质才有逆元 ,互质 即gcd(a,p)=1
下面介绍两种求法:
1.费马小定理(定义来自维基百科):
如果a不是p的倍数,这个定理也可以写成
题目: hdu1576 A/B
//费马小定理求逆元,p需为质数
#include<cstdio>
#include<iostream>
using namespace std;
typedef long long ll;
const int MOD=9973;
ll PowerMod(ll a,ll b,ll p){ //(a^b)%p
ll ans=1;
a=a%p;
while(b>0){
if(b&1) ans=(ans*a)%p;
b>>=1;
a=(a*a)%p;
}
return ans;
}
int main()
{
ll n,b;
int t;
scanf("%d",&t);
while(t--){
scanf("%d%d",&n,&b); //求(n/b)%MOD
printf("%lld\n",(n*PowerMod(b,MOD-2,MOD))%MOD);
}
return 0;
}
2.扩展欧几里得(定义来自维基百科):
已知整数a、b,扩展欧几里得算法可以在求得a、b的最大公约数的同时,能找到整数x、y(其中一个很可能是负数),使它们满足贝祖等式
。
好了,在明白上面的定理后我们开始分析乘法逆元:ax≡1 (mod p) 这个等式用中文描述就是 a乘一个数x并模p等于1,即 a%p*x%p=res,res%p=1;看上去就是同余定理的一个简单等式- -。那么问题来了。
为什么可以用费马小定理来求逆元呢?
由费马小定理 ap-1≡1 , 变形得 a*ap-2≡1(mod p),答案已经很明显了:若a,p互质,因为a*ap-2≡1(mod p)且a*x≡1(mod p),则x=ap-2(mod p),用快速幂可快速求之。
为什么可以用扩展欧几里得求得逆元?
我们都知道模就是余数,比如12%5=12-5*2=2,18%4=18-4*4=2。(/是程序运算中的除)
那么ax≡1 (mod p)即ax-yp=1.把y写成+的形式就是ax+py=1,为方便理解下面我们把p写成b就是ax+by=1。就表示x是a的模b乘法逆元,y是b的模a乘法逆元。然后就可以用扩展欧几里得求了。
知道逆元怎么算之后,那么乘法逆元有什么用呢?
做题时如果结果过大一般都会让你模一个数,确保结果不是很大,而这个数一般是1e9+7,而且这个数又是个素数,加减乘与模运算的顺序交换不会影响结果,但是除法不行。有的题目要求结果mod一个大质数,如果原本的结果中有除法,比如除以a,那就可以乘以a的逆元替代。(除一个数等于乘它的倒数,虽然这里的逆元不完全是倒数,但可以这么理解,毕竟乘法逆元就是倒数的扩展)。
题目:hdu1576 A/B
//扩展欧几里德求逆元 #include<cstdio> #include<iostream> using namespace std; typedef long long ll; const int MOD=9973; void exgcd(ll a,ll b,ll& d,ll& x,ll& y){ if(!b){ d=a; x=1; y=0; } else{ exgcd(b,a%b,d,y,x); y-=x*(a/b); } } ll inv(ll a,ll p){ //a模p的逆元,扩展欧几里德求逆元,a,p需互质 ll d,x,y; exgcd(a,p,d,x,y); return d==1?(x+p)%p:-1; //不互质的话,返回-1 } int main() { ll n,b; int t; scanf("%d",&t); while(t--){ scanf("%d%d",&n,&b); //求(n/b)%MOD printf("%lld\n",(n*inv(b,MOD))%MOD); } return 0; }