hdu2138 How many prime numbers (素数测试)

最新推荐文章于 2020-07-23 19:25:45 发布

原创最新推荐文章于 2020-07-23 19:25:45 发布 · 687 阅读

1 ·

CC 4.0 BY-SA版权

数论专栏收录该内容

21 篇文章

订阅专栏

本文介绍了一种高效的素数检测方法——米勒-拉宾算法，并提供了详细的C语言实现代码。该算法通过选取多个测试基数来判断一个数是否为素数，大大提高了检测准确率。

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

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2138

米勒拉宾算法

#include <stdio.h>

_int64 pow_mod(_int64 a,_int64 b,_int64 mod)
{
	_int64 ans=1;
	while(b)
	{
		if(b&1)
			ans=(ans*a)%mod;
		a=(a*a)%mod;
		b>>=1;
	}
	return ans;
}

int MillerBabbin(_int64 n,_int64 a)
{
	_int64 t;
	_int64 d=n-1;
	if(n==2||n==a)
		return 1;
	if((n&1)==0)//n为偶数
		return 0;
	while(!(d&1))
		d>>=1;
	t=pow_mod(a,d,n);
	while((d!=n-1)&&(t!=1)&&(t!=n-1))
	{
		t=t*t%n;
		d<<=1;
	}
	return ((t==n-1)||(d&1)==1);
}

int isPrime(_int64 n)
{
	int i;
	_int64 test[5]={2,3,5,7};
	for(i=0;i<4;++i)
	{
		if(!MillerBabbin(n,test[i]))
			return 0;
	}
	return 1;
}

int main()
{
	int n,ans;
	_int64 x;
	while(scanf("%d",&n)!=EOF)
	{
		ans=0;
		while(n--)
		{
			scanf("%I64d",&x);
			if(isPrime(x))
				ans++;
		}
		printf("%d\n",ans);
	}
	return 0;
}