CINTA作业三：同余、模指数、费尔马小定理、欧拉定理

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

int multiply(int a, int b)
{
	int r1 = 0, r2 = 0, s1 = 0, s2 = 0;
	while (b)
	{
		int mod1 = a % b;
		int div1 = a / b;
		a = b;
		b = mod1;
		int temp_r1 = r1;
		int temp_s1 = s1;
		r1 = r2;
		s1 = s2;
		r2 = temp_r1 - r1 * div1;
		s2 = temp_s1 - s1 * div1;
	}
	if (r1 >= b)
	{
		cout << "无解" << endl;
		return -1;
	}
	return r1;
}

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

int mod_exp(int x, int y, int p) //在主函数中指定
{
	int res = 1;//保留数
	while (y > 0)
	{
		if ((y & 1) == 1)res = (res * x) % p;
		y = y / 2;
		x = (x * x) % p;
	}
	return res;
}

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

5(2020)mod23=5((23−1)∗91)+18mod23=518mod23=6

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

根据欧拉定理：

则：2(100000)mod55=2(40∗2500)mod55=1 

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

7(1000)mod100=7(40∗25)mod100=1

因此为01