CINTA 作业三

1、实现求乘法逆元的函数，给定a和m，求a模m的乘法逆元，无解时请给出无解提示，并且只返回正整数。进而给出求解同余方程（ax = b mod m）的函数，即给定a，b，m，输出满足方程的x，无解给出无解提示。

//invs 函数
//输入：两个整数 a 和 m
//输出：a 在 模 m 下的乘法逆元，当无解的时候返回 0

//svlCeq函数
//输入：同余方程式 ax \equiv b (mod m) 中的 三个整型系数 a，b，m 和一个用于存储 x 的所有解的整型数组的地址
//输出：返回同余方程式的解的个数，解保存于数组 ans 中  

int egcd(int a, int b, int *r, int *s){
	if(b == 0) return a;
	
	int r0 = 1, s0 = 0, r1 = 0, s1 = 1;
	int q = 0, rest = 0;
	int t_r0,t_s0;
	while(b){
		q = a/b, rest = a%b;
		
		t_r0 = r0, t_s0 = s0;
		a = b, r0 = r1, s0 = s1;
		 
		b = rest, r1 = t_r0 - q*r1, s1 = t_s0- q*s1;
	}
	
	*r = r0, *s = s0;
	return a;		 
}

int invs(int a, int m){
	int g = 0, ivs = 0, s = 0;
	g = egcd(a, m, &ivs, &s);
	ivs %= m;
	if(g != 1) return 0;
	return ivs;
} 

int svlCeq(int a, int b, int m,int *ans){
	// 四个变量依次表示 方程 ar + ms = g 中的第一个系数，~ 中的第二个系数，a 和 m 的最大公因子，同余式的某一个解 x0 
	int r = 0, s = 0, g = 0, x0 = 0;
	g = egcd(a, -m , &r, &s);

	//当方程 ar + ms = g 不能整除 b 的时候，方程没有解，返回-1	
	if(b % g != 0) return -1;

	while(r <= 0) r += m;
	x0 = b*r/g ;
	//根据线性同余方程，如果 g|c，同余式 ax \equiv b (mod m) 恰好有 g 个不同的解 
	for(int i = 0; i < g; i++)
		ans[i] = ( x0 + i*m/g )% m;
 
	return g;
}

2、实现模指数运算的函数，给定x、y和m，求x的y次方模m。

int mod_exp(int x, int y, int m) {
	int res = 1;
	long long int ccl = 0;

	while (y > 0) {
		if (y & 1) {
			ccl = res * x;
			res = ccl % m;
		}
		y /= 2;
		ccl = x * x;
		x = ccl % m;
	}

	return res;
}

3、设p = 23和a = 5，使用费尔马小定理计算a^{2020} mod p?

解答：

由于 p = 23， 我们可以通过费尔马小定理的到 $5^{22}$ mod 23 与 1 同余，

由此我们可以得到 $5^{2020}\equiv 5^{18}$ mod 23。

在二进制下 18 = 10010，那么我们可以得到 $5^{18}\equiv 5^{2^4}\ast 5^{2^1}$ mod 23，

经过简单的运算我们可以算出结果等于 $6$。

4、使用欧拉定理计算2^{100000} mod 55。

解答：

由于 55 是合数，我们欧拉函数的性质得到 $\varphi (55)=\varphi (5)\ast \varphi (11)=40$。

则通过欧拉定理我们可以得到 $2^{100000}\equiv 2^{40}\equiv 1$ mod 55。

由此结果为 $1$ 。

5、手动计算7^{1000}的最后两个数位等于什么？

解答：

由题目描述我们可以认为题目在求 $7^{1000}$ mod 100 的结果。

通过手动计算我们可以得到 末尾两个数位是以 4 个数为周期循环出现的，他们依次为 07, 49, 43, 01。

经过简单推理我们可以得到最终的结果是 $01$。