- 1、本文档共30页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学基础与密码学部分实验重要内容
数学基础与密码学部分实验内容
信息安全数学基础实验部分
题目:定理:设n是一个正整数,如果对所有的素数p≤ QUOTE ,都有p?n,则n一定是素数。
注:古希腊数学家埃拉托斯散(Eratosthenes,公元前275—公元前194)发明了求
比某给定数小的素数的筛法技巧。
方法如下:
对于任意给定的正整数N,要求出所有不超过N的素数。我们列出N个整数,从中删除小于等于 QUOTE 的所有素数p1,…,pk的倍数。然后依次删除,
p1的倍数:2p1,…, QUOTE p1。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。pk的倍数:2pk,…, QUOTE pk
余下的整数(不包括1)就是所要求的不超过N的素数。
使用VC++编程语言编写一个可测定不超过1,000,000的素数判定程序。
程序设计思想:
有题目所给定的定理我们可以采用:首先将整数N下对应的小于 QUOTE 所有的素数尽,则通过 QUOTE 的倍数比不为素数,并将此种类型的整数提出后,最终只剩下素数输出到屏幕上。
程序流程图:
程序源代码:
#includestdio.h
int prime(int m)
{
int n,i;
n=m/2;
for(i=2;i=n;i++)
if(m%i==0)
return 0;
return 1;
}
int main(void)
{
int k;
for(k=2;k=1000000;k++)
if(prime(k)) printf(%d ,k);
return 0;
}
程序运行截图:
题目二 :
使用VC++编程语言设计实现一个算法程序库,要求包括以下部分:
1)欧几里德算法求a,b的最大公倍数;
2)扩展的欧几里德算法,求出gcd(a,b)和满足gcd(a,b)=ax+by的整数x和y;
3)求解模线性方程 ax ≡ b (mod n) 其中n0;
4)求解模线性方程组(中国余数定理);
5)模取幂运算,计算ab mod n (a,b1032);
6)Miller-Rabin随机性素数测试算法(要求判定n1016);
1)欧几里德算法求a,b的最小公倍数;
程序设计思想:
采用欧几里得思想的推论,两个整数的积相乘并除以他们的最大公约数可以求出这两个整数的最小公倍数。
程序流程图:
程序源代码:
#includestdio.h
int main()
{
int a,b,c;
printf(请依次输入两个数:\n);
scanf(%d %d,a,b);
int a_cup, b_cup, res; /*被除数, 除数, 余数*/
if (a 0 b 0)
{
a_cup = a;
b_cup = b;
res = a_cup % b_cup;
while (res != 0)
{
a_cup = b_cup;
b_cup = res;
res = a_cup % b_cup;
}
c=a*b/b_cup;
}
printf(%d\n,c);
return 0;
}
运行结果:
2)扩展的欧几里德算法,求出gcd(a,b)和满足gcd(a,b)=ax+by的整数x和y;
#includestdio.h
int x,y,q;
void extend_Eulid(int a,int b)
{
if(b == 0)
{
x = 1;y = 0;q = a;
}
else
{
extend_Eulid(b,a%b);
int temp = x;
x = y;
y = temp - a/b*y;
}
}
int main()
{
int a,b;
printf(Please input two numbers:);
scanf(%d %d,a,b);
if(ab)
{
int m=a;
a=b;
b=m;
}
else extend_Eulid(a,b);
printf(%d=(%d)*%d+(%d)*%d\n,q,x,a,y,b);
return 0;
}
3)求解模线性方程 ax ≡ b (mod n) 其中n0;
源代码:
#includeiostream
#includestring
#includetime.h
using namespace std;
#includeBig
您可能关注的文档
最近下载
- 高清晰全欧洲铁路网地图.pdf
- 组建创业团队.ppt VIP
- 常见的新生儿高频振荡通气(周伟).ppt
- 第九篇:同红军在一起(续)-初中语文八年级上册名著《红星照耀中国》导读系列课件.pptx VIP
- 融创首创武汉经开国际智慧生态城市149R2地块项目超高层避难层悬挑支模架专项施工方案.pdf
- 部编版道德与法治六年级上册8《我们受特殊保护》(教案).docx
- 统编版选择性必修1与岳麓版必修(Ⅰ)高中历史教科书比较研究——以“政治制度”单元为例.pdf
- 小学五年级上册道德与法治 独领风骚的古代技术创造 教学设计 .pdf
- 浙江省义乌市小商品出口贸易结构现状存在的问题及对策.docx
- 《专业认知实习(2)》实习教学大纲.docx VIP
文档评论(0)