CINTA6

1. 运用 CRT 求解：
   x ≡ 8 (mod 11)
   x ≡ 3 (mod 19)
   解：a=8,b=3,p=11,q=19,n=pq=209
   由egcd算法可知，$p{p}^{-1}\equiv 1(modq)，11p^{-1}\equiv 1(mod19)$，$p^{-1}=7$；$q{q}^{-1}\equiv 1(modp)，19q^{-1}\equiv 1(mod11)$，$q^{-1}=7$，令$y\equiv aq{q}^{-1}+bp{p}^{-1}(modn)，y=41$；

2. 运用 CRT 求解：
   x ≡ 1 (mod 5)
   x ≡ 2 (mod 7)
   x ≡ 3 (mod 9)
   x ≡ 4 (mod 11)
   解：对于前两个，$a=1,b=2,c=3,d=4,p=5,q=7,s=9,t=11,M=3465，b1=3465/5=693，b1b{1}^{-1}\equiv 1(mod5)，b{1}^{-1}=2；b2=3465/7=495，b2b{2}^{-1}\equiv 1(mod7)，b{2}^{-1}=3；b3=3465/9=385，b3b{3}^{-1}\equiv 1(mod9)，b{3}^{-1}=4；b4=3465/11=315，b4b{4}^{-1}\equiv 1(mod11)，b{4}^{-1}=8$。则$x=(1\ast 693\ast 2+2\ast 495\ast 3+3\ast 385\ast 4+8\ast 315\ast 4)mod(3456)=1731。$

3. 手动计算 $2000^{2019}(mod221)$，不允许使用电脑或者其他电子设备。[提示：这是一
   道看上去与中国剩余定理无关的计算题。]
   由中国剩余定理，$221=17\ast 13，2000=50\ast 40$，所以50↔（16，11），40↔（6，1），所以$（[16\ast 6mod17]，[11\ast 1mod13]）=(11,11)=11$，即$11^{2019}mod(221)$，由费尔马小定理：$11^{16}\equiv 1(mod17)，11^{12}\equiv 1(mod13)，(11^{2019}mod17,11^{2019}mod13)=(11^{3}mod17,11^{3}mod13)=(5,5)=5。$

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

在这里插入代码片#include<iostream>
using namespace std;
int inverse(int m,int b)
{
	int k = 1;
	while (1)
	{
		if ((m * k + 1) % b == 0)
		{
			return ((m * k + 1) / b);
			break;
		}
		k++;
	}
}
int main()
{
	cout << "输入同余方程个数及同余方程（格式为x mod c = a）" << endl;
	int num,M=1;
	int input1[10],b[10],ib[10],ans=0; int input2[10]; string mod; char x, e;
	cin >> num;
	for (int i = 0; i < num; i++)
	{
		cin >> x >> mod >> input1[i] >> e >> input2[i];
	}
	for (int i = 0; i < num; i++)
	{
		M *= input1[i];
	}
	for (int i = 0; i < num; i++)
	{
		b[i] = M / input1[i];
	}
	for (int i = 0; i < num; i++)
	{
		ib[i] = inverse(input1[i], b[i]);
	}
	//for (int i = 0; i < num; i++)
	//{
	//	cout << ib[i] << endl;
	//}
	for (int i = 0; i < num; i++)
	{
		ans += input2[i] * b[i] * ib[i];
	}
	ans %= M;
	cout << ans;
	system("pause");
	return 0;
}

测试结果