版权申明:本文为博主窗户(Colin Cai)原创,欢迎转帖。如要转贴,必须注明原文网址

  http://www.cnblogs.com/Colin-Cai/p/7290010.html 

  作者:窗户

  QQ:6679072

  E-mail:6679072@qq.com

  RSA最终加密、解密都要用到模乘的幂运算,简称模幂运算。

  回忆一下RSA,从明文A到密文B

  B=Ae1%N

  对B解密回到明文A,就是

  A=Be2%N

  其中,一般来说,加密公钥中的e1一般会比较小,取65537居多,但解密的时候,这个e2是一个非常非常大的数,显然,直接通过e2次模乘来解密是不现实的。

  为了让RSA的加密、解密成为现实,我们必须要找一个好的算法来做模幂运算。

  借上一节我设定的符号,以区别于传统上的幂的数学表示,

  定义a#b为a和b的模乘,

  定义a##n为n个a的模乘,或称a的n阶模乘。

  a = a##1,

  a##2 = a#a,

  a##4 = (a##2) # (a##2),

  ...

  a##2n = (a##2n-1) # (a##2n-1),

  从而,进行了n次模乘,分别获得了1、21 、22 ...2n 个a的模乘结果,

  而1、21 、22 ...2n 取其中任意一部分的和可以覆盖1~2n-1的所有整数。

  从而我们得到了一个算法来计算a##b:

  先不断通过每次得到的新结果自己同自己模乘,得到

  a、a##2、a##4...a##2n,使得满足2n≤b<2n+1

  再将b化为二进制表示,实际上也就是表示为2的各次幂相加形式,

  然后找到对应每个2的幂次a模乘结果,

  然后再把这些结果依次模乘,得到最终结果。

  打个比方,如果想求a##21,

  21用二进制表示为10101,所以21 = 16+4+1,

  第一步,求得a##2,a##4,a##8,a##16,

  第二步,取a,a##4,a##16,三者模乘就是最终需要的结果。

  

  当然,可以把第一步求a的2n阶模乘和第二步取需要的模乘融合在一起,这样就不需要存储每一个a的2n阶模乘结果,从而存储空间可以为常数级,而之前存储空间为线性级。

  该算法用bc描述如下:

  

#!/usr/bin/bc -q
define mod_mul(a1,a2,n)// a1*a2%n
{
return a1*a2%n;
}
define mod_exp(a,b,n)// a^b%n
{
while(b%==) {
a = mod_mul(a,a,n);
b /= ;
}
ret = a;
b /= ;
while(b!=) {
a = mod_mul(a,a,n);
if(b% == )
ret = mod_mul(a,ret,n);
b /= ;
}
return ret;
}
a=read();
b=read();
n=read();
print mod_exp(a,b,n),"\n"
quit

  该算法求a##b所做模乘次数:

  求各个a的2n阶模乘,所做模乘次数为log2b取整,也就是b的二进制的位数减1;

  取相应的2的正整数次幂的模乘结果再做模乘,所做模乘次数为b的二进制中1的个数减1。

  两者加一起为模乘次数。

  比如上面想求a##21,21用二进制表示为10101,所以该算法所需要模乘次数为(5-1)+(3-1)=6次。

  

  但此算法未必是最优的,我们来看看如下例子:

  想求a##441,

  441用二进制展开为110111001,于是模乘次数应该是(9-1)+(6-1)=13次。

  但我们可以用另外一种方式来求a##441,

  441 = 21 * 21,

  a##441 = (a##21) ## 21,

  我们求b = a##21需要6次模乘,

  再求b##21需要6次模乘,总共只需要12次模乘,比刚才少了1次。

  

  441是合数,我们再取个质数,239。

  239写成二进制是11101111,那么根据我们的算法应该做的模乘次数为(8-1)+(7-1)=13次。

  239 = 14*17+1,

  a##239 = ((a##14) ## 17) # a,

  先求b = a##14,需要5次模乘,

  再求c = b##1,需要5次模乘,

  最后再与a模乘,需要1次模乘,

  总共11次模乘,比13次要少。

  但是,即便是合数分解了,也未必得到更好的结果。

  例如49 = 7*7,

  分解前,49表示为二进制为110001,要做(6-1)+(3-1)=7次模乘,

  可是分解后,7表示为二进制制为111,

  总共要做(3-1)+(3-1)+(3-1)+(3-1)=8次模乘了。

  从而我们可以知道,我们给出的算法虽然是一个线性时间算法,但未必是最优算法,不过做到线性时间算法,从运算时间复杂度上来说,的确没有时间复杂度是这个算法的低阶无穷大级的算法了,从这个意义上,此算法已经“最优”了。

  如何做到最少次模乘?

  本问题为以下问题:

  (1)集合A初始为{1}

  (2)每一步从集合A中取两个数a和b,ab可相同,让c = a+b,再把c并入集合A,

    A = A∪{c}

  (3)输入正整数e,求A里面有元素e的最小步骤的过程。

  

  可惜此问题获得最优解似乎没有很好的算法,甚至远高于RSA可能基于的安全性——大数分解,但存在相对好的算法,从而可以用来改进我们的模幂算法。

  

  模幂算法是RSA的核心,不仅仅加密解密的时候需要,寻找密钥的时候也是需要的。

