lintcode 727. 中国剩余定理 chinese-remainder-theorem

python 中国剩余定理 chinese-remainder-theorem

1.题目

描述
给出两个数组 num[0..k - 1] 和 rem[0..k - 1]. 
在数组num[0..k - 1], 所有的元素都是互质的( gcd 为 1 ). 
我们需要找到满足下列条件的最小正数 x:
 x % num[0]    =  rem[0], 
 x % num[1]    =  rem[1], 
 .......................
 x % num[k-1]  =  rem[k-1] 
样例
样例 1:

输入:[3,4,5][2,3,1]
输出:11
解释:
11 是满足以下条件的最小值:
- 除以 3, 得到余数 2。
- 除以 4, 得到余数 3。
- 除以 5, 得到余数 1。
样例 2:

输入:[3,4,5,11][2,3,1,7]
输出:491
解释:
491 是满足以下条件的最小值:
- 除以 3, 得到余数 2。
- 除以 4, 得到余数 3。
- 除以 5, 得到余数 1。
- 除以 11, 得到余数 7

2.编程思路及代码

solution 1:首先想到的是使用除数中最大的值的倍数+余数,不断增大,直到满足所有要求。可惜Time Limit Exceeded,通不过。

class Solution:
    """
    @param num: the given array
    @param rem: another given array
    @return: The minimum positive number of conditions to meet the conditions
    """
    def remainderTheorem(self, num, rem):
        # write your code here
        mx=max(num)
        ind=num.index(mx)
        mr=rem[ind]
        num.pop(ind)
        rem.pop(ind)
        
        n=len(num)
        k=0
        while True:
            argmin=k*mx+mr
            j=0
            for i in range(n):
                j+=1
                if argmin%num[i]!=rem[i]:
                    break
                if j==n:
                    return argmin
            k+=1

solution 2:开始在网上搜索解决方法,终于找到了孙子定理/中国余数定理。
参考文件[1]和[2]给出了很好的介绍,这里只给出其精华部分。

引自参考文献1
问题1:计算一个整数 [公式] ,使得它满足除以3余2、除以5余3、除以7余2。

如果能够找到三个整数 x 1 , x 2 , x 3 x_1,x_2,x_3 x1,x2,x3 ,使得:

x 1 x_1 x1 除以3余2、除以5余0、除以7余0;

x 2 x_2 x2 除以3余0、除以5余3、除以7余0;

x 3 x_3 x3 除以3余0、除以5余0、除以7余2;

x = x 1 + x 2 + x 3 x=x_1+x_2+x_3 x=x1+x2+x3,易知 x x x 满足除以3余2、除以5余3、除以7余2。

进一步将问题转化为:

x 1 = 2 ∗ y 1 x 2 = 3 ∗ y 2 x 3 = 2 ∗ y 3 x_1=2*y_1\\ x_2=3*y_2\\ x_3=2*y_3 x1=2y1x2=3y2x3=2y3,使得:

y 1 y_1 y1 除以3余1、除以5余0、除以7余0;

y 2 y_2 y2 除以3余0、除以5余1、除以7余0;

y 3 y_3 y3 除以3余0、除以5余0、除以7余1;

则: x = 2 ∗ y 1 + 3 ∗ y 2 + 2 ∗ y 3 x=2*y_1+3*y_2+2*y_3 x=2y1+3y2+2y3

y 1 满足除以 5 余 0 、除以 7 余 0 y_1满足除以5余0、除以7余0 y1满足除以50、除以70,则 y 1 = 5 ∗ 7 ∗ k y_1=5*7*k y1=57k

y 1 m o d    3 = 1 ,则取 k = 2. y_1 \mod 3=1,则取k=2. y1mod3=1,则取k=2.

得到:
x = 2 ∗ 70 + 3 ∗ 21 + 2 ∗ 15 = 233 x=2*70+3*21+2*15=233 x=270+321+215=233

最后一步,x对所有除数的乘积取模运算,得到最小的满足解:

x = x m o d    ( 3 ∗ 5 ∗ 7 ) = 23 x=x \mod (3*5*7)=23 x=xmod(357)=23

本人的实现代码如下:

class Solution:
    """
    @param num: the given array
    @param rem: another given array
    @return: The minimum positive number of conditions to meet the conditions
    """
    def remainderTheorem(self, num, rem):
        # write your code here
        ls=[]
        k=1 
        for item in num:
            k*=item 
            
        for i in range(len(num)):
            elem=k//num[i]
            
            j=1 
            while elem*j%num[i]!=1: #找到最小满足的k值
                j+=1 
            ls.append(elem*j*rem[i])
        res=sum(ls)%k  
        return res

3.其他

如果存在余数相同的项,可通过将对应的除数相乘,这样少了一些项,但是好像也没节约时间。这里就不介绍了。

参考文献

[1] 知乎–中国剩余定理(CRT )
[2] 百度百科–孙子定理

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

外卖猿

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值