素数求法

  最近在学校ACM题库中碰到求素数的问题。由于那些题目对时间和空间都是有限制的，所以就搜索了下网上求素数的高效方法，并在这总结下！欢迎在评论中补充

一、常规定义法（判断n是否为素数）

只需用2~n-1中的每个数去除，都不能整除则为素数。

http://baike.baidu.com/view/1454797.htm?fr=aladdin

bool isprime(int n){
     int i;
     if(n==1)return false; //1既不是质数也不是合数
     if(n==2)return true;
     if(n%2==0)return false;
     for(i=3;i*i<=n;i+=2)
     if(n%i==0)return false;
     return true;
}

二、厄拉多塞筛选法

具体：在2的上面画一个圆圈，然后划去2的其他倍数；第一个既未画圈又没有被划去的数是3，将它画圈，再划去3的其他倍数；现在既未画圈又没有被划去的第一个数是5，将它画圈，并划去5的其他倍数……依次类推，一直到所有小于或等于N的各数都画了圈或划去为止。这时，表中画了圈的以及未划去的那些数正好就是小于N的素数。

就是一个画出素数，划去素数倍数的过程

在计算机中，筛法可以用给数组单元置零的方法来实现。具体来说就是：首先开一个数组：a[i]，i=1，2，3，…，同时，令所有的数组元素都等于下标值，即a[i]=i，当i不是素数时，令a[i]=0 。当输出结果时，只要判断a[i]是否等于零即可，如果a[i]=0，则令i=i+1，检查下一个a[i]。

#include<iostream>
using namespace std;

int main(){     //求1000内的素数 
    int n[1001];
    for(int i=2;i<=1000;i++)n[i]=i;
    for(int j=2;j<=1000/2;j++){
            if(n[j]!=0)
            for(int k=2;k*j<=1000;k++)
            n[k*j]=0;
            }//for
    for(int i=2;i<=1000;i++){
            if(n[i]!=0)
            cout<<n[i]<<" ";
            }
    system("pause");
    return 0;
}

三、6N±1法

任何一个自然数，总可以表示成为如下的形式之一：
                             6N，6N+1，6N+2，6N+3，6N+4，6N+5 (N=0，1，2，…)

除2和3外，只有6N+1或6N+5(6N-1)可能是素数,所以可以先判断一个数是否是6N±1的形式，如果不是，直接pass。

#include<iostream>
using namespace std;
bool isprime(int n){  //筛选后用第一种方法
     int i;
     if(n==1)return false; 
     if(n==2)return true;
     if(n%2==0)return false;
     for(i=3;i*i<=n;i+=2)
     if(n%i==0)return false;
     return true;
}

int main(){     //判断某个数是否为素数 
    int n;
    cin>>n;
    if(((n-1)%6==0)||((n+1)%6==0))	//筛选
    cout<<isprime(n);
    else if(n==2||n==3)
    cout<<true;
    else cout<<false;
    system("pause");
    return 0;
}