RSA简介(二)——模幂算法的更多相关文章

  1. RSA算法的C++string实现(模幂算法和欧几里得算法的使用)后附思路

    void resetNumA(string numAStr); //使用string重置numB void resetNumB(string numBStr); //将数组转换为字符串,用于输出 st ...

  2. RSA简介(四)——求逆算法

    此处所谓求逆运算,是指在模乘群里求逆. 第一节里提到互质的两个定义: (1)p,q两整数互质指p,q的最大公约数为1. (2)p.q两整数互质指存在整数a,b,使得ap+bq=1. 只要明白了欧几里得 ...

  3. RSA简介(三)——寻找质数

    要生成RSA的密钥,第一步就是要寻找质数,本节专讲如何寻找质数. 我们的质数(又称素数).合数一般是对正整数来讲,质数就是只有1和本身两个的正整数,合数至少有3个约数,而1既不是合数也不是质数. 质数 ...

  4. RSA算法原理——(2)RSA简介及基础数论知识

    上期为大家介绍了目前常见加密算法,相信阅读过的同学们对目前的加密算法也算是有了一个大概的了解.如果你对这些解密算法概念及特点还不是很清晰的话,昌昌非常推荐大家可以看看HTTPS的加密通信原理,因为HT ...

  5. Modular_exponentiation模幂运算

    https://en.wikipedia.org/wiki/Modular_exponentiation 蒙哥马利(Montgomery)幂模运算是快速计算a^b%k的一种算法,是RSA加密算法的核心 ...

  6. 二模13day1解题报告

    二模13day1解题报告 T1.发射站(station) N个发射站,每个发射站有高度hi,发射信号强度vi,每个发射站的信号只会被左和右第一个比他高的收到.现在求收到信号最强的发射站. 我用了时间复 ...

  7. 二模 (8) day2

    第一题: 题目描述: 有 n 个炸弹,有些炸弹牵了一根单向引线(也就是说引线只有在这一端能被炸弹点燃),只要引爆了这个炸弹,用引线连接的下一个炸弹也会爆炸.每个炸弹还有个得分,当这个炸弹被引爆后就能得 ...

  8. 非对称加密RSA、Elgamal、背包算法、Rabin、D-H、ECC(椭圆曲线加密算法)等。使用最广泛的是RSA算法

          非对称加密算法需要两个密钥:公开密钥(publickey)和私有密钥(privatekey).公开密钥与私有密钥是一对,如果用公开密钥对数据进行加密,只有用对应的私有密钥才能解密:如果用私 ...

  9. 银行卡号码校验算法(Luhn算法,又叫模10算法)

    有时候在网上办理一些业务时有些需要填写银行卡号码,当胡乱填写时会立即报错,但是并没有发现向后端发送请求,那么这个效果是怎么实现的呢. 对于银行卡号有一个校验算法,叫做Luhn算法. 一.银行卡号码的校 ...

随机推荐

  1. Linux系统vi模式下显示行号

    在命令模式下输入:set nu或者:set number都可以为vi设置行号,如果要取消的话,则输入:set nonu行号的设置是vi的环境设置,不会影响文本的内容.

  2. thinkphp中try catch的运用

    public function doedit_set(){ $info=$this->_post("info"); $id=$this->_post("id& ...

  3. JAVA容器结构图

  4. solr学习-基础环境搭建(一)

    目前网上关于solr6.+的安装教程很少,有些6.0之前的教程在应用到6.+的版本中出现很多的问题,所以特别整理出来这一片文章,希望能给各位码农一些帮助! 很少写些文章,如有不对的地方,还希望多多指导 ...

  5. 【Android Developers Training】 97. 序言:访问通讯录数据

    注:本文翻译自Google官方的Android Developers Training文档,译者技术一般,由于喜爱安卓而产生了翻译的念头,纯属个人兴趣爱好. 原文链接:http://developer ...

  6. 关于jquery全选反选 批量删除的一点心得

    废话不多说直接上代码: 下面是jsp页面的html代码: <table id="contentTable" class=""> <thead& ...

  7. CSS3学习系列之布局样式(一)

    多栏布局 使用float属性或position属性的缺点: 使用float属性或position属性进行页面中的简单布局的例子. <!DOCTYPE html> <html lang ...

  8. 把对象列表转化成json数据格式

    package JsonTest; import java.util.ArrayList; import java.util.List; public class test { public stat ...

  9. 2238"回文字串"报告

    题目描述: 回文串,就是从前往后和从后往前看都是一样的字符串.那么现在给你一个字符串,请你找出该字符串中,长度最大的一个回文子串. 输入描述: 有且仅有一个仅包含小写字母的字符串,保证其长度不超过50 ...

  10. C# Lambda表达式运用

    原文作者: C# Lambda表达式 原文作者2: Lambda表达式详解 Lambda表达式 "Lambda表达式"是一个匿名函数,是一种高效的类似于函数式编程的表达式,Lamb ...