【网络安全】【密码学】【北京航空航天大学】实验三、数论基础(下)【C语言实现】

实验三、数论基础(下)

一、实验内容

1、中国剩余定理(Chinese Remainder Theorem)

(1)、算法原理

m1, m2, … mk 是一组两两互素的正整数,且 M = m1 · m2 · … · mk 为它们的乘积, 则如下的同余方程组:
x == a1 (mod m1)
x == a2 (mod m2)

x == ak (mod mk)

对于模M有唯一的解 x = (M · e1 · a1 / m1 + M · e2 · a2 / m2 + … + M · ek · ak / mk) (mod M)
其中 ei 满足 M / mi · ei == 1(mod mi)

(2)、算法流程

本算法的大致流程如下图所示:

在这里插入图片描述

(3)、算法的代码实现(C语言)

#include <stdio.h>int reverse(int k, int m);  // 函数,返回k模m的逆元int main()
{int i;int r;       // 方程组中的方程个数 (不能超过100)int b[100];  // 余数数组int m[100];  // 模数数组int mul = 1;int M[100];  // M数组int M1[100];  // M'数组int x = 0;  // 方程组的根//	printf("%d", reverse(3, 7));  // 一行测试代码printf("请输入方程的个数:");scanf_s("%d", &r);  // 选用安全的输入函数,避免可能的栈溢出(攻击)printf("请输入 %d 个余数,之间以空格分隔:", r);for(i = 0;i < r;i ++){scanf("%d", &b[i]);}printf("请输入 %d 个模数,之间以空格分隔:", r);for(i = 0;i < r;i ++){scanf("%d", &m[i]);mul *= m[i];}for(i = 0;i < r;i ++){M[i] = mul / m[i];}for(i = 0;i < r;i ++){M1[i] = reverse(M[i], m[i]);}for(i = 0;i < r;i ++){x += M1[i] * M[i] * b[i];}x %= mul;printf("此同余方程组的解(模%d)是:", mul);printf("%d", x);return 0;
}int reverse(int k, int m)
{int i;for(int i = 1;i < m;i ++){if(k * i % m == 1){return i;}}return -1;
}

(4)、算法测试

测试点1:

x == 1 (mod 4)
x == 2 (mod 5)
x == 3 (mod 7)

运行时截图:

在这里插入图片描述
解为 x == 17 (mod 140)

测试点2:

x == 7 (mod 23)
x == 9 (mod 28)
x == 16 (mod 33)

运行时截图:

在这里插入图片描述

解为 x == 19189 (mod 21252)

测试点3:

x == 23 (mod 283)
x == 28 (mod 102)
x == 33 (mod 35)

运行时截图:
在这里插入图片描述

解为 x == 43888 (mod 1010310)

2、素性检测算法(Miller-Rabin’s Test for Primality)

(1)、算法原理
根据费马小定理,设 p素数a整数,且满足 (a, p) = 1, 则满足 a ^ (p - 1) = 1 (mod p), 以及二次探测定理:如果 p 是一个素数,且 0 < x < p, 且同余方程 x ^ 2 = 1 (mod p) 成立,那么 x = 1x = p - 1米勒·拉宾 Miller-Rabin 素性检测算法是基于以上两个定理的随机化算法,用于判断一个整数是合数还是素数。

(2)、算法流程

本算法的大致流程如下图所示:

在这里插入图片描述

(3)、算法的代码实现(C语言)

