筛法求素数

更好的阅读体验

素数是常用的一种神奇的数字，记得我读初中的时候就手动计算过一百到两百的素数，计算量还是不小的，现在想想还是蛮有趣的。今天就来讲一下怎么快速的计算素数。

以一百万为例。

最原始的方法就是逐个判断。

#define Max 1000000
bool prime[Max]={false};
int main ()
{
	
	for(int i=1;i<Max;i++)
	{
		bool flag=true;
		for(int j=1;j<i;j++)
			if(i%j==0)
				flag=false;
		if(flag)
			prime[i]=true;
	}
}

以上是最好写但是效率最低的一种写法。对时间稍微要求高一点就出问题了。下面来介绍一种比较快的计算方法，筛法。

我们都知道任何一个正整数都可以唯一质因数分解，假设一个合数分解成了很多个因子，那么显然他是最小因子的整数倍，根据这个原理我们可以将所有的合数被最小的因子分解掉，也就是这里所说的“筛掉”，这就是筛法的主要思想。

#define Max 1000000
bool prime[Max];//用来筛掉合数
int num[Max];//用来存储素数
int main()
{
	memset(prime,false,sizeof(prime));//false意味着是素数
	int count=0;//统计素数个数，并存储素数
	for(int i=2;i<Max;i++)
		if(!prime[i])
		{
			num[flag++]=i;//将素数存起来
			for(int j=i*i;j<Max;j+=i)
				prime[j]=true;//将有因子为i的合数筛掉
		}
}

但是上面的写法并不是最快的，最快的应该是欧拉线性筛，如下：

const int maxm = 1e7+10;
bool p[maxm];
int pri[maxm];
int init(int n) {
    for(int i = 2; i<=n; i++) {
        if(!p[i])
            pri[++pri[0]] = i;
        for(int j = 1; j<=pri[0]&&i*pri[j]<=n; j++) {
            p[i*pri[j]] = true;
            if(i%pri[j]==0)
                break;
        }
    }
    return pri[0];
}

仅为一个蒟蒻的个人之见，有任何问题欢迎交流，望大神勿喷。