中国剩余定理(模数互质)

最新推荐文章于 2022-05-17 23:21:57 发布

前方是否可导？

最新推荐文章于 2022-05-17 23:21:57 发布

阅读量192

点赞数

分类专栏：其他类型

本文链接：https://blog.csdn.net/weixin_44110537/article/details/107237976

版权

其他类型专栏收录该内容

22 篇文章 0 订阅

订阅专栏

from functools import reduce
import sympy
def CRT(n, a):
    sum = 0
    prod = reduce(lambda a, b: a * b, n)#reduce(function,n):将数据集合n(可以是列表,也可以是元组等)的前两个数据做function函数运算,然后将结果与n中的后一个元素做function运算,最后会返回一个结果.

    for n_i, a_i in zip(n, a):#zip(a,b)将a,b中的数据打包成元组的列表[(a1,b1),(a2,b2),(a3,b3)......]
        p = prod // n_i
        sum += a_i * sympy.invert(p, n_i) * p
    return int(sum % prod)
#,其中n是模数,a是余数.
#中国剩余定理可以得到一个最小的满足方程组的解.
'''
假设mult[i]表示除n[i]的其他n中的元素的积
sum(a[i]*invert(mult[i],n[i])*mult[i])%(mult[i]*n[i])其中i=0,1,2,3,4,5......
'''