#include <stdio.h>
#include <stdlib.h>typedef long long unsigned LLU;
typedef int BOOL;#define TRUE 1
#define FALSE 0// 长整数快速模乘算法
LLU quickMult(LLU a, LLU b, LLU c)
{LLU result = 0;while(b > 0) {if(b & 1)result = (result + a) % c;a = (a + a) % c;b >>= 1;}return result;
}// 长整数快速幂取模算法
LLU quickPower(LLU a, LLU b, LLU c) 
{LLU result = 1;while(b > 0) {if(b & 1)result = quickMult(result, a, c);a = quickMult(a, a, c);b >>= 1;}return result;
}// 米勒·拉宾素性检验算法(单次测试)
BOOL MillerRabinPrimeTest(LLU n) 
{LLU d, x, newX, a = 1;int i;for (i = 0; i < 4; i ++)a *= rand();a = a % (n - 3) + 2;  // 随机地选取一个a∈[2,n-2]int s = 0;  // s为d中的因子2的幂次数。d = n - 1;while ((d & 1) == 0) {   // 将d中的因子2全部提取出来。s ++;d >>= 1;}x = quickPower(a, d, n);for (i = 0; i < s; i ++) { // 进行s次二次探测newX = quickPower(x, 2, n);if (newX == 1 && x != 1 && x != n - 1)return FALSE;  // 用二次定理的逆否命题,此时n被确定为合数。x = newX; }if (x != 1)return FALSE;  // 用费马小定理的逆否命题判断,此时x=a^(n-1) (mod n),那么n确定为合数。return TRUE; //用费马小定理的逆命题判断。能经受住考验至此的数,大概率为素数。
}//经过连续特定次数的Miller-Rabin测试后,
//如果返回值为TRUE表示n为素数,返回值为FALSE表示n为合数。
BOOL isPrimeByMR(LLU n) 
{if((n & 1) == 0)return FALSE;int i;for (i = 0; i < 100; i ++)if(MillerRabinPrimeTest(n) == FALSE)return FALSE;return TRUE;
}// 主函数
int main()
{LLU n;printf("请输入待判断素性的整数:");scanf("%lld", &n);BOOL result;result = isPrimeByMR(n);printf("\n------判断中......------\n\n");if(result == TRUE)printf("%llu 是素数", n);elseprintf("%llu 是合数", n);return 0;
}

(4)、算法测试

