//素数判断有很多方法,试除法是平时常用的方法,筛选法我只举个埃拉托色尼筛选法和欧拉筛选,当然为了快还有一些其它的筛选方法,根据素数的分布来判断素数也是一种很有效率的方法,Miller-Rabin算法是针对很大的数来讨论的(当然小一点的数也可以判断)
//试除法判断素数
bool JudgePrime(int n)
{
if(n==2)
return true;
if(n%2==0)
return false;
for(int i=3;i*i<=n;i++)
if(n%i==0)
return false;
return true;
}
//埃拉托色尼筛选法--用空间换时间(申请全局变量最好)
bool prime[1000]; //用布尔数组是为了节约空间,一千万个数也就那么10M左右...
void IsPrime()
{
prime[0]=false;
prime[1]=false;
for(int i=2;i<1000;i++)
prime[i]=true;
for(int i=2;i<1000;i++){
if(prime[i]==true){
for(int j=i*2;j<1000;j+=i) //j=i*i数据容易溢出
prime[j]=false;
}
}
}
//欧拉筛选 //网上摘抄的,有空有需要再翻译一下
const int MAXN=3000001;
int prime[MAXN];//保存素数
bool vis[MAXN];//初始化
int Prime(int n)
{
int cnt=0;
memset(vis,0,sizeof(vis));
for(int i=2;i<n;i++)
{
if(!vis[i])
prime[cnt++]=i;
for(int j=0;j<cnt&&i*prime[j]<n;j++)
{
vis[i*prime[j]]=1;
if(i%prime[j]==0)//关键
break;
}
}
return cnt;//返回小于n的素数的个数
}
//根据素数的分布去判断
bool IsPrime(int n)
{
if(n==2 || n==3)
return true;
if(n%6!=1 && n%6!=5)
return false;
for(int i=5;i*i<=n;i+=6)
if(n%i==0 || n%(i+2)==0)//为什么判断i,i+2;应该是因为每一个合数它拆分最后都能拆成素数的乘积
return false;//我猜是这个原因,然而因为素数分布都在6倍数的两旁,所以判断i,i+2也就相当于判断素数
return true;
}
//Miller-Rabin素性测试
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <map>
using namespace std;
int a[5]={2,3,7,61,24251};
//如果选用2, 3, 7, 61和24251作为底数,那么10^16内唯一的强伪素数为46856248255981
long long Random( long long n ) //生成[ 0 , n ]的随机数,用随机数也可以
{
return ((double)rand( ) / RAND_MAX*n + 0.5);
}
long long FastExpMod(long long a,long long b,long long c) //快速幂取模 a^b%c
{
long long ans=1;
while(b){
if(b&1)
ans=(ans*a)%c;
a=(a*a)%c;
b>>=1;
}
return ans;
}
bool witness(int a,long long n) //正向迭代
{
long long p=n-1;
long long x=FastExpMod(a,p,n);
if(x==n-1) return true;
else
{
long long w;
do
{
p>>=1;
w=FastExpMod(a,p,n);
if(w==n-1) return true;
else if(w!=1) return false;
}
while((p&1)!=1);
if(w==1||w==n-1) return true;
else return false;
}
}
bool miller_rabin( long long n ) //检验n是否是素数
{
if(n == 2)
return true;
if(n < 2 || n % 2 == 0)
return false; //如果是2则是素数,如果<2或者是>2的偶数则不是素数
// for(int i = 1; i <=10; i++) //做times次随机检验
// {
// long long a = Random( n - 2 ) + 1; //得到随机检验算子 a
// printf("%d\t",a);
// if(!witness( a, n )) //用a检验n是否是素数
// return false;
// }
// return true;
for(int i=0;i<5;i++)
if(!witness(a[i],n))
return false;
return true;
}
int main()
{
int n;
while(scanf("%d",&n)!=EOF){
if(n==2 || n==3 || n==7 || n==61 || n==24251)
printf("yes\n");
else if(miller_rabin(n))
printf("yes\n");
else
printf("no\n");
}
return 0;
}
额,就这样.