【CINTA第6次作业】

最新推荐文章于 2024-06-28 14:48:01 发布

2201_75315239

最新推荐文章于 2024-06-28 14:48:01 发布

阅读量900

点赞数 25

文章标签：算法

本文链接：https://blog.csdn.net/2201_75315239/article/details/134771021

版权

CINTA第6次作业

在这里插入图片描述
解：
$\begin{aligned} &令a=8,b=3,p=11,q=19\\ &\because p,q互素\\ &\therefore 根据egcd算法解pp^{-1}\equiv1\mod q和qq^{-1}\equiv1\mod p\\ &得p^{-1}=7,q^{-1}=7\\ &\because n=pq=209\\ &\therefore x\equiv aqq^{-1}+bpp^{-1}\mod n\\ &带入数据解得x=41\\ \end{aligned}$

在这里插入图片描述
解：
$\begin{aligned} &令a_0=1,a_1=2,a_2=3,a_3=4,m_0=5,m_1=7,m_2=9,m_3=11\\ &\therefore M=m_0m_1m_2m_3=3465\\ &\therefore 记b_0=M/m_0=693,b_1=M/m_1=495,b_2=M/m_2=385,b_3=M/m_3=315\\ &\because b_0,m_0互素\\ &\therefore 根据egcd算法解b_0b_0^{-1}\equiv1\mod m_0,解得b_0^{-1}=2\\ &同理可得b_1^{-1}=3,b_2^{-1}=4,b_3^{-1}=8\\ &\therefore x\equiv a_0b_0b_0^{-1}+a_1b_1b_1^{-1}+a_2b_2b_2^{-1}+a_3b_3b_3^{-1}\mod M\\ &带入数据解得x=1731\\ \end{aligned}$

在这里插入图片描述
解：
$\begin{aligned} &\because n=221=13*17,且13、17为素数\\ &\therefore \mathbb{Z}_n^*\cong \mathbb{Z}_p^* \times\mathbb{Z}_q^*\\ &\therefore 2000\longleftrightarrow(11,11)\\ &\therefore (11,11)^{2019}=([11^{2019}\mod 13],[11^{2019}\mod 17])\\ &下面分别对11^{2019}\mod 13和11^{2019}\mod 17进行求解\\ &根据费尔马小定理,可得11^{12}\equiv1\mod 13,11^{16}\equiv1\mod 17\\ &\therefore 11^{2019}\mod 13=11^{168*12+3}\mod 13=11^3\mod 13=5\\ &同理11^{2019}\mod 17=11^{126*16+3}\mod 17=11^3\mod 17=5\\ &\therefore (11,11)^{2019}=([11^{2019}\mod 13],[11^{2019}\mod 17])=(5,5)\\ &\therefore 2000^{2019}\mod 221=5 \end{aligned}$

在这里插入图片描述

def exgcd(a, b):
    if b == 0:
        return 1, 0, a
    x, y, q = exgcd(b, a % b)
    x, y = y, (x - a // b * y)
    return x, y, q

def CRT(a, n):
    sum = 0
    p = 1
    for n_i in n:
        p *= n_i
    for a_i, n_i in zip(a, n):
        q = p // n_i
        sum += a_i * q * exgcd(q, n_i)[0]
    return sum % p
def main():
    a = [1, 2, 3, 4]
    n = [5, 7, 9, 11]
    ans = CRT(a, n)
    print("运用CRT解同余方程：")
    for i in range(len(a)):
        print("x≡", a[i], " mod ", n[i])
    print("x =", ans)
if __name__ == '__main__':
    main()