CINTA第六次作业

第十章

题1

题目：
运用 CRT 求解：
$x\equiv 8(\mathrm{mod}11)$
$x\equiv 3(\mathrm{mod}19)$
解答：
$19$ 模 $11$ 的乘法逆元为 $7$, $11$ 模 $19$ 的乘法逆元为 $7$,
利用 CRT 有：
$x\equiv 8\times 19\times 7+3\times 11\times 7\equiv 41(\mathrm{mod}201)$

题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)$
解答：
设上述第 $i$ 条同余方程为 x ≡ a i ( m o d m i ) , i ∈ { 1 , 2 , 3 , 4 } x\equiv a_i\pmod {m_i}, i\in\{1,2,3,4\} x≡ai​(modmi​),i∈{1,2,3,4}

$a_i$	$b_i=({\prod }_1^n{m}_i)/m_i$	$b_i$ 模 m 的乘法逆元
1	693	2
2	495	3
3	385	4
4	315	8

故
$$x\equiv \sum _{i=1}^n{a}_i{b}_i{b}_i^{-1}\equiv 1731(\mathrm{mod}3465)$$

题3

题目：手动计算 $2000^{2019}(\mathrm{mod}221)$，不允许使用电脑或者其他电子设备。[提示：这是一道看上去与中国剩余定理无关的计算题。]
解答：
利用群同构转化即可.

设 $n=pq$，$p>1$ 和 $q>1$ 是两个互素的正整数，则有：$Z_n^{\ast }\cong Z_p^{\ast }\times Z_q^{\ast }.$

由于 $221=13\times 17$，再由上述命题可知：$Z_{221}^{\ast }\cong Z_{13}^{\ast }\times Z_{17}^{\ast }$.
故：$2000\leftrightarrow (2000\mathrm{mod}13,2000\mathrm{mod}17)\leftrightarrow (11,11)$.
又 $(11,11{)}^{2019}=(11^{12\times 168+3},11^{16\times 126+3})=(11^3,11^3)=(5,5)$，其中第二个等号由费尔马小定理可知.
不难发现 $5\leftrightarrow (5,5)$，故答案为 $5$.

题7

题目：实现一个利用 CRT 求解同余方程的程序（Python 或者 C 语言都可以）。
代码如下：

//一共n条同余方程，同余方程为 x≡a(mod b)
int n, a[15], b[15];

//(套egcd模板)求x的乘法逆元
int mul_inverse(int x, int a)
{
    int r0 = 1, s0 = 0, r1 = 0, s1 = 1;
    while(a)
    {
        int tr = r0, ts = s0, t = x % a, q = x / a;
        r0 = r1, s0 = s1, x = a;
        r1 = tr - q * r1, s1 = ts - q * s1, a = t;
    }
    return r0;
}

//套crt公式
long long crt(int a[], int b[], int n)
{
    long long ans = 0;
    long long B = 1;//B是同余方程所有b的乘积
	for(int i = 0; i < n; i ++ ) B *= b[i];
    for(int i = 0; i < n; i ++ )
    {
        ans += a[i] * (B / b[i]) * mul_inverse(B / b[i] % b[i], b[i]) % B;
    }
    return ans % B;
}