快速指数取模运算与用扩展欧几里得算法求解最大公约数和求乘法逆元_第1页
快速指数取模运算与用扩展欧几里得算法求解最大公约数和求乘法逆元_第2页
快速指数取模运算与用扩展欧几里得算法求解最大公约数和求乘法逆元_第3页
快速指数取模运算与用扩展欧几里得算法求解最大公约数和求乘法逆元_第4页
快速指数取模运算与用扩展欧几里得算法求解最大公约数和求乘法逆元_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

实验1.1快速指数取模运算一、实验1.1源代码：#include"stdio.h"#include"stdlib.h"#include"iostream"usingnamespacestd;voidMode(inta,intb,intn){ intc=1; do{ if(a%2==0) { a=a/2; b=(b*b)%n; } else { a=a-1; c=(b*c)%n; } }while(a!=0); cout<<"取余结果为："<<c<<endl;}voidmain(){ inta,b,n; cout<<"输入格式范例b^amodn"<<endl; cout<<"........................................................."<<endl; cout<<"请输入指数a："<<endl; cin>>a; cout<<"请输入该基数b："<<endl; cin>>b; cout<<"请输入被除数n："<<endl; cin>>n; Mode(a,b,n); }实验效果图：实验1.2用扩展欧几里得算法求解最大公约数和求乘法逆元实验1.2源代码：#include<stdio.h>intextended_Gcd(inta,intb,int&x,int&y)//求最大公约数{if(b==0){x=1;y=0;returna;}else{intgcd=extended_Gcd(b,a%b,x,y);intt=x;x=y;y=t-(a/b)*y;returngcd;}}intextended_Ivn(intf,intd,int*result)//求乘法逆元{intx1,x2,x3,y1,y2,y3,t1,t2,t3,q;x1=y2=1;x2=y1=0;x3=(f>=d)?f:d;y3=(f>=d)?d:f;while(1) {if(y3==0) {*result=x3;//两个数不互素则result为两个数的最大公约数，此时返回值为零return0;}if(y3==1) {*result=y2;//两个数互素则resutl为其乘法逆元，此时返回值为1return1;}q=x3/y3;t1=x1-q*y1;t2=x2-q*y2;t3=x3-q*y3;x1=y1;x2=y2;x3=y3;y1=t1;y2=t2;y3=t3;}}intmain(){intx,y,z; inta,b; int*p,*q; p=&x;q=&y;z=0; printf("请输入两个数:");scanf("%d%d",&a,&b);if(extended_Gcd(a,b,*p,*q)==1) { extended_Ivn(a,b,&z);printf("%d和%d互素，乘法的逆元是：%d\n",a,b,z); }else { z=extended_

人人文库> 全部分类> 教育资料 > 中学教育

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论