1.写一个模指数运算函数Mod_Exp,输入a、b和m,输出a^b mod m,即a的b次方模m
#include<iostream>
using namespace std;
int Mod_Exp(int a,int b,int m)
{
int c = 1;
for (int i = 0; i < b; i++)
c *= a;
int n = c % m;
return n;
}
int main()
{
int a, b, m;
cout << "请输入a值、b值、m值:" << endl;
cin >> a >> b >> m;
cout << "a^b mod m的值为:" << endl;
cout<< Mod_Exp(a, b, m);
return 0;
}
运行结果:
2.写一个求乘法逆元的函数Mul_Inverse,输入a和m,求a模m的乘法逆元。提示,要求只输出正整数
#include<iostream>
using namespace std;
int gcd(int a,int m)
{
while (m != 0)
{
int temp = m;
m = a % m;
a = temp;
}
return a;
}
int Mul_Inverse(int a, int m)
{
int b = gcd(a, m);
if (b == 1)
{
int c = (m + 1) / a;
return c;
}
else cout << "逆元不存在!" << endl;
}
int main()
{
int a, m;
cout << "请输入a值和m值:";
cin >> a >> m;
cout << "a mod m 的乘法逆元为:";
cout << Mul_Inverse(a, m);
return 0;
}
运行结果
第四章习题:
1.设 p = 23 和 a = 3,使用费尔马小定理计算 a^ 2019 mod p?
解:
a^2019=a^(22*91+17)≡a^17 mod 23=3^17 mod 23 ≡3^1001≡3^8*3≡16
5. 请证明 13 整除 2 ^70 + 3^70。[提示:这是一道名为证明题的计算题。]
解:
由费马小定理:
2^12≡1 mod 13
3^12≡1 mod 13
2^70=2^(12*5+10)≡2^10≡9*4≡10 mod 13
3^70=3^(12*5+10)≡3^10≡9*9≡3 mod 13
2^70+3^70≡13≡0 mod 13
所以,得证
6. 使用欧拉定理计算 2 ^100000 mod 55
解:
⌀(55)=⌀(5)⌀(11)=4*10=40
因为,gcd(2,55)=1
所以,2^⌀(55)≡2^40≡1 mod 55
2^100000=2^(40*2500)≡1 mod 55
8. 手动计算 7 ^1000 的最后两个数位等于什么?
解:
计算最后两位数 即求7^1000 mod100
⌀(100)=⌀(4*25)=⌀(2^2*5^2)=⌀(2^2)⌀(5^2)=(2^2-2)*(5^2-5)=40
因为,gcd(7,100)=1
7^1000≡7^(40*25)=(7^40)^25≡(7^⌀(100))^25≡1 mod 100
所以,7^1000 mod100 的最后两位数等于01