python实现中国剩余定理

定理内容：

题目要求：
（-1） 输入list_m,list_a
（0） 写出欧几里得算法，定义出gcd(a :int,b:int)函数：用来求取两个整数的最大公约数

（1） 首先要判断输入的list_m 中所有的元素是否两两互素，如果不是，就不能 用中国剩余定理（CRT).

（2）满足条件，第二步求出list_m中所有元素之积 m

（3）求list_M ; list_M[i] = m // list[i]

（4）求list_M_inverse; list_M_inverse[i] * list_M[i] ≡ 1 （mod list_m[i])

（5） 求list_x[i] = list_a[i] * list_M[i] * list_M_inverse[i]

（6）求出方程组的解X ≡ ∑ list_x[i] (mod m )

# 求两个数的最大公约数gcd函数
def gcd(a,b):
    if b == 0:
        return a
    else :
        return gcd(b,a % b)

 判断一个列表任意两个数是否两两互质
 def compare(list):
     for i in range(0,len(list)):
         flag = 1
         for j in range(i+1,len(list)):
             if gcd(list[i],list[j]) != 1:
                 print('不能直接使用中国剩余定理！')
                exit()
#如果满足条件，就会继续执行，否则退出程序


# 求出输入的m1,m2,..,mk 的乘积m
def product_m(list):
    m = 1
    for i in list:
        m *= i
    return m


# 求M1,M2,..,MK 的值 Mj = m / mj 并返回一个名为shang的列表
def get_divsion(list,m):
    div = []
    for i in list:
        div.append(m // i)
    return div


def get_inverse(a, m):  # 求一个数a 的逆  再模m 的值  这个函数返回的是一个值不是列表
    if gcd(a, m) != 1:
        return None
    u1, u2, u3 = 1, 0, a
    v1, v2, v3 = 0, 1, m
    while v3 != 0:
        q = u3 // v3
        v1, v2, v3, u1, u2, u3 = (u1 - q * v1), (u2 - q * v2), (u3 - q * v3), v1, v2, v3
    return u1 % m


# 求Xj   算法为：Xj = (M * M_INVERSE * a) % mj
def get_x(M:int ,M_inverse:int ,a :int,m:int):
    product_x = (M * M_inverse * a) % m
    return product_x



# 算出最终答案X = X1+X2+...Xk
def get_solution(list_m, list_a):
                                                        # compare(list_m)
    m = product_m(list_m)

    list_M = get_divsion(list_m, m)

    list_M_inverse = []
    list_X = []
    total = 0

    for i in range(0,len(list_M)):
        list_M_inverse.append(get_inverse(list_M[i], list_m[i]))

    for i in range(len(list_M)):
        list_X.append(get_x(list_M[i],list_M_inverse[i],list_a[i],m))

    for x in list_X:
        total += x

    return total % m


# 测试数据
#list_a = [2,3,2]
#list_m = [3,5,7]
#print(get_solution(list_m, list_a))

 

# 调用get_solution()函数即可使用中国剩余定理
# get_solution()函数要传入得是两个列表list_a,list_m
# 读取与输入list_a,list_m;并将其变为整型的数据在

代码的注释原本写的很多，但是那个py文件内容被我误删了，我又懒得重写注释，就随便写了点！
！！本文仅供自己积参考，如有错误请别联系我，请自行解决！！