实验三RSA算法代码的实现.docVIP

实验三 RSA算法代码实现 代码实现结果: 四、实验要求 1、在RSA算法中，求解私钥d可以采用扩展欧几里得算法和线性代入法两种方法，但这两种方法的计算效率很低（即计算d要花费很多时间，特别是数次乘积运算占用了大量计算时间，使得RSA时间复杂度增加），请你通过网络或自己设计一种及一种以上更为高效的求解d的方法，并写出详细计算过程。 2、用VC调试实验步骤给出的RSA加解密程序，结合原理分析程序中密钥是如何产生的？如何加密的？如何解密的？并在程序中给于解释。 3、你还能采用别的方法自主设计RSA公钥密码算法吗？如果能请列出相关代码。 【参考答案】 在RSA算法中，求解私钥d可以采用以下方法： 扩展欧几里得算法(模幂乘法运算) 线性代入法(模幂乘法运算) 平方-乘法：按平方计算模指数运算 一种二进制查表方幂模快速计算方法：该方法根据指数二进制形式结构进行分组计算，记忆典型操作数，通过查表法大量减少了乘法次数。 采用分解计算和内联汇编方法加速RSA算法的执行 设计规划密钥空间 【代码分析】 #includeiostream #includevector #includestring using namespace std; int Euclid(int a, int n)//na,求n和a的最大公因数为1 { int x,y,r; x=n;y=a; for (int i=0;;) { if(y==0) return x;//a=0 if(y==1) return y;//a=1 r=x%y;//a≠0,1,r=n%a x=y; y=r; } } double extenEuclid(double a,double n)//利用扩展的Euclid计算a^-1 mod n 的乘法逆元 { double x1=1,x2=0,x3=n,y1=0,y2=1,y3=a,Q; double t1,t2,t3; for(int i=0;;) { if(y3==0) { return x3; coutno reverseendl;//a=0,无乘法逆元 } if(y3==1) return y2;//a=1,a%n=1 Q=int(x3/y3);//Q=[n/a] t1=x1-Q*y1;//t1=x1=1 t2=x2-Q*y2;//t2=Q t3=x3-Q*y3;//t3=n-Qa,相当于d原理中的1=de-xfn x1=y1;x2=y2;x3=y3; y1=t1;y2=t2;y3=t3;//返回d的结果 } } bool Rabin(int a,int n)//Miler-Rabin素性测试算法对一个给定的大数进行测试 { vectorint b;//判断一个变量b是否是素数 unsigned int N=n-1; for(int i=0,j=1;;i++) { if(jN) break; if((Ni)(unsigned int)1) b.push_back(1); else b.push_back(0); j*=2; }//将n-1表示成二进制形式 for(int k=0;kb.size();k++) coutb[k]; coutendl; int d=1,x=0; for(int ii=b.size()-1;ii=0;ii--) { x=d; d=(d*d)%n; if(d==1x!=1x!=n-1) return false; if(b[ii]==1) d=(d*a)%n; }//b为素数 if(d!=1) return false;//不为素数 return true; } double quickindex1(double a,double m, double n)//实现a^m mod n 的运算 { vectorint b; unsigned double N=m;//公钥m=e for(int ii=0,j=1;;ii++) { if(jN) break; if((Nii)(unsigned double)1) b.push_back(1); else b.push_back(0); j*=2; }//将m表示成二进制形式 double c=0,d=1; for(int i=b.size()-1;i=0;i--) { c*=2; d=(d*d)-int((d*d)/n)*n; if(b[i]==1) { c+=1; d=(d*a)-int((d*a)