CINTA 第六次作业

$1.运用CRT求解：$
$$x\equiv 8(\mathrm{mod}11)$$
$$x\equiv 3(\mathrm{mod}19)$$
$解：运用CRT得$
$$x=(8\ast 19\ast 19^{-1}+3\ast 11\ast 11^{-1})(\mathrm{mod}11\ast 19)$$
$其中，19\ast 19^{-1}\equiv 1(\mathrm{mod}11),11\ast 11^{-1}\equiv 1(\mathrm{mod}19),得$
$19^{-1}=7,11^{-1}=7$
$$x=(8\ast 19\ast 7+3\ast 11\ast 7)(\mathrm{mod}11\ast 19)=1295(\mathrm{mod}209)$$
$得，x=41$

$2.运用CRT求解：$
$$x\equiv 1(\mathrm{mod}5)$$
$$x\equiv 2(\mathrm{mod}7)$$
$$x\equiv 3(\mathrm{mod}9)$$
$$x\equiv 4(\mathrm{mod}11)$$
$解：运用CRT得$
$x=(1\ast (7\ast 9\ast 11)\ast (7\ast 9\ast 11{)}^{-1}+2\ast (5\ast 9\ast 11)\ast (5\ast 9\ast 11{)}^{-1}+3\ast (5\ast 7\ast 11)\ast (5\ast 7\ast 11{)}^{-1}+4\ast (5\ast 7\ast 9)\ast (5\ast 7\ast 9{)}^{-1})(\mathrm{mod}5\ast 7\ast 9\ast 11)$
$其中，(7\ast 9\ast 11)\ast (7\ast 9\ast 11{)}^{-1}\equiv 1(\mathrm{mod}5),(5\ast 9\ast 11)\ast (5\ast 9\ast 11{)}^{-1}\equiv 1(\mathrm{mod}7)，$
$(5\ast 7\ast 11)\ast (5\ast 7\ast 11{)}^{-1}\equiv 1(\mathrm{mod}9)，(5\ast 7\ast 9)\ast (5\ast 7\ast 9{)}^{-1}\equiv 1(\mathrm{mod}11)$
$得(7\ast 9\ast 11{)}^{-1}=2，(5\ast 9\ast 11{)}^{-1}=3，(5\ast 7\ast 11{)}^{-1}=4，(5\ast 7\ast 9{)}^{-1}=8$
$$x=(1\ast 693\ast 2+2\ast 495\ast 3+3\ast 385\ast 4+4\ast 315\ast 8)(\mathrm{mod}3465)$$
$得，x=1731$

$3.手动计算2000^{2019}(\mathrm{mod}221)。$

$解：221=13\ast 17$
$根据费马小定理，a^{p-1}\equiv 1(\mathrm{mod}p)，所以有$
$$2000^{12}\equiv 1(\mathrm{mod}13)$$
$$2000^{16}\equiv 1(\mathrm{mod}17)$$
$而2019\equiv 3(\mathrm{mod}12)，2019\equiv 3(\mathrm{mod}16),即$
$$2000^{2019}\equiv 2000^3(\mathrm{mod}13)$$
$$2000^{2019}\equiv 2000^3(\mathrm{mod}17)$$
$2000=2\ast 10^3，而10(\mathrm{mod}13)=-3$
$所以，2000^3(\mathrm{mod}13)=2\ast (-3{)}^3(\mathrm{mod}13)=-54(\mathrm{mod}13)=11$
$2000=20^2\ast 5，而20(\mathrm{mod}17)=3$
$2000^3(\mathrm{mod}17)=(3{)}^2\ast 5(\mathrm{mod}17)=45(\mathrm{mod}17)=11$
$即$
$$2000^{2019}\leftrightarrow (11,11)$$
$2000^{2019}=11$

$7.实现一个利用CRT求解同余方程的程序（Python或者C语言都可以）。$

def egcd(a, b):
    if b == 0:
        return (a, 1, 0)
    else:
        gcd, x, y = egcd(b, a % b)
        return (gcd, y, x - (a // b) * y)

def CRT(equations):
    N = 1
    for a, m in equations:
        N *= m

    x = 0
    for a, m in equations:
        N_i = N // m
        _, s, _ = egcd(N_i,m)
        x += a * s * N_i

    return x % N

if __name__ == "__main__":
    equations = []
    n = int(input("请输入方程个数: "))
    for i in range(n):
        a = int(input("请输入 a{}: ".format(i+1)))
        m = int(input("请输入 m{}: ".format(i+1)))
        equations.append((a, m))

    result = CRT(equations)
    print("方程的解为:", result)