筛选素数的几种方法

判断一个数是否为素数：

 

int prime(int n)
{
	if(n==1)
		return 0;
	if(n==2||n==3)
		return 1;
	for(int i=2;i*i<=n;i++)
		if(n%i)
			return 0;
	return 1;
}

这个相信大家都写过，下面介绍一种更快的判断方法：

bool Is_prime(int n)
{
    if(n==1) return false;
    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) return false;
    return true;
}

 证明：令x≥1，将大于等于5的自然数表示如下： ......6x−1,6x,6x+1,6x+2,6x+3,6x+4,6x+5,6(x+1),6(x+1)+1

观察一波我们可以发现，6x、6x+2、6x+3、6x+4、6(x+1)都不是质数，也就是说质数必定等于6x-1或6x+1。

(很粗略orz 可能不是很严谨)

下面介绍三种筛选素数的方法：

(1)Eratosthenes筛法

int vis[MAX];
memset(vis,0,sizeof(vis));
for(int i=2;i<=n;i++)
	for(int j=i*2;j<=n;j+=i)
		vis[j]=1;
vis[1]=1;

思想：一个数的倍数一定不是素数。时间复杂度在O(nlgn)。

下面有一张动图可以帮助你理解。

(2)对Eratosthenes筛法的改进：

int MAX=sqrt(1000000+0.5);
	for(int i=2;i<=MAX;i++)
		if(!vis[i])
			for(int j=i*i;j<=1000000;j+=i)
				vis[j]=1;
vis[1]=1;

那个1000000大家可以任意设置。

(3)欧拉筛素数：

int prime[MAX];
int vis[MAX];
void Prime()
{
    int k=-1;
    for(int i=2;i<=MAX;i++)
    {
        if(!vis[i])//是素数
            prime[++k]=i;
        for(int j=0;j<=k&&i*prime[j]<=MAX;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)
                break;
        }
    }
}

 首先，任何合数都能表示成多个素数的积。所以，任何的合数肯定有一个最小质因子。我们通过这个最小质因子就可以判断什么时候不用继续筛下去了。这也是i%prime[j]==0；break 这句话的作用。举个例子，当i=4时，prime[0]=2，prime[1]=3，内层循环首先筛掉了4*prime[0]=8，这时候4%2==0，应该跳出循环，如果不跳出循环的话，下一次会筛掉4*prime[1]=12，然而这样是没有必要的，因为12会在i=6时被筛掉。也正因为这个语句，欧拉筛素数不会有重复的情况。在线性复杂度内可以完成。

测试以后发发现方法(2)和方法(3)的时间效率非常接近，方法(3)仅仅快一点点。