求a/b mod m

2019-04-13 20:48发布

今天做题遇到需要求a/b mod m的值,这里b能整除a,将方法和大家分享,我是ACM初学者,找了大半天才找到方法,大牛们就见笑啦。
1)求b关于模m的逆元p;
什么叫逆元:
设m为正整数,a为正整数,如果存在a' 使得:
   a * a' = 1(mod m)
  成立,则a叫模m的可逆元,a' 叫a模m的逆元。
设b关于模m的逆元为p,则(b*p)mod m=1. b有逆元的条件是gcd(b,m)==1,显然
如果m为素数肯定有逆元。

求逆元的方法:
1)扩展欧几里得算法,extEuclid(b,m,x,y),b的系数x就是b关于模m的逆元;该
算法就不叙述了
2)如果m为素数,则b的(m-2)次方 mod m的值就是b关于模m的逆元,因为根据费马
小定理,b的(m-1)次方 mod m=1.
2)(a*p)mod m的值即为所求。

这里证明一下该结论的正确性:
因为 (b*p)mod m=1,所以(a/b)mod m = (a/b*b*p)mod m = (a*p)mod m.