测试点1:
判断1000023是素数还是合数。(答:合数

运行时截图:

在这里插入图片描述

测试点2:
判断1000033是素数还是合数。(答:素数

运行时截图:

在这里插入图片描述
测试点3:
判断100160063是素数还是合数。(答:合数

运行时截图:

在这里插入图片描述测试点4:
判断1500450271是素数还是合数。(答:素数

运行时截图:

在这里插入图片描述

说明:算法为概率性判断,即可能将合数错判为素数(对计算机来说,已在极短的时间内完成了100次重复的MR测试,故该错判的概率极低),但绝无可能将素数错判为合数。

二、参考文献

1、《密码编码学与网络安全——原理与实践(第七版)》(Cryptography and Network Security, Principles and Practice, Seventh Edition),【美】威廉 斯托林斯 William Stallings 著,王后珍等 译,北京,电子工业出版社,2017年12月。

2、《密码学实验教程》,郭华 刘建伟等 主编,北京,电子工业出版社,2021年1月。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/344239.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

使用new pm写一个pass

范例来自LLVM Techniques, Tips, and Best Practices Clang and Middle-End Libraries llvm ir到ir是由一个个pass处理的&#xff0c;从一个ir到另一个ir会改变一些东西 书里面就是说想要给指针变量添加一个noalias属性 书里面使用插件的形式&#xff0c;但是不知道怎么我搞不…

MetaGPT前期准备与快速上手

大家好&#xff0c;MetaGPT 是基于大型语言模型&#xff08;LLMs&#xff09;的多智能体协作框架&#xff0c;GitHub star数量已经达到31.3k。 接下来我们聊一下快速上手 这里写目录标题 一、环境搭建1.python 环境2. MetaGpt 下载 二、MetaGPT配置1.调用 ChatGPT API 服务2.简…

手机扫一扫文档、证件表格就能打印,不用专门找扫描软件啦

现在各大品牌的手机扫一扫功能都可以满足扫描文档、扫描表格的功能&#xff0c;平常如果需要将纸质文件扫描成电子档&#xff0c;拿起手机就可以实现。文件签名后再扫描成PDF版本&#xff0c;相信大家都用得不少了。 今天就用华为手机的扫一扫给大家示范一下如何将文档、证件、…

设计模式—行为型模式之状态模式

设计模式—行为型模式之状态模式 状态&#xff08;State&#xff09;模式&#xff1a;对有状态的对象&#xff0c;把复杂的“判断逻辑”提取到不同的状态对象中&#xff0c;允许状态对象在其内部状态发生改变时改变其行为。 状态模式包含以下主要角色&#xff1a; 环境类&am…

Kafka配置Kerberos安全认证及与Java程序集成

Background 本文主要介绍在 Kafka 中如何配置 Kerberos 认证&#xff0c;以及 java 使用 JAAS 来进行 Kerberos 认证连接。本文演示为单机版。 所用软件版本 查看 Kerberos 版本命令&#xff1a;klist -V 软件名称版本jdk1.8.0_202kafka2.12-2.2.1kerberos1.15.1 1、Kerberos …

【STM32】STM32学习笔记-串口发送和接收(27)

00. 目录 文章目录 00. 目录01. 串口简介02. 串口相关API2.1 USART_Init2.2 USART_InitTypeDef2.3 USART_Cmd2.4 USART_SendData2.5 USART_ReceiveData 03. 串口发送接线图04. USB转串口模块05. 串口发送程序示例06. 串口发送支持printf07. 串口发送支持printf_v208.09.10. 01.…

Linux网络基础及bonding实际操作

1.查看linux基础的网络配置 网关&#xff1a;route -n ip地址&#xff1a;ifconfig 或 ip a DNS服务器&#xff1a;cat /etc/resolv.conf 主机名&#xff1a;hostname 路由&#xff1a;route -n 网络连接状态&#xff1a;ss 或 netstat 2.临时修改网卡名称 3.永久修…

在全志T113-i平台上实现H.265视频解码步骤详解

H.265&#xff0c;也被称为HEVC(HighEfficiency Video Coding)&#xff0c;作为H.264的继任者&#xff0c;提供了更好的视频压缩和更高的视频质。H.265通过引入更多先进的编码技术&#xff0c;如更强大的运动估计和更高效的变换编码&#xff0c;对比H.264进行了改进。这些改进使…

宏景eHR 多处 SQL注入漏洞复现

0x01 产品简介 宏景eHR人力资源管理软件是一款人力资源管理与数字化应用相融合,满足动态化、协同化、流程化、战略化需求的软件。 0x02 漏洞概述 宏景eHR view、trainplan_tree.jsp等接口处存在SQL注入漏洞,未经过身份认证的远程攻击者可利用此漏洞执行任意SQL指令,从而窃…

Python web自动化测试框架搭建(功能接口)——通用模块

1、通用模块&#xff1a; config.conf: 公共配置文件&#xff0c;配置报告、日志、截图路径&#xff0c;以及邮件相关配置 [report] reportpath E:\workspace\WebAutomation\src\functiontest\Report\2017-07-18 screen_path E:\workspace\WebAutomation\src\functiontest\R…

vue上传文件加进度条,fake-progress一起使用

el-upload上传过程中加进度条&#xff0c;进度条el-progress配合fake-progress一起使用&#xff0c;效果如下&#xff1a; 安装 npm install fake-progress 在用到的文件里面引用 import Fakeprogress from "fake-progress"; 这个进度条主要是假的进度条&#xff…

数字信号处理 唐向宏著 pdf +课后答案 免费下载

数字信号处理——原理、实现与仿真 pdf 唐向宏著 &#xff0b;课后答案 杭州电子科技大学 费劲心思在网上花钱买的&#xff0c;共享给大家 永久链接&#xff1a;https://wwi.lanzoup.com/b0140pf4f 密码&#xff1a;aflj 里面除了有原书PDF&#xff0c;还有课后题答案