埃拉托斯特尼筛法

埃拉托斯特尼筛法

质数又称素数。指在一个大于1的自然数中，除了1和此整数自身外，没法被其他自然数整除的数。怎么判断n以内的哪些数是质数呢？

埃拉托斯特尼筛法

厄拉多塞是一位古希腊数学家，他在寻找素数时，采用了一种与众不同的方法：先将2－N的各数放入表中，然后在2的上面画一个圆圈，然后划去2的其他倍数；第一个既未画圈又没有被划去的数是3，将它画圈，再划去3的其他倍数；现在既未画圈又没有被划去的第一个数是5，将它画圈，并划去5的其他倍数……依次类推，一直到所有小于或等于N的各数都画了圈或划去为止。这时，表中画了圈的以及未划去的那些数正好就是小于 N的素数。
如下图所示：

而其实迭代系数i不需要遍历到n-1为止，只需到√(n-1)即可。反证法：

• 如果n-1不是质数，那么n-1可以化解成两个整数因子相乘，n-1=d1×d2。
• 如果d1和d2均大于√(n-1)，则有：n-1＝d1×d2 > √(n-1)×√(n-1)＝n-1
• 则n-1必有因子d1或d2小于等于√(n-1)，从而n-1可以被小于或等于√(n-1)的某整数的遍历到。
• 小于n-1的整数如果是质数必然会被遍历到。

如果N是合数，则一定存在大于1小于N的整数d1和d2，使得N=d1×d2。如果d1和d2均大于√N，则有：N＝d1×d2>√N×√N＝N。而这是不可能的，所以，d1和d2中必有一个小于或等于√N。

以100为例，我们先创建一个100个数字的数组。 
先使用最小的素数2，将所有2的倍数(除2本身)标记为合数。 
接下来2+1的数是3，此时检查3是不是素数，检查标记，发现没有被标记为合数(因为不是2的倍数)，所以再将所有3的倍数标记为合数。 
下一个数是4，发现他已经被标记为合数，所以他可以表示小于4的素数的乘积2*2，所以4的倍数必定含有因子2，所以所有4的倍数已经全部被标记过，直接跳过4。 
下一个数是5，没有被标记为合数，把所有小于100的5的倍数标记为合数 
………这样一直计算到sqrt(100)，即10。 
那么为什么不标记大于10的数例如11呢？因为所有的倍数已经被标记过了，例如22，33，44，55…分别有因子2，3，2，5，大于10的倍数，例如11*11已经超过max了，参见最上面的推论 
注意这里有一个优化点，很多书籍上或者教程上都没有说出来，只要标记大于本身的倍数就行了，例如5，只要标记5*5,5*6,5*7…为合数，因为5*2,5*3,5*4…已经被之前出现的数的倍数标记过了

代码实现

const int maxn=1e3;
int a[maxn];
void isprime(int n)//埃拉托斯特尼筛法
{
	fill(a,a+n,1);//将所有的数都认为是素数初始化为 1， 
    a[0]=a[1]=0//0和1显然不是素数，标记为0
	for(int i=2;i*i<=n;i++)
	{
		if(a[i])//如果i是素数，显然 i 的倍数都不是素数， 
		{
			for(int j=i*i;j<=n;j+=i)这里直接忽略了小于i的倍数，因为之前肯定已经出现过了
				a[j]=0;//不是素数标记为0  
		}
	}
}