CINTA——2

CINTA——作业2

编程题

1、写一个模指数运算函数Mod_Exp，输入a、b和m，输出a^b mod m，即a的b次方模m

最简单的方法是重复地a乘b mod c，具体实现代码如下：

#include<iostream>
using namespace std;
int Mod_Exp(int a,int b,int c) {
	int p = pow(a, b);
	int d = p % c;
	return d;
}
int main() {
	int a, b, c;
	cout << "请输入a,b,c" << endl;
	cin >> a >> b >> c;
	cout << "a的b次方模c的结果是"<<Mod_Exp(a, b, c) << endl;
	return 0;
}

还有一种方法是利用中间数值，具体实现代码如下（使时间复杂度从O(n)变为O(log（n）)）：

//模指数运算——递归
#include<iostream>
using namespace std;
int rec_mod_exp(int x,int y,int p) {
	if (y == 0) return 1;
	int z = rec_mod_exp(x, y / 2, p);
	if ((y&1) == 0) //按位与
	{ return z * z % p; }
	else return x * z * z % p;
}
int main() {
	int a, b, c;
	cout << "请输入a,b,c" << endl;
	cin >> a >> b >> c;
	cout << "a的b次方模c的结果是" << rec_mod_exp(a,b,c)<<endl;
	return 0;
}

//模指数运算——迭代版本
#include<iostream>
using namespace std;
int mod_exp(int a, int b, int c) {
	int res = 1;
	while (b > 0) {
		if ((b&1) == 1) //按位与
		{
			res = (res * a) % c;
		}
		b = b / 2;
		a = (a * a) % c;
		return res;
	}
}
int main() {
	int a, b, c;
	cout << "请输入a,b,c" << endl;
	cin >> a >> b >> c;
	cout << "a的b次方模c的结果是" << mod_exp(a, b, c) << endl;
	return 0;
}

2、写一个求乘法逆元的函数Mul_Inverse，输入a和m，求a模m的乘法逆元。提示，要求只输出正整数。

#include<iostream>
using namespace std;
//扩展欧几里得算法
	int egcd(int a, int b, int &r, int &s)//通过引用类型来变相输出多个变量值
	 {
		if (b == 0)
		{
			r = 1;
			s = 0;
			return a;
		}
		int x = egcd(b, a % b, r, s);
		int temp = s;
		s = r - (a / b) * s;
		r = temp;
		return x;//返回最大公约数
	}
//计算乘法逆元（需要调用扩展欧几里得算法）
int Mul_Inverse(int a,int m) {
	int x, y;
	int temp = egcd(a, m, x, y);//调用欧几里得算法
	if (temp != 1) //最大公约数不是1，说明不互素
	{ cout << "乘法逆元不存在" << endl;
	exit(-1);
	}
	else {
		x = (x % m + m) % m;
	}
	//输出只能是正整数
	if (x < 0 || x == 0) { cout << "乘法逆元不存在" << endl; exit(1); }
	return x;
}
int main() {
	//a mod m
	int a, m;
	cout << "请输入a,m（a mod m）" << endl;
	cin >> a >> m;
	cout << "乘法逆元是:" << Mul_Inverse(a, m) << endl;
	return 0;
}

计算题

1、设p = 23 和a = 3，使用费尔马小定理计算$a^{2019}$ mod p

费尔马小定理：
设p 是素数，设a 为任意整数，且a$\not\equiv$ 0(mod p)，则
$a^{p-1}$ $\equiv$ 1 (mod p)

5、请证明13 整除$2^{70}$ + $3^{70}$。[提示：这是一道名为证明题的计算题。]

欧拉定理：
欧拉Phi 函数ϕ 对任意正整数n，返回1 到n 之间且与n 互素的整数个数，即：
ϕ(n) = |{b:1 $\le$ b < n且gcd(b,n)=1 }|

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

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