洛谷 P1495 【模板】中国剩余定理（CRT）/ 曹冲养猪

本来期末复习时间不多了，想挪用点代码完成洛谷上的作业题。但是这题搜到的的代码就是一直有一个测试点过不了关。加上c++一年半没碰，Java不会，只能用自己弱弱的python＋b站视频搞出来一个代码。很开心，终于过了。视频链接和代码如下。

乘法逆元的概念解释：【【乘法逆元】模运算下除以一个整数，就相当于乘以这个整数的乘法逆元】https://www.bilibili.com/video/BV1kv4y1P7WV?vd_source=89a434dcd0263bf94b2f9f7ca0f5c7f7
中国剩余定理解题方法：【【中国剩余定理】如何求解一次同余方程组】https://www.bilibili.com/video/BV1XW4y1k7us?vd_source=89a434dcd0263bf94b2f9f7ca0f5c7f7
代码：

def mod_inverse(a, m):
    # 扩展欧几里得算法求乘法逆元
    def egcd(a, b):
        if a == 0:
            return (b, 0, 1)
        else:
            g, y, x = egcd(b % a, a)
            return (g, x - (b // a) * y, y)

    g, x, _ = egcd(a, m)
    if g != 1:
        raise Exception('Modular inverse does not exist')
    else:
        return x % m

n = int(input())
a = [0]*n  # 初始化a列表
b = [0]*n  # 初始化b列表
product_a = 1  # 初始化a的累积乘积
for i in range(n):
    a[i], b[i] = map(int, input().split())
    product_a *= a[i]  # 累积a[i]

results = 0
for i in range(n):
    c_i = product_a // a[i]  # 计算c[i]
    inv_c_i = mod_inverse(c_i, a[i])  # 计算c[i]在模a[i]下的乘法逆元
    results += c_i * inv_c_i * b[i] # 累加
results = results % product_a  # 求results除以product_a的余数
print(results)

 炫耀一下：