【信息安全数学基础】Fermat素性检测

 1 extern "C"
 2 {
 3 #include<miracl.h>
 4 #include<mirdef.h>
 5 }
 6 #include <stdio.h>
 7 
 8 int main() 
 9 {
10     FILE *fp;
11     int i = 0, k, probability = 1;//k是安全参数，pr是概率
12     printf("请输入安全参数k：");
13     scanf("%d",&k);
14     big m, r, g, a, b, m_4;
15     miracl *mip = mirsys(4000, 10);//初始化MIRACL系统,4000位的10进制数
16     m = mirvar(0);//big类型初始化函数
17     m_4 = mirvar(0);
18     r = mirvar(0);
19     g = mirvar(0);
20     a = mirvar(0);
21     b = mirvar(1);
22     fp = fopen("20.txt", "r+");
23     mip->IOBASE = 10; 
24     cinnum(m, fp);//将大数字符串转换成大数
25     fclose(fp);
26     //奇数m,若取任意整数2<=a<=m-2，gcd（a，m）=1，使得a^(m-1)=1（mod m），则m有1/2的概率为素数
27     for (i = 0;i<k;i++) 
28     {
29         probability = probability * 2;
30         decr(m, 4, m_4);//m_4=m-4
31         bigrand(m_4, a);//产生一个小于m-4的大数随机数a
32         printf("随机数a=");
33         cotnum(a, stdout);
34         incr(a, 2, a);//a=a+2，这样可以保证a>=2且a<=m-2
35         decr(m, 1, g);//g=m-1，作为计算模的指数
36         powmod(a, g, m, r);//r=a^g （mod m），如果r=1，那么m就有二分之一的概率为素数
37         egcd(a, m, g);//g=gcd(a,m)
38         if (mr_compare(g, b) || mr_compare(b, r))//如果gcd(a,m)不等于1或模值r不等于1，不是素数，跳出
39 {
40             printf("gcd(a.m)不等于1或r不等于1\n");
41             break;
42         }
43         else
44             printf("最大公约数为1，模值为1，因此输入有%d/%d可能为素数\n",probability - 1, probability);
45     }
46     printf("由此可得所输入的数字：\n");
47     cotnum(m, stdout);//输出m
48     if (i == k)
49         printf("有%d/%d的可能性输入的是素数\n",probability - 1, probability);
50     else 
51         puts("不是素数\n");
52     mirexit();
53     return 0;
54 }