【CINTA06】CINTA第六次作业

【CINTA06】CINTA第六次作业

题1

$运用CRT求解：x\equiv 8(mod11)x\equiv 3(mod19)$

$$\begin{gathered} 解：a=8,p=11,b=3,q=19,n=pq=209 \\ {p}^{-}{1}^{=}7,q^{-}{1}^{=}7 \\ x=(aq{q}^{-}{1}^{+}bp{p}^{-}{1}^{)}modn=41 \end{gathered}$$

题2

$运用CRT求解：x\equiv 1(mod5)x\equiv 2(mod7)x\equiv 3(mod9)x\equiv 4(mod11)$

$$\begin{gathered} 解：M=5\ast 7\ast 9\ast 11=3465 \\ 记b0=M/5=693，b1=M/7=495，b2=M/9=385，b3=M/11=315 \\ bibi{}^{-}{1}^{=}1mod{m}_i \\ b0{}^{-}{1}^{=}2,b1{}^{-}{1}^{=}3，b2{}^{-}{1}^{=}4，b3{}^{-}{1}^{=}8 \\ x=\sum aibibi{}^{-}{1}^{m}odM=1731 \end{gathered}$$

题3

$手动计算20002019(mod221)，不允许使用电脑或者其他电子设备。[提示：这是一道看上去与中国剩余定理无关的计算题。]$

$解：221=13\ast 1720002019《=》（11，7）x=11mod13,x=7mod17a=11,p=13,b=7,q=17p^{-}{1}^{=}4,q^{-}{1}^{=}10x=(aq{q}^{-}{1}^{+}bp{p}^{-}{1}^{)}mod221=2420002019(mod221)=24$

题7

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

$求逆元思路：由费马小定理可知：p^q-1^{=}1modq{p}^q-2^{=}{p}^{-}{1}^{m}odq{p}^{-}{1}^{=}{p}^q-2^{m}odq$

#include<iostream>
using namespace std;

int inversion(int x, int y)
{
	int a = 1;
	if (y == 2)
	{
		if (x % 2 != 0)
			return a;
	}
	for (int i = 0; i <y-2 ; i++)
	{
		a = a * x % y;
	}
	return a;
}

int CRT(int a,int p,int b,int q)
{
	int s = inversion(p, q);
	int r = inversion(q, p);
	int x = a * q * r + b * p * s;
	x = x % (p * q);
	return x;
}


int main()
{
	int a, p, b, q;
	cout << "请输入a,p,b,q(p,q互素）：" << endl;
	cin >> a >> p >> b >> q;
	int x = CRT(a, p, b, q);
	cout << "求得：";
	cout << x << endl;
}