C++ GMP库实现任意位数的Montgomery Modular Multiplication

最新推荐文章于 2022-12-22 09:21:11 发布
最新推荐文章于 2022-12-22 09:21:11 发布
阅读量1.9k 收藏 3
点赞数 1
分类专栏： 软件安装 文章标签： gmp 蒙哥马利模乘
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_37333947/article/details/78329691
版权

嗯刚刚终于搞定了大数运算真是让我异常感动于是我又选择了来记录一下这次历程。

本来是想用C++实现一个MMM来测试verilog写出来的对不对的，刚开始觉得很容易啊不就把计算过程一写就好了，然而写到32位的时候求n'的时候需要用到r * r-1就发现溢出了……然后就不得不寻求大数运算的方法。


然而由于我太懒了不想自己动手写big integer，自己这方面的知识又很匮乏不知道有没有什么现成的工具，然后就搜了一下发现了GNU的GMP库然后就准备试一下了。


然而感觉网上很少配置GMP库的相关教程？或者也可能是我的打开方式不对，总之谷歌和度娘都搜不到太多有用的，只搜到了唯一一篇windows下配置GMP库的博客，还好良心博主的方法肥肠好用，终于不像之前装软件那样痛苦不堪。以下也就再写一下GMP库的安装方法，和GMP的C++ interface的使用（就用我的MMM代码啦）吧。


-------------------------------------------------------------------------------------以上是废话以下是正文的分割线------------------------------------------------------------------------------------------------


一、GMP库在Windows下的配置

参考文章：http://blog.csdn.net/u012629110/article/details/51220727

Emmm这篇文章里面已经讲得非常详细啦我就懒得再写一遍了（明明刚刚还想重新写的然而太懒了orz），配置过程中没出现什么大问题，就是本来昨天下午

