CINTA第六次作业

---

---

一、运用CRT求解：x ≡ 8 (mod 11）； x ≡ 3（mod 19）

利用中国剩余定理，解法如下：
1、记a=8，b=3，p=11，q=19和n=pq=209
2、使用egcd算法求解p^-1和q^-1;
它们分别使得pp^-1 ≡1（mod q）和qq^-1 ≡1(mod p);
也就是1111^-1≡1(mod 19)得11^-1=7；
同样1919^-1 ≡1(mod 11)得19^-1-=7;
3、令y≡aqq^-1+bpp^-1(mod n)
设x=aqq^-1+bpp^-1
x=8(mod 11)+ 3(mod 19)=8197+3117=1064+231=1295
y=1295(mod 209)=41
4、验证可知y是正确解

二、运用CRT求解：x ≡ 1 (mod 5）；x ≡ 2 (mod 7）；x ≡ 3 (mod 9）；x ≡ 4 (mod 11）

利用中国剩余定理—推广版，解法如下：
1、记a1=1，a2=2，a3=3，a4=4；
m1=5，m2=7，m3=9，m4=11；
M=m1m2m3m4= 3465
2、使用egcd算法求解b1^-1,b2^-1,b3^-1,b4^-1
b1=m2m3m4=693
b2=m1m3m4=495
b3=m1m2m4=385
b4=m1m2m3=315
b1b1^-1≡1(mod 5)
b2b2^-1≡1(mod 7)
b3b3^-1 ≡1(mod 9)
b4b4^-1≡1(mod 11)
解得：b1^-1=2,b2^-1=3,b3^-1=4,b4^-1=8
3、x=$$\sum _{i=0}^{n-1}{a}_i{b}_i{b}_i^{-1}(\mathrm{mod}M)$$``
所以x=16932+24953+338584+43158(mod 3465)=1731

三、手动计算2000²⁰¹⁹(mod 221)，不允许使用电脑或其他电子设备

211=13*17
13和17都为素数
$$Z_{13}^{\ast }\times Z_{17}^{\ast }$$
2000↔(11,11)
所以：
(11,11)²⁰¹⁹=([11²⁰¹⁹mod 13],[11²⁰¹⁹mod 17])
11²⁰¹⁹mod 13=11^12*168+3mod13=11³mod 13 =5；
11²⁰¹⁹mod 17=11^16*126+3 mod 17=11³mod 17=5；
综上：
(11,11)²⁰¹⁹=（5,5）=5

七、实现一个利用CRT求解同余方程的程序

#include <iostream>
using namespace std;
const int N=1e5+5;
int n,m[N],a[N]; 
 
int egcd(int a,int b,int &x,int &y)
{
	if(b=0)
	{
		x=1;
		y=0;
		return a;
	}
	int d = egcd(b,a%b,x,y);
	int t=x;
	x=y;
	y=t-a/b*x;
	return d;
 } 

int CRT()
{
	int ans=0,M=1,x,y;
	for(int i=1;i<=n;i++)
	{
		M*=m[i];           //M为m1*m2*...mn 
	}
	for(int i=1;i<=n;i++)
	{
		int bi=M/m[i];     //bi=除了bi外的其他b相乘 
		egcd(bi,m[i],x,y);   // 
		ans=(ans+a[i]*x*bi)%M;
	}
	return (ans+M)%M; 
}


int main(){
	cout<<"请输入n值为：";
	cin>>n;
	cout<<"请依次输入m的值为："<<endl; 
	for(int i=1;i<=n;i++)
	{
		cin>>m[i];
	}
	cout<<"请依次输入a的值为："<<endl; 
	for(int i=1;i<=n;i++)
	{
		cin>>a[i];
	}
    cout<<"结果为："<<CRT();
    return 0;
}

总结

有错误待改正