欧拉筛(线性筛)优化素数

普通筛法(埃氏筛)

我们一般要得到1~n以内的素数，用的都是埃拉托斯特尼筛法(后文都用埃氏筛 代替)。

那么，什么是埃氏筛?埃氏筛的本质从2开始删去素数本身倍数，向后找到的第一个数字一定是素数。
证明如下：
设已找到第n个素数，删去此数自身倍数后找到剩下的第一个数字L，知L之前有且仅有n个素数，且都无法整除L，即L无法被小于自身的所有素数整除，推出L是素数（L就是第n+1个素数）。照此由第一个素数2往后递推即可。

即从2开始，我们把素数的所有倍数(小于n)都剔除掉，再用下一个素数，也就是3筛，把3留下，把3的倍数剔除掉；接下去用下一个质数5筛，把5留下，把5的倍数剔除掉；不断重复下去…

#include <stdio.h>
int a[100000001]={0};
int main()
{
	int n,i,j,k=0;
	scanf("%d",&n);
	for(i=2;i<=n;i++)
	{
		if(a[i]==0)   //为0表示它是素数 
		{
			for(j=2;i*j<=n;j++)
			{
				a[i*j]=1;
			}
			k++;
		}
	}
	printf("共%d个素数\n",k);
	return 0;
}

但是埃氏筛有一个问题，有一些数会被多次筛掉，如6会先被2筛掉，又会被3筛掉，这里就有重复运算的问题，那么我们有没有一种一个数只被筛一次的筛法？这就是欧拉筛出现的原因。

欧拉筛(线性筛)

欧拉筛的定义是每个数只能被最小的质因数筛去，这样就能保证每个合数只被筛一次。

在实现形式上，埃氏筛是直接用素数遍乘自然数删去自身倍数；线性筛是先建立一个素数数列，然后用每个自然数遍乘数列中的素数删去素数倍数。
如：
第一步：建立标记数列（从2到n的自然数）和素数数列；

第二步：第一个数字是2（标记为素数），把素数2放进素数数列，标记2×2为非素数；

第三步：第二个数字是3（标记为素数），把素数3放入素数数列，标记3×2、3×3为非素数；

第四步：第三个数字是4（标记为合数），标记4×2为非素数，2为4的质因数，break；

第五步：第四个数字是5（标记为素数），把素数5放入素数数列，标记5×2、5×3、5×5为非素数；

第六步：第五个数字是6（标记为合数），标记6×2为非素数，2为6 的质因数，break；

第七步：第六个数字是7（标记为素数），把素数7放入素数数列，标记7×2、7×3、7×5、7×7为非素数；
…

//用线性筛法找到给定整数n以内的所有素数
#include<stdio.h>
int a[100000001]={0},p[50000000];
// a数组做自然数列表,p数组保存素数 
int main(){
	int n,i,j,k=0;
	scanf("%d",&n);
	for(i=2;i<=n;i++){
		if(a[i]==0){
			k++;
			p[k]=i;   //素数存入,注意是从1开始 
		}
        for(j=1;j<=k&&i*p[j]<=n;j++){// k保证不会引用没存储素数的p[j] ,i*p[j]设置运算的结束条件 
			a[i*p[j]]=1;               //将i乘以所有非最小质因数的素数的数标记为合数 
			if(i%p[j]==0){
			break;       //当i的最小质因数出现时，跳出循环 
		    }
		}
	}
	printf("共%d个素数\n",k);
	return 0;
}

时间复杂度大概是从O(nloglogn)缩减到O(n)