1、RSA 算法的实现摘 要本文设计的是一套完整实用的 RSA 文件加密解决方案，并具体编码实现。本文采用费马小定理测试素数，使用 Montgomery 加快大数模乘运算，用 C+实现 RSA 加密算法类库，并在 32 位 windows 平台封装成组件。在.Net 平台引用此组件，实现可以对任意文件进行 RSA 加密操作的窗体应用程序。经过加密的文件以及密钥文件都是文本文件。本文首先给出关键类类图、整个应用程序的结构描述文档，然后对关键模块流程图、详细的接口文档进行阐述，并给出关键的实现代码，最后对应用程序进行测试，对测试结果进行分析研究，进而对应用程序进行改进，对关键算法进行尽可能的优化，最终
2、得到一个在 windows运行的可以用指定密钥对任意文件进行 RSA 加密并可解密的完整应用程序，和一些相关的可移植组件。关键词： RSA；文件加密；Montgomery；费马定理Implement of RSA AlgorithmAbstractIn this paper, a solution of encrypting file with RSA algorithm and the codes of this system are introduced. Fermat theory is used to test prime number. Montgomery is used to c
3、ut short the time of modular multiplication of large number. The class library of RSA is implemented in C+, and packaged to component on the platform of 32 bits windows. On the platform of .Net, the application is implemented with reference of this component and can encrypt any file with RSA. Both e
4、ncrypted files and key files are text files. In this paper, core class figures and the framework are first introduced. Then the flow of core modules and detail interfaces are stated and the kernel codes are showed also. Finally, it analyzes the result of test, then optimizes core algorithm. In the c
5、onclusion, an entire application which can encrypt any files with RSA algorithm using given key and some transplanted components are implemented.Key words: RSA ; File Encryption ; Montgomery ; Fermat目 录论文总页数：35 页1 引言 11.1 课题背景 11.2 RSA 算法介绍与应用现状 11.3 RSA 应用于文件加密的分析 21.3.1 文件加密使用 RSA 的可行性 .21.3.2 文件加
6、密使用 RSA 的意义 .32 RSA 文件加密软件的设计与实现 .42.1 需求分析与总体设计 .42.1.1 功能分 析 .42.1.2 工程方案选择 .42.2 各部分的设计与开发 .52.2.1 实现 RSA 加密算法的 C+核心类库 52.2.2 封装 C+核心类库的 DLL 组件 .252.2.3 引用 DLL 的.Net 类与实现文件操作功能的窗体应用程序 263 软件整体测试与分析改进 .273.1 编写测试各项性能需要的精确计时类 273.2 测试数据与分析改进 273.2.1 密钥生成测试 273.2.2 数据输入输出测试 283.2.3 加密解密测试 29结 论 31参考
7、文献 32附 录 33致 谢 34声 明 35第 1 页 共 35 页1 引言1.1 课题背景RSA 公钥加密算法是第一个既能用于数据加密也能用于数字签名的算法。它易于理解和操作，也十分流行。算法的名字以发明者的姓氏首字母命名：Ron Rivest, Adi Shamir 和 Leonard Adleman。虽然自 1978 年提出以来，RSA 的安全性一直未能得到理论上的证明，但它经历了各种攻击，至今（2006 年）未被完全攻破。随着越来越多的商业应用和标准化工作，RSA 已经成为最具代表性的公钥加密技术。VISA、MasterCard、IBM、Microsoft 等公司协力制定的安全电子交
8、易标准（Secure Electronic Transactions，SET）就采用了标准 RSA算法，这使得 RSA 在我们的生活中几乎无处不在。网上交易加密连接、网上银行身份验证、各种信用卡使用的数字证书、智能移动电话和存储卡的验证功能芯片等，大多数使用 RSA 技术。当今公钥加密更广泛应用于互联网身份认证，本课题将公钥加密算法 RSA应用于小型文件加密。将任意文件加密成文本的解决方案，使其使用更加灵活。整个工程的分层设计，给引用移植和后续开发带来便利。1.2 RSA 算法介绍与应用现状RSA 算法可以简单叙述如下：取素数 p，q，令 n=pq.取与(p-1)(q-1)互素的整数 e，由方
9、程 de=1 (mod (p-1)(q-1)解出 d，二元组(e,n)作为公开密钥，二元组(d,n)作为私有密钥b=ae mod n， c=bd mod n.附录中给出了证明 a=c (mod n).RSA 公开密钥加密算法自 20 世纪 70 年代提出以来，已经得到了广泛认可和应用。发展至今，电子安全领域的各方面已经形成了较为完备的国际规范。RSA 作为最重要的公开密钥算法，在各领域的应用数不胜数。RSA 在硬件方面，以技术成熟的 IC 应用于各种消费类电子产品。RSA 在软件方面的应用，主要集中在 Internet 上。加密连接、数字签名和数字证书的核心算法广泛使用 RSA。日常应用中，有
10、比较著名的工具包 Open SSL(SSL，Security Socket Layer,是一个安全传输协议，在 Internet 上进行数据保护和身份确认。Open SSL 是一个开放源代码的实现了 SSL 及相关加密技术的软第 2 页 共 35 页件包，由加拿大的 Eric Yang 等发起编写的。相关详细介绍见http:/www.openssl.org/about/ )。Open SSL 应用 RSA 实现签名和密钥交换，已经在各种操作系统得到非常广泛的应用。另外，家喻户晓的 IE 浏览器，自然也实现了 SSL 协议，集成了使用 RSA 技术的加密功能，结合 MD5 和 SHA1，主要用于
11、数字证书和数字签名，对于习惯于使用网上购物和网上银行的用户来说，几乎天天都在使用 RSA 技术。1.3 RSA 应用于文件加密的分析1.3.1 文件加密使用 RSA 的可行性通过 1.2 节的论述，不难看出 RSA 当今的应用多在于数字签名和证书等方面。之所以只应用于这些短小数据的加密解密，是因为 RSA 算法加密极慢，速度是 DES 对称密钥加密速度的千分之一左右。正是因为这样，把 RSA 应用于普通文件加密的想法一直被忽略。通常文件被想象成大数据块，但是实际上在日常应用中，有些极其重要的文本资料是并不太大的，比如因担心遗忘而用普通文本记录的银行帐号和密码、不应被陌生人知道的重要电话号码、几
12、千字节大的重要小图片等。虽然 RSA 加密运算的速度十分慢，但是在 PC 性能越来越好的今天，对于几千字节的数据进行一次几百位密钥的 RSA 加密，所消耗的时间应该是可以接受的。下面结合大数运算程序的调试，从理论上简单的分析消耗时间。在一台普通配置的 PC 机上对一个整数进行幂模运算，因为公开密钥的 e 通常取的较小，所以指数取一个小整数，比如 C353，模一个 70 字节长的整数(140 位十六进制，大数单元以线性组方式实现，对应到 RSA 算法中，这相当于约 560bit的 n)，调试一个函数测试，按初等数论中的知识对程序进行算法优化，最终在一台配置为 AMD Athron2800+，外频
13、 333MHZ，物理内存 512MB 的 PC 上测试需要约 45 毫秒时间。如果按这种速度，逐字节对 1KB 的数据进行同样的运算，所消耗的时间理论上为 45 毫秒的 1024 倍即约 45 秒。这个时间并不是非常长。其实从一个简单的角度来说，既然 RSA 用于数字签名可行，那就完全可以用于同样大小的普通文件。对于较大的文件，如果分成与数字签名同样大小的段(这里假设数字签名较短，不分段一次计算加密完成)，分开的各段逐一进行加密运算，那所需要的时间也只是按文件大小线性的增长。通常数字签名为几十字节，加密运算并不需要很长的等待，这就说明对于几百字节或一两 K 字节大小的文件来说，如果进行 RSA
14、 加密，并不会是非常漫长的工作。当然，如果文件更大，加密就显得十分漫长了。比如按前面叙述的 45 毫秒大数运算程序推理，加密 1M 字节大小的文件需要约 1 天的时间。所以，要在普通 PC 用几百位以上的长密钥 RSA 加密文件，文件不能过大，一般可以接受的上限是几 KB。第 3 页 共 35 页如果要在较短时间内加密大文件，需要缩短密钥长度以减小运算量，这将带来安全性隐患。本文的第 3 章将根据实际调试好的软件，测试给出具体的时间消耗数据。例如，在一台配置为 AMD Athron2800+，外频 333MHZ，物理内存 512MB 的PC 上测试实现的软件，以 560bit 的 n 逐字节加
15、密一个 1KB 大小的文件需要 55秒。通常记录如银行帐号密码等重要数据的文本文件大小不足百字节，加密只需要数秒钟。所以对于小型文件，进行较长密钥的 RSA 加密是完全可行的。1.3.2 文件加密使用 RSA 的意义如 1.3.1 节所述，小型文件加密可以使用 RSA。比如，因担心遗忘而用普通文本记录的银行帐号和密码、不应被陌生人知道的重要电话号码、几千字节大的重要小图片等。可行的方法未必是必要的，本小节讨论何种文件适合用非对称密钥加密，即 RSA 加密文件的意义所在。对于前面叙述的带有重要信息的小型文本和二进制数据的维护，如果不加密，将无法放心的保存在计算机上，尤其是连网的或机房里的公共计算
16、机。如果借助功能强大的大型多用户数据保护程序维护几个小型文件，显得十分烦琐，好比杀鸡用牛刀。如果采用对称密钥加密，即加密解密的密钥相同，只适合部分情况。在某些情况下，使用对称密钥加密文件，交流使用不够方便。比如，张三由于某种原因，需要将自己的某个文件在公共计算机上留给李四，而不希望别人看到内容。如果采用对称密钥加密，张三和李四提前约好一个密码就可以。但是如果张三想要在同一台公共计算机上再留一个秘密文件给王五，而不希望别人看到，就要和王五另外约定一个密码。如果需要在这台公共计算机上留十个文件给不同的人，自己就要记和十个人约定好的密码，这样以来交流起来不够方便，因为对于张三，要自己维护太多的密钥。
17、非对称密钥(公开密钥方式)恰好解决这样的问题。只要大家都在这台计算机或这台计算机可以访问到的地方，留下自己的公开密钥，一切就变的容易解决了。张三要留给李四的文件，就用李四的公开密钥加密，要留给王五的文件，就用王五的公开密钥加密。李四和王五只要把留给自己的文件用自己的私有密钥解密，就可以得到留给自己的文件了。显然，非对称密钥体制更适合多用户交流，而将这种加密方式直接应用于文件加密，使我们在公开场合的交流更加灵活方便。综上所述，使用前面叙述的方式加密文件有两点重要意义：应用非对称密钥加密任意文件，使非对称密钥的应用不仅仅局限于互联网络。非对称加密后的数据变换成文本，使得我们可以通过几乎任何方式安全
18、传递任意文件，比如在只有 http 的环境使用 xml 方式。第 4 页 共 35 页2 RSA 文件加密软件的设计与实现2.1 需求分析与总体设计2.1.1 功能分析经过 1.3.2 节的论述，我们可以将对软件的要求总结如下： 可以按要求的位数生成非对称密钥。 可以用指定密钥以 RSA 算法加密任意一个文件，加密生成的数据为纯文本。 可以装载加密过的文件，并用指定的密钥解密还原出原文件。 提示信息完整、操作舒适、图形界面雅观按上述描述，给出 Use Case 和 Statechart 如图 2-1。图 2-1 本项目的 Use Case 和 Statechart根据以上分析，一般来说，需要进
19、行编码的程序有 RSA 密钥生成 RSA 加密解密 任意文件的读取 各环节必要的数据编码转换 图形操作界面。2.1.2 工程方案选择综合考虑复用性、可维护性和执行效率，较妥当的方法是分层设计。核心的 RSA 算法由 C+类库实现，针对用户所在的操作系统封装成本地化组件。其他各功能如文件操作、数据编码转换和图形界面等，由托管代码借助虚拟机平台标准库的功能快速开发实现(本文针对选用.Net 上的 C#论述，选用 java 由 JNI或其他方式调用本地组件，设计模式上是完全类似的)。这种开发方式，核心功能集中在最底层，在不断的封装中针对具体环境对组件功能不断扩充，任意一个层面的封装都可以被直接应用到
20、其他项目，比如在 Web 使用以前为某窗体程第 5 页 共 35 页序写的组件、给嵌入式设备交叉编译算法库等。但是每一层都需要依赖底层的所有组件。图 2-2 形象的说明了分层设计给复用带来的好处。图 2-2 综合考虑复用性、可维护性和执行效率的分层设计选用这种设计方案，上层使用 C#，底层算法使用 C+，可以由一个 Visual Studio 解决方案管理，给调试带来极大的方便。整个工程分四层，实现 RSA 加密算法的 C+核心类库、封装 C+核心类库的 DLL 组件、引用 DLL 的.Net 类、实现文件操作功能的.Net 窗体应用程序。2.2 节详细介绍各部分的设计与开发。考虑到工作量，本
21、软件加解密数据没有严格遵从 RSA 标准 PKCS #1，而是在满足设计要求的前提下，以一种尽可能简单的方式实现加密和解密。2.2 各部分的设计与开发2.2.1 实现 RSA 加密算法的 C+核心类库1. 大数存储和四则运算根据 RSA 算法的要求，为了实现大数的各种复杂运算，需要首先实现大数存储和基本四则运算的功能。当今开源的大数运算 C+类有很多，多用于数学分析、天文计算等，本文选用了一个流行的大数类型，并针对 RSA 算法和本项目的具体需要对其进行了扩充和改进。下面简单介绍大数存储和四则运算的实现原理。最先完成的功能是大数的存储，存储功能由 flex_unit 类提供。和普通的类型一样，
22、每一个大数对应一个 flex_unit 的实例。类 flex_unit 中，用一个无符号整数指针 unsigned * a 指向一块内存空间的首地址，这块内存空间用来存储一个大数，所以可以说，大数是被存储在一个以 unsigned 为单元的线性组中。在方法 void reserve( unsigned x )中通过 C+的 new 来给 a 开辟空间，当 flex_unit 的第 6 页 共 35 页实例中被存入比当前存储的数更大的数时，就会调用 reserve 来增加存储空间，但是当 flex_unit 的实例中被存入比当前存储的数更小的数时，存储空间并不会自动紧缩，这是为了在运算的时候提高
23、执行效率。结合指针 a，有两个重要的无符号整数来控制存储，unsigned z 和 unsigned n，z 是被分配空间的单元数，随数字变大不断增大，不会自己紧缩，而 n 是当前存储的大数所占的单元数，组成一个大数的各 unsigned 单元的存入和读出由 set、get 方法完成，变量 n 是只读的。类型 unsigned 在 32 位机是 32 位的，所以对于 flex_unit 这个大数类来说，每个大数最大可以达到 2*32 个字节长，这已经超过了 32 位机通常的最大内存容量，所以是足够进行 RSA 所需要的各种运算的。图 2-3 形象的说明了大数存储类 flex_unit 对大数的
24、管理。*a unsigned 类型的指针大数占 n 个单元开辟了 z 个单元大的内存内存空间图 2-3 flex_unit 对大数的管理在 flex_unit 的存储功能基础上，将其派生，得到 vlong_value，在vlong_value 中实现四则运算函数，并实现强制转换运算符 unsigned，以方便大数类型和普通整数的互相赋值。当大数被强制转换为 unsigned 时，将取其最低四字节的值。四则运算实现的原理十分简单，都是按最基本的算术原理实现的，四则运算过程的本质就是按一定数制对数字的计算，比如相加，就是低位单元对齐，逐单元相加并进位，减法同理。而乘除法和取余也都是按照竖式运算的原
25、理实现，并进行了必要的优化。虽然实现了四则运算函数，但是若是程序里的运算都要调用函数，显得烦琐而且看起来不美观，所以我们另写一个类vlong，关联(Associate，即使用 vlong_value 类型的对象或其指针作为成员)vlong_value，在 vlong 重载运算符。这样，当我们操作 vlong 大数对象的时候，就可以像使用一个简单类型一样使用各种运算符号了。之所以将 vlong_value的指针作为成员而不是直接构造的对象，也是为了提高执行效率，因为大型对象的拷贝要消耗不少机器时间。第 7 页 共 35 页2. 大数幂模与乘模运算Montgomery 幂模算法在实现了 vlong
26、 类型后，大数的存储和四则运算的功能都完成了。考虑到RSA 算法需要进行幂模运算，需要准备实现这些运算的方法。所以写一个 vlong的友元，完成幂模运算功能。幂模运算是 RSA 算法中比重最大的计算，最直接地决定了 RSA 算法的性能，针对快速幂模运算这一课题，西方现代数学家提出了很多的解决方案。经查阅相关数学著作，发现通常都是依据乘模的性质，先将幂模运算化简为乘模运算。nbnabamod)()od(m)(通常的分解习惯是指数不断的对半分，如果指数是奇数，就先减去一变成偶数，然后再对半分，例如求 D= ，E=15 ，可分解为如下 6 个乘模运算。CEnCod21m32623 nod74CC14
27、4556归纳分析以上方法，对于任意指数 E，可采用如图 2-4 的算法流程计算 。第 8 页 共 35 页开始D=1;P=C mod nE0?E 为奇数? nPDmod)(E=E-1 nPmod)(E 为偶数?E=E/2YesNoResult=D结束YesYesNoNo图 2-4 幂模运算分解为乘模运算的一种流程按照上述流程，列举两个简单的幂模运算实例来形象的说明这种方法。 求 的值17mod25开始 D = 1 P = 2 mod 17 = 2 E = 15E 奇数 D = DP mod n = 2 P = PP mod n = 4 E= (E-1)/2 =7E 奇数 D = DP mod
28、n = 8 P = PP mod n = 16 E= (E-1)/2 =3第 9 页 共 35 页E 奇数 D = DP mod n = 9 P = PP mod n = 1 E= (E-1)/2 =1E 奇数 D = DP mod n = 9 P = PP mod n = 1 E= (E-1)/2 =0最终 D = 9 即为所求。 求 的值13mod28开始 D = 1 P = 2 mod 13 = 2 E = 8E 偶数 D = 1 P = PP mod n = 4 E = E/2 =4E 偶数 D = 1 P = PP mod n = 3 E = E/2 =2E 偶数 D = 1 P =
29、 PP mod n = 9 E = E/2 =1E 奇数 D = DP mod n = 9 P = 不需要计算 E = (E-1)/2 =0最终 D = 9 即为所求。观察上述算法，发现 E 根据奇偶除以二或减一除以二实际就是二进制的移位操作，所以要知道需要如何乘模变量，并不需要反复对 E 进行除以二或减一除以二的操作，只需要验证 E 的二进制各位是 0 还是 1 就可以了。同样是计算 ，下面给出从右到左扫描二进制位进行的幂模算法描述，设中nCEmod间变量 D,P,E 的二进制各位下标从左到右为 u,u-1,u-2,0。Powmod(C,E,n) D=1;P=C mod n;for i=0
30、to u do if(Ei=1)D=D*P(mod n); P=P*P(mod n);return D;有些文献将上述算法称为平方乘积二进制快速算法，例如参考文献中的基于 RSA 算法的一种新的加密核设计 ，其实这种算法本质上和图 2-4 的流程完全一致，只是把根据指数奇偶分开的减一和除以二合并成对指数二进制各位的判断而已。在本软件的代码中采用直接扫描 vlong 二进制各位的办法。剩下的问题就是乘模运算了。提高乘模运算的速度是提高模幂运算速度的关键。一般情况下，n 是数百位乃至千位以上的二进制整数，用普通的除法求模而进行乘模运算是不能满足速度的要求的。为此，Montgomery 在 1983
31、 年提出了一种模加右移的乘模算法（主要著作发表于 1985 年） ，从而避免了通常求模算法中费时的除法步骤。本软件仅仅是应用 Montgomery(蒙哥马利)算法，算第 10 页 共 35 页法的具体推导证明需要颇多数论知识，不在本文的讨论范围内，如需了解可参见蒙哥马利的相关著作。下面简单描述 RSA 中常用的 Montgomery(蒙哥马利)算法供参考理解源程序。选择与模数 n 互素的基数 R=2k，n 满足 2k1 n=n) x -= n;return x;exp(C,E,n) /蒙哥马利幂模D=R-n;P=C*R mod n;i=0;第 11 页 共 35 页while(true)if(
32、E 的当前二进制位 Ei=1)D=M(D*P); /从低位到高位检测二进制位i+=1;if(i=E 的二进制位数)break;P=M(P*P);return D*R-1 (mod n);在具体的实现中，对应 monty 类的 mul 和 exp 方法。全局函数 modexp 初始化 monty 对象并调用其 exp 方法，使用的时候直接调用 modexp 即可。3. 寻找素数Eratosthenes 筛选与 Fermat 素数测试首先要说明的是，事实上，当今的计算机还不足以聪明到立刻计算生成一个很大的随机素数。一般来说，要得到 100%准确的大素数，都是通过查已经计算好的素数表的方式。但是素数
33、表的方式给 RSA 的安全性带来隐患，因为攻击者如果得到了密钥生成时所使用的素数表，攻破 RSA 加密的难度将会大大降低。本程序起初使用素数表的方式，后来考虑到安全性问题，生成密钥的方式改为随机计算生成。这样，短时间内如果要得到一个 100%准确的大素数是很困难的，只能以尽可能高的概率得到一个大素数。经过 2.2.1.1 和 2.2.1.2 小节，所有的大数运算功能都准备完毕，在此基础上，本工程将寻找素数的功能置于类 Prime_factory_san 之中。外部只要调用本类实例的成员 vlong find_prime( vlong inp;i+)unsigned p = pli;unsign
34、ed r = start % vlong(p);if (r) r = p - r;while ( r SS )br = 0;r += p;这里利用 start 对各小素数因子 p 求模的办法，得到当前 p 在素数搜索范围内的最小倍数在 b中的对应位置，将其剔除后，不断后移 p 个位置，将这个小素数因子 p 在搜索范围内的所有倍数全部剔除，如图 2-5 所示。在完成对所有小素数因子的类似操作后，他们的倍数在搜索范围内的位置标记 br被全部标记为 0。实际上这就是 Eratosthenes 筛选法。图 2-5 在素数搜索范围内剔除小素数因子 p 的倍数接下来，对可能为素数的数(即标记数组 b中值为
35、 1 的元素对应的数)进行素数测试。数论学家利用费马小定理研究出了多种素数测试方法，本程序使用一种最简单的方式，直接应用费马小定理。取一个与 p 互素的整数 A，对于大素数 p 来说应该满足 Ap-1mod p=1，但是我们把 p 代入一个大整数，满足这个关系的数不一定是素数。这时我们改变 A，进行多次测试，如果多次测试都通过，这个数是素数的概率就比较大。按这种原理，我们编写素数测试函数如下。int is_probable_prime_san( const vlong /测试次数const unsigned anyrep = 2,3,5,7 ; /测试用的底数for ( unsigned i=
36、0; irep; i+=1 )if ( modexp( anyi, p-vlong(1), p ) != vlong(1) ) return 0;/modexp 是幂模函数，按上一小节叙述的算法编码。/这里 modexp 计算 anyip-1mod p。return 1;测试通过，程序就判定这个数为找到的素数，将找到的素数返回给上层程序使用。在这里其实有一个不可忽视的问题，就是得到一个测试通过的合数。对于这种情况，RSA 算法加密解密是否还可以实现，是一个需要从数学角度论证的问题。因为得到素数的概率很高，经过一整天的生成密钥和加密操作，没有发现失败的密钥， 所以本文暂没有对这个问题进行讨论。实
37、际得到素数的流程：(1) 先得到一个随机的大整数 N 当作寻找的起点(2) 确定一个寻找范围的大小 SS，把(N,N+SS)范围内的小素数倍数去掉，即前面叙述的古希腊某人发明的筛选法小素数因子从开始取，取几百个( 论文中将小素数因子个数记为 NP)(3) 对范围内没有去掉的数逐一进行素数测试，一个数如果通过测试次数达到一定标准，就判为素数(4) 如果范围内没找到素数，就令 N=N+SS，回到(2)继续寻找用以上算法，直到以某成功概率得到素数为止综上所述，总结素数寻找的流程，如图 2-6 所示。第 14 页 共 35 页开始按 start 参数初始化标记数组 bSS ; i=0计数 iSS?bi
38、=1?判定为素数?start+=1;i+=1结束返回素数寻找结果 startYesYesYesNoNoNo图 2-6 函数 find_prime 寻找素数的流程框图得到了大素数，即 RSA 算法中的 p、q，我们就可以计算出密钥，进行加密等操作了。4. 二元一次不定方程在 RSA 算法中，往往要在已知 A、M 的情况下，求 B 的最小值，使得 (AB) mod M = 1。即相当于求解 B、N 都是未知数的二元一次不定方程 AB-MN=1 的最小整数解。而针对不定方程 ax-by=1 的最小整数解，古今中外都进行过详尽的研究，西方有著名的欧几里德算法，即一种辗转相除法，中国有秦九韶的“大衍求一
39、术” 。欧几里德算法是一种递归算法，较容易理解。下面举例说明用欧几里德算法求解二元一次不定方程的最小整数解。给定不定方程 11x-49y=1，求最小的 x（1） 11 x - 49 y = 1 49 mod 11 = 5（2） 11 x - 5 y = 1 11 mod 5 = 1第 15 页 共 35 页（3） x - 5 y = 1 5 mod 1 = 0逆向代入：令 y=0 代入（ 3）得 x=1令 x=1 代入（ 2）得 y=2令 y=2 代入（ 1）得 x=9x=9;y=2 即为所求。程序中，全局函数 vlong modinv( const vlong Byteb=new ByteS
40、ystem.Math.Max(RSAprimeplen1,RSAprimeplen2);s1=“;s2=“;for(int i=0;iRSAprimeplen1;i+)Byte tmp=System.Convert.ToByte(254.0*rnd.NextDouble();if(tmp!=0)bi=tmp;else bi=1; s1=wujunjie_rsa.FromASCIIByteArray(b);for(int i=0;iRSAprimeplen2;i+)Byte tmp=System.Convert.ToByte(254.0*rnd.NextDouble();if(tmp!=0)bi
41、=tmp;else bi=1;s2=wujunjie_rsa.FromASCIIByteArray(b);return 1;第 17 页 共 35 页加密过程：图 2-8 载入待加密的文本图 2-9 准备加密文本第 18 页 共 35 页图 2-10 加密后生成的文本图 2-11 加密过程完成相关代码：private void menuItem10_Click(object sender, System.EventArgs e)/公钥加密if(wujunjie_rsa.charlist.Count=0)emptymsg em=new emptymsg(this); em.Show();retu
42、rn;第 19 页 共 35 页Stream myStream ;SaveFileDialog saveFileDialog1 = new SaveFileDialog();saveFileDialog1.Filter = “Hex text files (*.hextxt)|*.hextxt|All files (*.*)|*.*“ ;saveFileDialog1.FilterIndex = 1 ;saveFileDialog1.RestoreDirectory = true ;if(saveFileDialog1.ShowDialog() = DialogResult.OK)if(myS
43、tream=saveFileDialog1.OpenFile() != null)textBox1.Text+=“rn正在对读入的文件进行处理，请稍候：）rn“;System.Threading.Thread.Sleep(500);HighResolutionTimer timer = new HighResolutionTimer(); timer.Start();using (StreamWriter sw = new StreamWriter(myStream) sw.WriteLine(“# RSA.HexText“);sw.WriteLine(“#_“);Byte b=new Byt
44、ewujunjie_rsa.RSAstep;wujunjie_rsa.result_hexstrings.Clear();progressBar1.Minimum=0;progressBar1.Maximum=wujunjie_rsa.charlist.Count;for(int i=0;iwujunjie_rsa.charlist.Count;i+=System.Convert.ToInt32(wujunjie_rsa.RSAstep)for(int j=0;jwujunjie_rsa.RSAstep;j+)第 20 页 共 35 页bj=System.Convert.ToByte(wuju
45、njie_rsa.charlisti+j);string s;wujunjie_rsa.RSA_san_en(b,wujunjie_rsa.RSAstep);s=wujunjie_rsa.get_result_hexstring();wujunjie_rsa.result_hexstrings.Add(s);progressBar1.Value=i+1;for(int i=0;iwujunjie_rsa.result_hexstrings.Count;i+)string hs=System.Convert.ToString(wujunjie_rsa.result_hexstringsi);if
46、(hs=null|hs=“)sw.WriteLine(“0“);else sw.WriteLine(hs);sw.WriteLine(“#_“);sw.Write(“# “);sw.WriteLine(DateTime.Now);wujunjie_rsa.result_hexstrings.Clear();myStream.Close();timer.Stop(); textBox1.Text+=“rn消耗时间：“+timer.ElapsedTime+“rn“;textBox1.Text+=“rn处理完成，新生成文件rn“+saveFileDialog1.FileName+“rn“;第 21
47、页 共 35 页progressBar1.Value=0; 解密过程：图 2-12 载入加密后生成的文本图 2-13 解密已加密的文本第 22 页 共 35 页图 2-14 解密生成解密文件图 2-15 解密过程完成相关代码：private void menuItem8_Click(object sender, System.EventArgs e)/私钥解密if(wujunjie_rsa.hextxtlist.Count=0)emptymsg em=new emptymsg(this); em.Show();第 23 页 共 35 页return;Stream myStream ;SaveF
48、ileDialog saveFileDialog1 = new SaveFileDialog();saveFileDialog1.Filter = “Text files (*.txt)|*.txt|All files (*.*)|*.*“ ;saveFileDialog1.FilterIndex = 2 ;saveFileDialog1.RestoreDirectory = true ;if(saveFileDialog1.ShowDialog() = DialogResult.OK)if(myStream=saveFileDialog1.OpenFile() != null)textBox
49、1.Text+=“rn正在对十六进制文本进行处理，请稍候：）rn“;System.Threading.Thread.Sleep(500);HighResolutionTimer timer = new HighResolutionTimer(); timer.Start(); using (BinaryWriter bn = new BinaryWriter(myStream)progressBar1.Minimum=0;progressBar1.Maximum=wujunjie_rsa.hextxtlist.Count;for(int i=0;iwujunjie_rsa.hextxtlist.Count
