CINTA6

小猴猴有大newnew

已于 2023-11-19 22:14:49 修改

阅读量145

点赞数 2

文章标签：美女

于 2023-11-19 22:08:57 首次发布

本文链接：https://blog.csdn.net/szy20040312/article/details/134493920

版权

运用 CRT 求解：
x ≡ 8 (mod 11)
x ≡ 3 (mod 19)
解：a=8,b=3,p=11,q=19,n=pq=209
由egcd算法可知， $pp^{−1} ≡ 1 (mod q)，11p^{-1}≡ 1 (mod 19)$ ， $p^{-1}=7$ ； $qq^{−1} ≡ 1 (mod p)，19q^{-1}≡ 1 (mod 11)$ ， $q^{-1}=7$ ，令 $y ≡ aqq^{−1} + bpp^{−1}(mod n)，y =41$ ；
运用 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，b1b1^{-1}≡ 1 (mod 5)，b1^{-1}=2；b2=3465/7=495，b2b2^{-1}≡ 1 (mod 7)，b2^{-1}=3；b3=3465/9=385，b3b3^{-1}≡ 1 (mod 9)，b3^{-1}=4；b4=3465/11=315，b4b4^{-1}≡ 1 (mod 11)，b4^{-1}=8$ 。则 $x = (1 * 693 * 2 + 2 * 495 * 3 + 3 * 385 * 4 + 8 * 315 * 4) m o d (3456) = 1731 。$
手动计算 $2000^{2019} (mod 221)$ ，不允许使用电脑或者其他电子设备。[提示：这是一
道看上去与中国剩余定理无关的计算题。]
由中国剩余定理， $221 = 17 * 13 ， 2000 = 50 * 40$ ，所以50↔（16，11），40↔（6，1），所以 $（ [16 * 6 m o d 17] ， [11 * 1 m o d 13] ） = (11, 11) = 11$ ，即 $11^{2019}mod(221)$ ，由费尔马小定理： $11^{16}≡ 1 (mod 17) ，11^{12} ≡ 1 (mod 13)，( 11^{2019}mod 17, 11^{2019}mod 13 ) = ( 11^{3}mod 17, 11^{3}mod 13 )=(5,5)=5。$
实现一个利用 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;
}