中国剩余定理及扩展欧几里得算法

中国剩余定理

已知m1、m2、m3是两两互质的正整数,求最小正整数x,使它被m1、m2、m3除所得余数分别为C1、C2、C3 。
孙子定理的思想便是先分别找出被其中数mi除余1而被另二数整除的数Mi(i=1,2,3),则所求的数之一便是   C1M1+C2M2+C3M3;
若欲求的是最小的符合要求的数,则将上面的得数减去m1、m2、m3的整数倍(0,1,2,…)即可. 
在古算题中,m1=3,m2=5,m3=7;C1=2,C2=3,C3=2;M1=70,M2=21,M3=15.   
其中 : M1=70=3×23+1=5×7×2;   
M2=21=5×4+1=3×7×1;   
M3=15=7×2+1=3×5×1;
而 C1M1+C2M2+C3M3=2×70+3×21+2×15=233   
∵233>2×3×5×7=2×105,故所求最小数为        233-2×105=23   

孙子定理可以推广到对任意n个数mi的情形,n≥2,n∈N。
附:中国剩余定理的五种解法:http://blog.sina.com.cn/s/blog_a6f9a3b60101favb.html

两个定理

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
快速模幂运算: 在实现快速模幂运算时,需要注意到一个简单的公式: a ^ b % m = (a % m) ^ b % m 这个公式可以极大地减少计算量,加快计算速度。具体实现时,可以通过不断将指数b除以2,并将底数a平方,直到指数b为0,最后将结果对模数m取余即可得到快速幂结果。代码实现如下: ```python def fast_power(base, exponent, mod): result = 1 base = base % mod while exponent > 0: if exponent % 2 == 1: result = (result * base) % mod exponent = exponent // 2 base = (base * base) % mod return result ``` 扩展欧几里得: 扩展欧几里得算法是用于求解两个整数a和b的最大公约数gcd以及一组满足 ax+by=gcd 的整数解x和y。在实现该算法时,需要注意到以下性质: - gcd(a, b) = gcd(b, a % b) - 若ax+by=gcd(a,b),则bx'+(a%b)y'=gcd(b,a%b),其中 x'=y, y'=x-(a//b)y 从而可以递归调用该算法,直到b=0为止。最后返回gcd以及最后一次迭代中求得的x和y即可。代码实现如下: ```python def extended_euclidean_algorithm(a, b): if b == 0: return a, 1, 0 else: gcd, x, y = extended_euclidean_algorithm(b, a % b) return gcd, y, x - (a // b) * y ``` 中国剩余定理中国剩余定理是用于求解一组同余方程组的解。在实现该算法时,需要注意到以下性质: - 若m1,m2互质,则同余方程组x≡a1(mod m1),x≡a2(mod m2)的解为x=(a1m2M2+a2m1M1)(mod M),其中M=m1m2M1=M/m1M2=M/m2,且M1x1≡1(mod m1),M2x2≡1(mod m2) - 若m1,m2不互质,则同余方程组x≡a1(mod m1),x≡a2(mod m2)无解 可以逐个求解每个同余方程,最后通过合并得到一个更小的同余方程组来求解。具体实现时,可以使用扩展欧几里得算法求解M1M2的乘法逆元,从而计算出x1和x2,再代入上述公式即可得到最终解。代码实现如下: ```python def chinese_remainder_theorem(a, m): M = 1 for mi in m: M *= mi result = 0 for i in range(len(a)): Mi = M // m[i] gcd, x, y = extended_euclidean_algorithm(Mi, m[i]) result += a[i] * Mi * x result %= M return result ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值