素数筛选：普通，埃氏筛，欧拉筛

如何找到1-N间所有质数。介绍三种方法，从易到难，从慢到快

1.普通筛法，就是直接用质数的定义来做

代码：

#include<bits/stdc++.h>
using namespace std;
int main()
{
	int i,j,n;
	bool p=true;
	cin>>n;
	for(i=2;i<=n;i++)
	{
		p=1;
		for(j=2;j*j<=i;j++)
		{
			if(i%j==0)//能被其它数整除，说明不是质数 
			p=0;
		}
		if(p)
		{
			cout<<i<<" ";
		}
	}
}

2.埃氏筛

上一个代码的思路是找到一个数，再去检验它。

那我们能不能主动出击，主动找到所有合数？

于是我们可以将所有已知质数的倍数都给删去，剩下的最小数一定是质数。

再不断推进，最终筛去范围内所有质数。

#include<bits/stdc++.h>
using namespace std;
int panduan[1000005];
void ai_shi(int n)
{
	int i,j;
	for(i=2;i<=n;i++)
	{
		if(!panduan[i])//用质数 
		{
			for(j=2;j*i<=n;j++)
			{
				panduan[i*j]=1;//去掉合数 
		    } 
		}
	}
}
int main()
{
	int i,j,n;
	cin>>n;
	memset(panduan,0,sizeof(panduan));
	ai_shi(n);
	for(i=2;i<=n;i++)
	{
		if(!panduan[i])
		{
			cout<<i<<" ";
		}
	}
 } 

以上两种方法都比较简单，其实大部分题目用埃氏筛足以解决。但是丧心病狂的洛谷怎么可能这么简单

于是就有第三种筛法：欧拉筛

就是埃氏筛的升级。

在埃氏筛中，假如我们要筛去数字15

则上一代码的子函数ai_shi（） 中，

当i=3，j=5时，15被筛了一次

当i=5，j=3时，15又被筛了一次

重复的一次浪费了时间，所以我们能不能保证，每个数都被它的最小质因数筛去呢？

欧拉筛比较复杂，先上代码

#include<bits/stdc++.h>
using namespace std;
int j=0;//记录质数个数 
int prime[10000005];//建立一个质数集 
int pan_duan[1000005];//这个集合的用处与埃氏筛无异 
void ou_la(int n)
{
	int i,p;
	for(i=2;i<=n;i++)
	{
		if(!pan_duan[i])//没被筛去的是质数 
		{
			prime[j++]=i;
		}
		for(p=0;i*prime[p]<=n;p++)/* i*pan_duan[p]<n 是一个小优化*/
		{
			pan_duan[i*prime[p]]=1;
			
			if(i%prime[p]==0)
			{
				break;
			}
		}
	}
}
int main()
{
	int n,i;
	cin>>n;
	memset(pan_duan,0,sizeof(pan_duan));
	ou_la(n);
	for(i=0;i<j;i++)
	{
		cout<<prime[i]<<" ";
	}
}

整个代码的精髓在于第15到第24行

解决两个问题：

1.第17行中，为什么要用每个已知的质数依次与i想乘来找出合数。

2.第19行（划重点），为什么 i%prime[p]==0 就要跳出循环。

问题一

欧拉筛所谓的用最小质因数来筛，并不是指找到一个质数，然后再找出他的倍数，不然与埃氏筛无异。

我们要找的，是与最小质因数相对的那个数，用它来乘以每个质数

比如18=2*9=3*6=9*2。

那么，我们要做的就是找到i=9时，再拿i与质数集中的2相乘，这样筛去的18。

对比欧拉筛，是先找到2，再依次筛去2*2,2*3,2*4...直到2*9筛去了18。

而欧拉筛之所以要选择这种方法，就是因为问题二的答案

问题二：

i%prime[p]有两种情况

1.！=0

此时说明从2开始到质数prime[p]都不是i的最小质因数，那么i*prime[p]的最小质因数就一定是prime[p]这个质数。

2.==0

此时，我们可以确定prime[p]是i的最小质因数。那么i*prime[p]的最小质因数也是i

但是i*prime[p+1]的最小质因数就不再是prime[p+1]了

再用18作为例子

当i=18时

pan_duan[18*prime[1]]=1;   即：18*2=36不是质数且最小质因数是2

p++之后

pan_duan[18*prime[2]]=1    即：18*3=54不是质数且最小质因数是3

54的最小质因数是2。

第二句话明显是错误的。

筛去54的应该是当i=27时，用2*27=54才对。

因此直接用 break 跳出循环，避免重复找数。

———————————————————————————————————————————

说了那么一堆，其实就是保证每个数都被其最小质因数筛去

54=2*3*3*3，那么它只能被2筛去，不能被三筛去。

所以在欧拉筛中 54 只能等于 2*27。

解决了这两个问题，欧拉筛能简单理解了。