素数测试方法

方法一

埃拉托斯尼斯筛选法：正整数n是素数，当且仅当她不能被任何一个小于根号n的素数整除

解释：一个数可以分解为一个素数乘以一个整数相乘的形式，素数介于[1,sqrt(n)], 另一个整数介于[sqrt(n),n].

代码实现：筛选小于N的所有素数并存储

bool isprime[N];///判断是否为素数
int prime[N];///顺序存储得到的素数
int nprime;///记录素数的个数
void doprime()
{
    int i,j;
    nprime=0;
    memset(isprime,1,sizeof(isprime));
    isprime[1]=0;
    for(i=2;i<N;i++)
    {
        if(isprime[i])
        {
            nprime++;
            prime[nprime]=i;
            for(j=i*i;j<N;j+=i)///素数的倍数不是素数
                isprime=0;
        }
    }
}

方法二

任何一个自然数总可以表示成如下形式之一

6N,6N+1,6N+2,6N+3,6N+4,6N+5;(N=0，1，2，3.。。。)，显然，只有形如6N+1,6N+5的自然数有可能是素数，除了2和3之外，所有的素数都可以表示成6N+1或6N-1的形式

代码实现如下：筛选小于MAXN的所有素数并存储

 int prime[MAXN];///顺序存放筛选得到的素数
 int nprime;///素数的个数
 bool isprime(int k)
 {
      if(k%2==0)
        return false;
     for(int i=3;i*i<k;i+=2)
        if(k%i==0)
        return false;
     return true;
 }
 void doprime()
 {
     nprime=0;
   prime[++nprime]=2;
     prime[++nprime]=3;
     for(int i=1;i<MAXN;i+=6)
        for(int j=-1;j<=1;j+=2)
        if(isprime(i+j))
        prime[++nprime]=i+j;
 }