某电的信息安全数学基础综合实验
在vs2017下使用miracl库做大素数的素性检测算法
环境搭建请参考https://blog.csdn.net/qq_42450533/article/details/102493504
直接贴代码了
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include "miracl.h"
int main()
{
FILE *fp;//大数文件
char fpname[100];//fpname为大数文件的文件名
int k,num;//k为安全参数,num为循环的次数
double possible;//possible为待检测数为素数的可能性
miracl *mip = mirsys(500, 10);//初始化miracl内存空间,为500位的10进制数
big m = mirvar(0);//m为从大数文件中读入的大数
big b = mirvar(0);//b为随机生成的大数
big g = mirvar(0);//g为b和m的最大公约数
big r = mirvar(0);//r为b的m-1次方 (mod m)
big hhh = mirvar(0);//hhh为一个中间数
big one = mirvar(1);//one为1
big two = mirvar(2);//two为2
printf("输入大数文件完整路径\n");
scanf("%s",fpname);
if ((fp = fopen(fpname, "r")) == NULL)
{
printf("fail to open the file\n");
system("pause");
exit(0);
}
else {
cinnum(m, fp);//从fp文件中把大数赋值给big型变量m
//system("pause");
cotnum(m, stdout);//把 m 打印出来
printf("\n\n");
printf("reading success\n");//成功读取
printf("\n");
printf("input k:");//输入安全传输k
scanf("%d", &k);
printf("\n");
//开始进行素数判断
for (num = 1; num <= k; num++)
{
//开始生成一个随机数b,要求2<=b<=m-2
decr(m, 1, hhh);//hhh=m-1
while (1)
{
bigrand(hhh, b);//生成随机数b,使b<hhh,即b<m-1
if (compare(b,two)>=0)break;//随机数b>=2时跳出循环加入下一部分
}
//开始计算g 并判断
egcd(m, b, g);//g=(m,b)
if (compare(g, one))//若g=a,返回0,若g>a,返回1
{
printf("m与随机数的最大公约数不为1,m为合数\n");
system("pause");
exit(0);
}
//开始计算r, 并判断
powmod(b, hhh, m, r);//r=b的hhh(m-1)次方 (mod m)
if (compare(r, one))//若r=a,返回0,若r>a,返回1
{
printf("m与随机数生成的 g 不为1,m为合数\n");
system("pause");
exit(0);
}
}
//经过k轮检验,m不是合数
possible = 1 - 1 / (pow(2, k));//计算概率
printf("m为素数的概率是 %lf ,安全参数为 %d \n\n", possible,k);
system("pause");
fclose(fp);
mirexit();
}
return 0;
}
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。