密码学——加密演算法邓安文电子教案 第9章 质数与大整数算术.pptVIP

1.* ? 2006 第9章 质数与大整数算术 返回总目录 第9章 质数与大整数算术 教学目的 掌握大整数算术 了解Montgomery算术 了解Miller-Rabin质数测试 了解Agrawal-Kayal-Saxena算法 了解强质数和DSA质数 了解Java中的BigInteger Class套件 ? 大整数的加减乘法 本章内容 ? 大整数的除法 ? Montgomery算术 ? Miller-Rabin质数测试 ? Agrawal-Kayal-Saxena算法 ? 公开密钥密码的质数 ? Java的BigInteger Class ? 大整数算术与数论套件及软件 大整数的加减乘法 9.1 大整数的加减乘法 unsigned int carrier=0; for(unsigned short i=0;iL;i++) { unsigned int w=(unsigned int)x[i]+ (unsigned int)y[i]+carrier; if(w0xffff) { carrier=1; z[i]=(unsigned short)(w-0x10000); } else { carrier=0; z[i]=(unsigned short)w; }; }; z[L]=(unsigned short)carrier; 减法 算法 x-y=d unsigned int carrier=0; for(unsigned short i=0;iL;i++) { int w=(unsigned int)x[i]-(unsigned int)y[i] -carrier; if(w0) { carrier=1; d[i]=(unsigned short)(w+0x10000); } else { carrier=0; d[i]=(unsigned short)w; }; }; 乘法 算法 x×y=m unsigned int carrier=0; for (unsigned short i=0;i2*L-1;i++) { unsigned int w=carrier; for(unsigned short j=0;j=i;j++) w=w+(unsigned int)x[j]*(unsigned int)y[i-j]; if(w0xffff) { carrier=(unsigned short)w/0x10000; m[i]=(unsigned short)(w-carrier16); } //REM: carrier16=carrier*0x10000 else { carrier=0; m[i]=(unsigned short)w; }; }; m[2*L-1]=(unsigned short)carrier; Karatsuba乘法 算法Karatsuba乘法 KarastubaMultiply(x,y) { L=max(x.bitlength,y.bitlength); if(L1) { LL=ceil(L/2); x=x[1]*2?LL+x[0];//将x,y分割成2半 y=y[1]*2?LL+y[0]; A=KarastubaMultiply(x[0],y[0]) B=KarastubaMultiply(x[0]+x[1],y[0]+y[1]); C=KarastubaMultiply(x[1],y[1]); return(C,B-A-C,A);//C*4?LL+(B-A-C)*2?LL+A } else return x*y; //L=1 } 大整数的除法 9.2 大整数的除法 division(x,y) { x=(x[n],x[n-1],...,x[1],x[0]); y(y[n-1],y[n-2],...,y[1],y[0]); q_test=min(x[n]*(2?16)+x[n-1] / y[n-1],2?16-1); temp=q_test*y; while(temp=x) { q_test--; temp=temp-y; }; q=q_test; remainder=x-q*y; return(q,remainder) }; 除法 //主要算法部分: a=(a[n+m-1],...,a[1],a[0]); b=(b[n-1],...,b[1],