Miller_Rabbin 素数测试

最新推荐文章于 2024-11-01 17:22:42 发布

原创最新推荐文章于 2024-11-01 17:22:42 发布 · 721 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#Prime MillerRabbin

ACM 专栏收录该内容

7 篇文章

订阅专栏

本文介绍了一个基于Miller_Rabbin算法实现的素数测试程序。该程序使用C语言编写，通过快速幂模运算来判断一个数是否为素数，并采用随机数生成提高测试准确性。时间复杂度为O((logn)^3)。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

/***********************************************
名称：Miller_Rabbin素数测试程序
作者：ForrestGump
时间：2013-05-05
时间复杂度：O((logn)^3)
版本：v1.0
***********************************************/
#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#define s 100
double random()//生成[0~1]之间的均匀随机数
{
	return (double)rand()/RAND_MAX;
}
int random(int m)//生成[1~m-1]之间的均匀随机数
{
	return (int)(random()*(m-2)+1.5);
}
__int64 quick_mod(__int64 a,__int64 n,__int64 m)//快速求a的b次方mod c
{
   __int64 d,t;
    d=1;
    t=a;
    while (n>0)
    {
        if (n%2==1)//模拟二进制
            d=(d*t)%m;
        n/=2;
        t=(t*t)%m;
    }
	return d;
}

bool Miller_Rabbin(__int64 n)//素数测试
{
	for(int i=0;i<s;i++)
	{
		__int64 a=(__int64)random(n);
		if(quick_mod(a,n-1,n)!=1)
		{
			return false;
		}
	}
	return true;
}
int main()
{
	__int64 n;
	while(scanf("%I64d",&n)!=EOF)
	{
		srand(time(NULL));//初始化随机数种子
	    if(Miller_Rabbin(n)&&n!=1)
		    printf("true\n");
	    else printf("false\n");
	}
}