素数判定总结

1.对于百万级别，判断单个数是否为素数，用埃拉托斯尼斯筛法打一个判断是否为素数的表预处理一下。

代码：

const int N=2000000;

bool isprime[N];

void doprime()
{
	long long i,j;
	for(i=1;i<N;i+=2)
		isprime[i]=true;
	isprime[1]=false;
	isprime[2]=true;
	for(i=2;i<N;i++)
	{
		if(isprime[i])
		{
			for(j=i*i;j<N;j+=i)
				isprime[j]=false;
		}
	}
}

2.在区间【L,R】内查询素数，其中L<=R,R在百万级别，用埃拉托斯尼斯筛法打一个素数表预处理一下。

代码：

const int N=2000000;

int nprime,prime[N]; 
bool isprime[N];

void doprime()
{
	nprime=0;
	long long i,j;
	for(i=1;i<N;i+=2)
		isprime[i]=true;
	isprime[1]=false;
	isprime[2]=true;
	for(i=2;i<N;i++)
	{
		if(isprime[i])
		{
			prime[nprime++]=i;
			for(j=i*i;j<N;j+=i)
				isprime[j]=false;
		}
	}
} 

3.判断一个【1,2^31】内的数是否为素数，根据基本素数判别法。打一个1~50 000的素数表，看该数字是否存在一个1~50 000的素因子。http://blog.csdn.net/w20810/article/details/43196133

代码：

#include <cstdio>  
#include <iostream>  
#include <cmath>  
using namespace std;  
  
#define N 50001  
  
bool isprime[N];  
int prime[N],nprime;  
  
void doprime()  
{  
    long long i,j;  
    nprime=0;  
    for(i=0;i<N;i++)  
        isprime[i]=true;  
    isprime[1]=0;  
    for(i=2;i<N;i++)  
    {  
        if(isprime[i])  
        {  
            prime[++nprime]=i;  
            for(j=i*i;j<N;j+=i)  
                isprime[j]=false;  
        }  
    }  
}  
  
bool isp(int n)  
{  
    int i,k=(int)sqrt(double(n));  
    for(i=1;prime[i]<=k;i++)  
        if(n%prime[i]==0)  
            return false;  
    return true;  
}  

4.对给定区间【L，R】查询素数，其中L<=R，R<2^31。 先打一个1~50 000的素数表，用这个素数表筛选区间【L，R】的素数。http://blog.csdn.net/w20810/article/details/43313261

代码：

    #include <cstdio>  
    #include <string.h>  
    #include <iostream>  
    #include <cmath>  
    using namespace std;  
      
    typedef long long int64;  
      
    bool isprime[1000005];  
    int prime[50001],nprime,use[1000005];  
      
    void doprime()  
    {  
        long long i,j;  
        prime[1]=2;  
        nprime=1;  
        for(i=1;i<50000;i+=2)  
            isprime[i]=true;  
        for(i=2;i<50000;i++)  
        {  
            if(isprime[i])  
            {  
                prime[++nprime]=i;  
                for(j=i*i;j<=50000;j+=i)  
                    isprime[j]=false;  
            }  
        }  
    }  
      
    bool solve(int64 L,int64 R,int64 &a,int64 &b,int64 &c,int64 &d)  
    {  
        int64 i,j,s;  
        for(i=1;i<nprime;i++)  
        {  
            s=L/prime[i];  
            while(s<=1 || s*prime[i]<L) s++;  
            for(j=s*prime[i];j<=R;j+=prime[i])  
                if(j-L>=0)  
                    isprime[j-L]=false;  
        }  
        int np=0;  
        for(i=0;i<=R-L;i++)  
        {  
            if(isprime[i])  
            {  
                use[np++]=i+L;  
            }  
        }  
        if(np<=1)  
            return false;  
              
        int dmax=-1,dmin=2147483647;  
        for(i=0;i<np-1;i++)  
        {  
            if(use[i+1]-use[i]<dmin)  
            {  
                dmin=use[i+1]-use[i];  
                a=use[i];  
                b=use[i+1];  
            }  
            if(use[i+1]-use[i]>dmax)  
            {  
                dmax=use[i+1]-use[i];  
                c=use[i];  
                d=use[i+1];  
            }  
        }  
        return true;  
    }  
      
    int main()  
    {  
        doprime();  
        int64 a,b,c,d;  
        int64 L,R;  
        while(scanf("%lld%lld",&L,&R)!=EOF)  
        {  
            if(L<=1) L++;  
            memset(isprime,true,sizeof(isprime));  
            if(solve(L,R,a,b,c,d))  
                printf("%lld,%lld are closest, %lld,%lld are most distant.\n",a,b,c,d);  
            else  
                printf("There are no adjacent primes.\n");  
        }  
        return 0;  
    }  

5.判断一个【1，2^63】内的数是否为素数，用Miller_rabin算法判断是否为素数。

    #include <iostream>  
    #include <cstdio>  
    #include <cmath>  
      
    using namespace std;  
    #define TIMES 5  
    typedef long long LL;  
      
    LL random(LL n)  
    {  
        return (LL)((double) rand()/RAND_MAX*n+0.5);  
    }  
      
    LL multi(LL a,LL b,LL mod) //a*b%mod  
    {  
        LL ans=0;  
        while(b)  
        {  
            if(b)  
                ans=(ans+a)%mod;  
            b>>=1;  
            a=(a<<1)%mod;  
        }  
        return ans;  
    }  
      
    LL kpow(LL a,LL n,LL mod) //a^n%mod  
    {  
        LL ans=1;  
        while(n)  
        {  
            if(n&1)  
                ans=(ans*a)%mod;  
            n>>=1;  
            a=(a*a)%mod;   
        }  
        return ans;  
    }  
      
    bool witness(LL a,LL n)  
    {  
        LL m=n-1;  
        LL j=0;  
        while(! (m&1))  
        {  
            j++;  
            m>>=1;  
        }  
        LL x=kpow(a,m,n);  
        if(x==1 || x==n-1)  
            return false;  
        while(j--)  
        {  
            x=x*x%n;  
            if(x==n-1)  
                return false;  
        }  
        return true;  
    }  
      
    bool miller_rabin(LL n)  
    {  
        if(n<2) return false;  
        if(n==2) return true;  
        if(!(n&1)) return false;  
        for(LL i=1;i<=TIMES;i++)  
        {  
            LL a=random(n-2)+1;  
            if(witness(a,n))  
                return false;  
        }  
        return true;  
    }  
      
    int main()  
    {  
        LL x;  
        while(cin>>x)  
        {  
            if(miller_rabin(x))  
                cout<<"yes\n";  
            else  
                cout<<"no\n";  
        }  
        return 0;  
    }