信息安全数学基础（18）模重复平方计算法

已于 2024-09-30 00:23:39 修改
阅读量710 收藏 8
点赞数 23
分类专栏： 信息安全数学基础 文章标签： 信息安全 网络空间安全 模重复平方计算法
于 2024-09-19 00:30:00 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_73399576/article/details/142336819
版权

前言 

       模重复平方计算法（Modular Exponentiation by Repeated Squaring）是一种用于高效计算形如 abmodn 的数值的算法，其中 a、b、n 均为整数。该算法通过减少乘法次数和利用模运算的性质来显著提高计算效率，特别适用于大数幂运算。

一、基本原理

       模重复平方算法的基本思路是将指数 b 转换为二进制形式，并利用指数的幂次规律与模运算的乘法分配律来逐步计算 abmodn。具体地，它将 b 表示为 b=b0​+2b1​+22b2​+⋯+2kbk​（其中 bi​ 为 0 或 1），然后分别计算 a2imodn（i=0,1,2,…,k），并根据 bi​ 的值决定是否将这些结果相乘并取模 n。

二、算法步骤

  1. 将指数 b 转换为二进制形式：从右至左（即从低位到高位）依次确定 b 的每一位是否为 1，并记录下这些位的位置。

  2. 初始化结果和底数的幂：令结果 res=1，底数的幂 base=a（注意，这里 base 初始化为 a，但实际上在算法执行过程中，我们会不断计算 base 的更高次幂并取模 n）。

  3. 遍历指数的二进制表示：从左至右（即从高位到低位）遍历指数的二进制数，对于每一位 bi​（i 从 0 开始）：

    • 将 base 平方并取模 n：base=base2modn。
    • 如果 bi​=1，则将 res 乘以 base 并取模 n：res=res×basemodn。

  4. 返回最终结果：遍历结束后，res 即为 abmodn 的值。

三、算法复杂度

       模重复平方算法的时间复杂度为 O(logb)，其中 b 是指数的位数。这是因为算法只需要遍历指数的二进制表示中的每一位，而指数的二进制表示的长度与 b 的位数成正比。相比于朴素的幂次计算方法（需要 b−1 次乘法），模重复平方算法大大提高了计算效率。

四、示例

   假设要计算 2^13mod7，首先将指数 13 转换为二进制数 1101：

  • 初始时，res=1，base=2。
  • 遍历到最低位（b0​=1），base=22mod7=4，res=res×basemod7=1×4mod7=4。
  • 遍历到次低位（b1​=0），base=42mod7=2，不更新 res。
  • 遍历到第三位（b2​=1），base=22mod7=4，res=res×basemod7=4×4mod7=2。
  • 遍历到最高位（b3​=1），base=42mod7=2，res=res×basemod7=2×2mod7=4。
  • 最终，res=4，即 213mod7=6（注意，这里的最后结果应为 6，但计算过程中存在笔误，正确过程应为 res 最终保持为 2，然后 2×4mod7=6）。

结论

       模重复平方计算法是一种高效计算模幂的算法，它通过减少乘法次数和利用模运算的性质来优化计算过程。该算法在密码学、大数运算等领域有着广泛的应用。

 结语 

迎难而上

让生命因奋斗而精彩

！！！

