快速指数取模运算与用扩展欧几里得算法求解最大公约数和求乘法逆元.docxVIP

一、实验1.1源代码： #includc stdio.h #include stdlib.h #include uiostreamM using namespace std; void Modc(int a, int b, int n) { int c=l; do{ if(a%2==0) { a=a/2; b=(b*b)%n; } else { a=a-1; c=(b*c)%n; } }whilc(a!=0); cout?取余结果为:M?c?endl; } void main () { int a, b, n; cout?输入格式范例 bAa mod n?endl; cout? u?endl; cout?H请输入指数a： n?endl; cin?a; coutvv请输入该基数b： ?endl; cin?b; coutvv请输入被除数n： ?endl; cin?n; Mode(a,b,n); 二、实验效果图: 回 S3 nod n 请输入*旨薮a： 37 请输入该基数b： 30 请输入被除数n： 77 取余结果为：2 Press any key to continue nr 实验1.2用扩展欧几里得算法求解最大公约数和求乘法逆元 一、实验1.2源代码： #include stdio.h int extended_Gcd(int a,int b, int x, int y) 〃求最大公约数 { if (b == 0) { x = 1; y = o； return a; } else { int ged = extended_Gcd(b, a%b, x, y); int t = x; x 二 y; y = t - (a / b) * y; return ged; } } int extended_Ivn(int f, int d, int ^result) 〃求乘法逆元 { int xl, x2, x3, yl, y2, y3, tl, t2, t3, q; xl = y2 = 1; x2 = yl =0; x3 = (f= d)?f:d; y3 = (f = d) ? d : f; while (1) { if (y3 == 0) { Result = x3; //两个数不互索则result为两个数的最大公约数，此时返回值 为零 return 0; if(y3 == 1) ^result = y2; 〃两个数互素则resutl为其乘法逆元，此吋返冋值为1 return 1; q = x3 /y3; tl = xl - q*yl; t2 = x2 - q*y2; t3 = x3 - q*y3; xl = yl; x2 = y2; x3 = y3; yl =tl; y2 = t2； y3 = t3； int main() { int x, y, z ; int a, b; int *p, *q; p = x; q = y; z = 0; printf(”请输入两个数:”); scanf(,,%d%du, a, b); if (extended_Gcd(a,b? *p, *q) == 1) { extended_Ivn(a, b, z); printf(%d 和％＜1 互素,乘法的逆元是:%d\n, a, b, z); } else { z=extended_Gcd(a,b, *p, *q); printf(n%d和％＜1不互素，最大公约数为:%d\nM, a, b, z); } return 0; 二、实验效果图: El C:\Users\Sorrel\Desktop\^f^$；^^\Debug\^^l.2.exe ■3 ?C:\Users\Soirel\Desktop\新建文件夹\Debug\^验l?2.exe” 1 = If a If 23 1 请输入两个数「1234 4321 1234和4321互素，乘法的逆兀是：-1082 二 Press any key to continue ■3 C:\Users\Sorrel\Desktop\Wr^^{i^\Debug\^^1.2.exe, 劇磁,诙元是:剧 Press any key to continue