中国剩余定理与扩展中国剩余定理

最新推荐文章于 2024-05-30 18:05:01 发布

葛嘣

最新推荐文章于 2024-05-30 18:05:01 发布

阅读量133

点赞数

文章标签： python 算法 Powered by 金山文档

本文链接：https://blog.csdn.net/m0_53437775/article/details/129654434

版权

一、中国剩余定理

代码模板

n = int(input())
# m 是 模 数，r 是 余 数
m,r = [0] * n,[0] * n  
s = 1
t = [0] * n

for i in range(n):
    m[i],r[i] = map(int,input().split())
    s *= m[i]
for i in range(n):
    t[i] = s // m[i]

# 扩 展 欧 几 里 得 求 逆 元
def exgcd(a,b):
    if b == 0:
        return 1,0,a
    x,y,gcd = exgcd(b,a % b)
    x,y = y,x - (a // b) * y
    return x,y,gcd
def inv(a,b):
    x,y,gcd = exgcd(a,b)
    return x % b

res = 0
for i in range(n):
    res += (t[i] * inv(t[i],m[i]) * r[i]) % s
    res %= s
print(res)

例题

https://www.luogu.com.cn/problem/P1495

原理

略

二、中国剩余定理扩展定理

代码模板

n = int(input())
m,r = [0] * n,[0] * n
for i in range(n):
    m[i],r[i] = map(int,input().split())
def exgcd(a,b):
    if b == 0:
        return 1,0,a
    x,y,gcd = exgcd(b,a % b)
    x,y = y,x - (a // b) * y
    return x,y,gcd
def excrt():
    if n == 1:
        return r[0]
    for i in range(n):
        x,y,gcd = exgcd(m[0],m[i])
        if (r[i] - r[0]) % gcd != 0:
            return -1

        # 其 实 就 是 计 算 得 到 特 解
        x = x * ((r[i] - r[0]) // gcd)
        # 同 时 这 个 数 有 可 能 为 负 数，
        # 要 映 射 成 最 小 正 整 数
        x = x % (m[i] % gcd)
  
        r[0] = x * m[0] + r[0]
        # 将 x = r[0](mod m[0])
        # 与 x = r[i](mod m[i])
        # 合 并 成 新 的 x = r[0](mod m[0])
        m[0] = m[0] * m[i] // gcd  # 最小公倍数
        r[0] %= m[0]
    return r[0]
print(excrt())