最低0.47元/天 解锁文章
wenyq7
关注
  • 1
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
专栏目录
C++ GMP大数库头文件
05-02
C++大数库gmp头文件，大数运算分析的必备文件
C++ GMP库
wsawsdl1的博客
12-26 5081
简介 大整数运算库。 官网地址：https://gmplib.org/ 使用 初始化 mpz_t a, b, c;   mpz_init(a);   mpz_init_set_str(b, "2000000000000000000000000000", 10);   mpz_init_set_str(c, "3000000000000000000000000000", 10);  输出 ...
C++如何安装gmp库
weixin_35749796的博客
12-22 1685
要在 C 语言中安装 gmp 库，您需要执行以下步骤： 下载 gmp 库的源代码。您可以从 gmp 库的官方网站(https://gmplib.org/)下载最新版本的源代码。 解压下载的压缩文件。您可以使用任何文件解压工具，如 WinRAR 或 7-Zip，将压缩文件解压到您希望安装 gmp 库的目录。 进入 gmp 库的源代码目录。使用命令行工具(如 cmd 或 bash)进入到 gmp...
编写程序实现逆波兰运算_Montgomery算法介绍和实现
weixin_39744240的博客
11-22 751
小白个人杂记录，点上方蓝字即可关注~1介绍本次主要介绍Montgomery算法原理，并给出部分算法的具体实现。2Montgomery algorithmMontgomery algorithm，蒙哥马利算法，1985年由美国数学Peter.L.Montgomery提出，是一种快速执行模运算的方法，主要包括模约减、模乘和模指算法。2.1经典模乘算法计算x*y mod m，一般通用的方法是先...
New Hardware Architectures for Montgomery Modular Multiplication Algorithm
06-13
在此期间，Montgomery提出的蒙哥马利模乘算法因其在模指数运算中的基础作用而备受关注。模指数运算是RSA公钥密码系统的核心，同时也用于椭圆曲线密码系统及其他多种方法中。 #### 主要贡献 本文的主要贡献在于提出...
Modular Multiplication Without Trial Division.pdf
06-15
Modular Multiplication Without Trial Division。By Peter L. Montgomery。蒙哥马利模乘的论文。
Analyzing and Comparing Montgomery Multiplication Algorithms
06-13
“Analyzing and Comparing Montgomery Multiplication Algorithms”（分析与比较蒙哥马利模乘算法）这篇文章旨在深入探讨并对比不同的蒙哥马利模乘算法实现方法。蒙哥马利模乘算法是一种在计算机科学中广泛应用的...
gmp_montgomery
02-10
КурсоваяработаПАЗИ1（蒙哥马利曲线GMP_lib） Инструкцияпосборке： Еслибиблиотеканеустановлена（m4）（возможнопотребуетс...
tommath.zip C++大整数库
12-08
本文将深入探讨tommath库的功能、应用及其在C++中的实现。 tommath库，通常以libtommath的形式出现，是一个轻量级且功能强大的C语言库，主要设计用于处理超过标准整型范围的大整数。它提供了一系列全面的API（应用...
蒙哥马利算法VS2010源码
07-09
写了个蒙哥马利快速幂模算的源码，VS2010的源码，有注释 可以借此来理解这个算法的内容
GMP大数库实现大整数模以及Miller Rabin素数测试算法
12-15
GMP大数库的中文使用手册，以及已经编译好的GMP大数库，仅适用于VC6.0，并有自己写的生成随机大素数，大整数模运算，以及Miller Rabin素数测试算法。
Modular multiplication of polynomials
笨鸟
07-18 927
#include #include using namespace std; int pd(int sum[],int ls,int h[],int lh)//判断sum是否大于h { if(ls>lh) return 1; if(ls<lh) return -1; if(ls==lh) { int i; for(i=ls-1;i>=0;i--) { if(su
poj 1060 Modular multiplication of polynomials
07-17 973
Modular multiplication of polynomials Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 4032   Accepted: 1842 Description Consider polynomials whose coeffici
GMP学习：初始化mpz数组
suki570的博客
02-14 2291
GMP库中有函数 mpz_array_init ，但是注释了 This is an obsolete function. Do not use it. mpz_array_init 的问题在于它永远不会释放分配的内存。 看到网上没有什么有关的资料，因此记录一下使用 GMP 初始化数组的方法： #include <stdio.h> #include <gmp.h> int...
GMP手册笔记
chenxizhan1995的博客
11-20 9181
GMP 6.1.2 使用GMP时，不要使用本文档未说明的函数、宏、数据类型，否则保证你的程序不与未来版本的GMP兼容。 基础概念 头文件和库 使用GMP所需要的一切声明都汇集到了头文件gmp.h里，此头文件同时适用于C和C++。GMP中有些函数接受FILE*参数，只有当同时包含了stdio.h时，这些函数才是可用的。 #include <stdio.h> #include <g...
MbedTLS中的Montgomery算法实现解析（一）
qq_43324702的博客
11-11 1227
MbedTLS中的Montgomery算法实现解析（一） 蒙哥马利算法模乘介绍 蒙哥马利模乘算法主要是为了进行大数运算a×b mod n，并且针对CPU进行了并行优化，来快速计算模乘。 在介绍蒙哥马利模乘之前，先让我们来了解蒙哥马利约减（Montgomery reduction） 蒙哥马利约减（Montgomery reduction） 首先，让我们来考虑一个简单的问题： a mod n 如果a是一个2048位的整数，n是一个1024位的整数，那我们要如何来计算结果呢？ 如果我们直接采用相除的方式，不论在时
MbedTLS中的Montgomery算法实现解析（三）
qq_43324702的博客
11-13 1028
MbedTLS中的Montgomery算法实现解析（三） 在前面的两篇博客中，我们详细讨论了Montgomery算法的数学原理以及MbedTLS中是如何快速求得-N-1(mod n)。在这一篇博客中，我们将探讨MbedTLS中的蒙哥马利约减（Montgomery reduction）和蒙哥马利模乘（Montgomery multiplication）是如何实现的。 在阅读本博客之前，建议先掌握MbedTLS的大数表示方式 蒙哥马利约减Montgomery reduction 在我们的数学推导中，蒙哥马利约减
快速乘 + Montgomery modular multiplication
Scar_Halo
10-05 827
O（logn）:O（logn）:O（logn）: ll mult(ll a, ll b, ll p) { a %= p; b %= p; ll ans = 0; while(b) { if( b&amp;amp;amp;amp;amp;amp;amp; 1) ans = (ans + a) % p; a = (a + a) % p; b &amp;amp;amp;amp;amp;amp;gt;&amp;amp;
使用openssl库实现RSA 对比Montgomery模幂和Binary模幂
最新发布
05-05
RSA是一种非对称加密算法，公开密钥可以用于加密，私有密钥可以用于解密。在RSA算法中，公开密钥是由两个大素数的乘积N和另一个整数e组成的。私有密钥是由一个整数d组成的。RSA算法的加密和解密过程都涉及到了大数模幂运算。在实现RSA算法时，Montgomery模幂和Binary模幂是两种常用的算法。 Montgomery模幂算法是一种快速模幂算法，它通过将模数转换为另一个数来加速模幂运算。在Montgomery模幂算法中，模数是一个2的幂次。该算法需要预处理一些参数，但一旦预处理完成，每次模幂运算的时间复杂度都是O(log n)级别的，其中n是模数的位数。 Binary模幂算法是另一种常用的模幂算法，它是一种递归算法。在Binary模幂算法中，指数被表示为二进制数，然后从低位到高位进行幂运算。该算法的时间复杂度也是O(log n)级别的，但是它不需要预处理任何参数。 下面是使用openssl库实现RSA算法的示例代码，其中包括Montgomery模幂和Binary模幂两种算法的实现： ``` #include <openssl/rsa.h> #include <openssl/pem.h> // 使用Montgomery模幂实现RSA加密 int rsa_encrypt_montgomery(const unsigned char *in, int inlen, unsigned char *out, RSA *rsa, int padding) { int rsa_len = RSA_size(rsa); unsigned char *buf = (unsigned char *)malloc(rsa_len); // 将模数转换成Montgomery形式 BN_CTX *ctx = BN_CTX_new(); BN_MONT_CTX *mont_ctx = BN_MONT_CTX_new(); BN_MONT_CTX_set(mont_ctx, rsa->n, ctx); // 对输入数据进行加密 int ret = RSA_public_encrypt(inlen, in, buf, rsa, padding); // 对加密结果进行Montgomery还原 BN_from_montgomery(buf, buf, mont_ctx, ctx); // 将结果拷贝到输出缓冲区 memcpy(out, buf, rsa_len); // 释放内存 BN_MONT_CTX_free(mont_ctx); BN_CTX_free(ctx); free(buf); return ret; } // 使用Binary模幂实现RSA加密 int rsa_encrypt_binary(const unsigned char *in, int inlen, unsigned char *out, RSA *rsa, int padding) { int rsa_len = RSA_size(rsa); unsigned char *buf = (unsigned char *)malloc(rsa_len); // 对输入数据进行加密 int ret = RSA_public_encrypt(inlen, in, buf, rsa, padding); // 将结果拷贝到输出缓冲区 memcpy(out, buf, rsa_len); // 释放内存 free(buf); return ret; } // 测试RSA加密 int main() { RSA *rsa = RSA_generate_key(2048, RSA_F4, NULL, NULL); unsigned char in[] = "hello world"; unsigned char out_montgomery[256]; unsigned char out_binary[256]; int outlen_montgomery = rsa_encrypt_montgomery(in, strlen((char *)in), out_montgomery, rsa, RSA_PKCS1_PADDING); int outlen_binary = rsa_encrypt_binary(in, strlen((char *)in), out_binary, rsa, RSA_PKCS1_PADDING); printf("Montgomery: %d\n", outlen_montgomery); printf("Binary: %d\n", outlen_binary); RSA_free(rsa); return 0; } ``` 从上述代码中可以看出，Montgomery模幂算法需要进行模数转换和Montgomery还原，而Binary模幂算法不需要。在实际应用中，Montgomery模幂算法通常比Binary模幂算法更快，因为它避免了一些不必要的模运算。但是，Montgomery模幂算法的预处理过程需要一些额外的时间和空间，因此在处理小数据时，Binary模幂算法可能更适合。
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值