IT 青年
关注
专栏目录
模重复平方算法
魏尔肖的博客
04-24 8585
今天在做diffie_hellman加密中需要求一个数很多次方模另一个数，如 44^36mod97,这时用模重复平方算法就可以快速计算出结果。#include <iostream> using namespace std;size_t repeatMod(size_t base, size_t n, size_t mod) { size_t a = 1; while(n)
信息安全数学基础-模重复平方计算法课件
03-20
《信息安全数学基础——模重复平方计算法》 在信息安全领域，数学基础扮演着至关重要的角色，其中模重复平方计算法（Modular Exponentiation）是一种常见的算法，广泛应用于加密算法如RSA、 Diffie-Hellman等。这个...
模重复平方计算法
热门推荐
yxtxiaotian的专栏
09-07 2万+
一、背景： RSA算法里头经常要用到“求x的n次方模m”这样的过程，通常使用O(log(n))的模重复平方算法来实现，提高效率。在《数论》和《信息安全数学基础》中也会介绍。 二、在介绍模重复平方算法之前，我们需要了解模的基本运算，如下（来着百度百科）： 基本概念 给定一个正整数    ，任意一个整数    ，一定存在等式    ； 其中    、
《信息安全数学基础》 张金全 课后习题答案
12-24
### 《信息安全数学基础》张金全 课后习题答案解析 #### 第1章作业参考答案 **题目1:** 设$a$为自己的学号，$b=210$，求整数$s$、$t$，使得$sa+tb=(a,b)$。 **解析:** 这道题目考查了如何利用欧几里得算法（辗转...
信息安全数学基础M2.pptx
03-16
《信息安全数学基础》第二章主要探讨的是同余理论，这是密码学和信息安全中的核心概念。同余是整数算术中的一个重要概念，它涉及到整数的除法和余数。 同余的基本概念可以用数学符号来表示，如果两个整数a和b满足a ...
信息安全中模重复平方算法和欧几里得算法
06-21
在信息安全领域，模重复平方算法（Modular Exponentiation，也称为快速幂运算）和欧几里得算法（Euclidean Algorithm）是两个重要的数学工具，它们在加密算法、数字签名和密钥交换等方面发挥着关键作用。下面将详细...
网络空间安全数学基础笔记1
08-03
模重复平方计算法是安全数学的重要算法，用于解决同余式的解，包括模重复平方计算法的定义、模重复平方计算法的应用等内容。 RSA加密是安全数学的重要应用，用于解决加密问题，包括RSA加密的定义、RSA加密的原理、...
模重复平方算法 C语言
06-05
信息安全求原根时要用到的模重复平方算法，是书上第一种算法的C语言实现
信息安全数学基础——模重复平方计算法（两种方法实现C+JAVA）
m0_46625346的博客
10-30 9485
前言 正在学习信息安全数学基础这门课程，老师让编程实现模重复平方计算法。 先介绍一下模重复平方算法 一、模重复平方算法 图片来自大学MOOC陈恭亮老师的视频中，稍有模糊。。。 下面我们用C语言来实现这个算法。我们可以思考一下，这个算法一共可以分为几个小部分？ 首先，我们需要一个十进制转化为二进制的函数，主要用于将指数转化为二进制形式； 我们还需要知道，我们的循环计算需要执行多少次，所以我们需要一个函数来确定指数的位数； 最后一个函数就是我们的模重复平方计算函数了。当然我们还可以来定义一个打印函
递归实现模重复平方计算法
lkzhang的博客
05-05 1354
求一个大数的较高次幂模运算，可以将其指数分解成二进制的形式进行分解。 举一个简单的例子： 我们要求501^13(mod 667)，这时13可以分解成2进制的形式，13= 2^3 + 2^2 + 2^0。 （此图是陈恭亮《信息安全数学基础的书上图，侵删） 使用递归实现的代码如下： #include <stdio.h> #include <stdlib.h> int mod...
模重复平方计算法的C实现
qq_43100483的博客
10-23 1973
模重复平方计算法（C语言版本） \quadRSARSARSA算法是197819781978年由R.RivestR.RivestR.Rivest、A.ShamirA.ShamirA.Shamir和L.AdlemanL.AdlemanL.Adleman提出的一种用数论构造的、也是迄今为止理论上最为成熟完善的公钥密码密码体制，该体制已得到广泛的应用。 算法描述 1.密钥的产生 （1）选两个保密的大素数ppp和qqq; （2）计算n=p∗qn=p*qn=p∗q，ψ(n)=(p−1)(q−1)\psi(n)=(p-1
RSA | 模重复平方算法 | C语言实现
JacobTsang的博客
06-22 2422
#include <stdio.h> /****************************************************************************** 模重复平方计算法 博文网址 http://blog.csdn.net/yxtxiaotian/article/details/52464496 幂模算法（模重复平方算法），时间复杂度O(...
实验二 模重复平方法
qq_45755193的博客
12-10 2438
模重复平方法
信息安全数学基础（二）：欧拉函数
qq_44109982的博客
12-17 1473
这篇博客将复习欧拉函数的定义及性质，然后给出三个关于欧拉函数的重要定理，最后介绍一种加速平方运算的方法（模重复平方计算法） 一.欧拉函数的性质 定理：设m1,m2是互素的两个正整数，如果k1,k2分别遍历模m1和模m2的简化剩余系，则 m2k1+m1k2遍历模m1,m2的简化剩余系 所以对于欧拉函数： 设m,n是互素的两个正整数，则f(m * n)=f(m)*f(n) 欧拉函数值的计算： f(m)=m * Π (1-1/pk) 对于素数的幂次，其实就是f(pn)=pn-pn-1 二.欧拉定理 设m是大于1的
密码学-模幂运算
m0_51531740的博客
05-11 8572
模幂运算是对给定的整数p、n、a，计算，这个运算在密码学中应用极为普遍，RSA、ElGamal、DH交换等重要密码方案中都涉及模幂运算。 1、算法原理 模重复平方计算法： 算法步骤： ①将n表示为二进制，，用一个数组保存n的各个数位，n共有r+1位。 ②循环计算（用数组b保存所有中间结果）：先循环计算一系列结果：a^2mod p，a^4mod p，a^8mod p，…，再把指数n的二进制表示中取值为1的位对应的a的幂相乘，即可得到最终结果。 2、代码...
用C语言代码实现RSA算法的参数生成，私钥和公钥的计算，以及其中利用模重复平方、蒙哥马利算法、中国剩余定理加速RSA加密和解密，代码需要注释
最新发布
05-30
RSA算法是一种非对称加密算法，其安全性基于大质数的难以分解性。该算法包含参数生成、私钥和公钥计算、加密和解密四个步骤。以下是用C语言实现RSA算法的代码，包括注释。 ```c #include<stdio.h> #include<stdlib.h> #include<math.h> #include<time.h> // 定义大数类型 typedef unsigned long long u64; // 定义结构体存储RSA参数 typedef struct rsa_key_struct { u64 p; // 大素数p u64 q; // 大素数q u64 n; // 模数n=p*q u64 phi; // 欧拉函数值phi=(p-1)*(q-1) u64 e; // 公钥e u64 d; // 私钥d } RSA_KEY; // 判断是否为素数 int is_prime(u64 n) { if (n == 2) return 1; if (n == 1 || n % 2 == 0) return 0; u64 i, m = sqrt(n); for (i = 3; i <= m; i += 2) { if (n % i == 0) return 0; } return 1; } // 生成大素数 u64 generate_prime() { u64 n; do { // 生成随机数 n = (u64)rand() << 32 | rand(); // 取奇数 n |= 1; } while (!is_prime(n)); return n; } // 求最大公约数 u64 gcd(u64 a, u64 b) { u64 r; while (b != 0) { r = a % b; a = b; b = r; } return a; } // 求扩展Euclid算法的逆元，即a*x≡1(mod b)的最小正整数解x u64 extended_euclid(u64 a, u64 b) { u64 x1 = 1, y1 = 0, x2 = 0, y2 = 1; u64 x, y, q, r; while (b != 0) { q = a / b; r = a % b; x = x1 - q * x2; y = y1 - q * y2; a = b; b = r; x1 = x2; y1 = y2; x2 = x; y2 = y; } return x1; } // 快速模幂运算，即a^b(mod n) u64 modexp(u64 a, u64 b, u64 n) { u64 c = 1; while (b != 0) { if (b & 1) c = (c * a) % n; a = (a * a) % n; b >>= 1; } return c; } // 蒙哥马利算法，将加密和解密转换为模重复平方运算 u64 montgomery_transform(u64 a, u64 n, u64 r, u64 n_inv) { u64 t = (a * r) % n; return (t * n_inv) % r; } // 中国剩余定理加速RSA加密和解密 u64 chinese_remainder_theorem(u64 c, RSA_KEY *key) { u64 mp = c % key->p; // c mod p u64 mq = c % key->q; // c mod q u64 dp = key->d % (key->p - 1); // d mod (p-1) u64 dq = key->d % (key->q - 1); // d mod (q-1) u64 q_inv = extended_euclid(key->q, key->p); // q关于p的逆元 u64 x1 = modexp(mp, dp, key->p); // mp^dp(mod p) u64 x2 = modexp(mq, dq, key->q); // mq^dq(mod q) u64 h = (q_inv * (x1 - x2 + key->p)) % key->p; // (x1-x2)q_inv(mod p) return (x2 + h * key->q) % key->n; // x2+hq(mod n) } // 生成RSA参数 void generate_rsa_key(RSA_KEY *key) { u64 p, q, phi, n, e, d; // 初始化随机数生成器 srand((unsigned int)time(NULL)); // 生成两个大素数p和q do { p = generate_prime(); q = generate_prime(); } while (p == q); // 计算模数n和欧拉函数值phi n = p * q; phi = (p - 1) * (q - 1); // 选择公钥e，要求1<e<phi且e与phi互质 do { e = (u64)rand() % (phi - 2) + 2; } while (gcd(e, phi) != 1); // 计算私钥d，要求e*d mod phi=1 d = extended_euclid(e, phi); if (d < 0) d += phi; // 存储RSA参数 key->p = p; key->q = q; key->n = n; key->phi = phi; key->e = e; key->d = d; } int main() { RSA_KEY key; u64 message = 123456789, ciphertext, plaintext; // 生成RSA参数 generate_rsa_key(&key); // RSA加密 ciphertext = montgomery_transform(message, key.n, key.n + 1, -key.n_inv % key.n); ciphertext = modexp(ciphertext, key.e, key.n); // RSA解密 plaintext = chinese_remainder_theorem(ciphertext, &key); // 输出结果 printf("message: %lld\n", message); printf("ciphertext: %lld\n", ciphertext); printf("plaintext: %lld\n", plaintext); return 0; } ``` 以上代码实现了RSA算法的参数生成、私钥和公钥的计算、利用模重复平方、蒙哥马利算法、中国剩余定理加速RSA加密和解密。其中，生成大素数使用了Miller-Rabin素性检验算法，判断是否为素数使用了Trial Division算法，求最大公约数使用了辗转相除法，求扩展Euclid算法的逆元使用了扩展Euclid算法，快速模幂运算使用了蒙哥马利幂算法。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

IT 青年

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值