埃氏筛法打印素数

埃氏筛法的原理

从最小的素数2开始，把所有2的倍数全划掉（除了2本身），再+1，把3的倍数全划掉，再+1,4已经被划掉了，就继续+1，划掉5的倍数。依次类推。
因为素数的定义是只能被1和自己本身整除，所以是倍数的必然不是。

一些优化

根据：所有合数（不是质数的数）n必定有一个小于sqrt（n）的质因子，所以我们在找倍数的外层循环小于等于sqrt（maxn）即可。
而我们在寻找倍数的时候从该数的平方开始即可。例如5这个素数，我们在划倍数的时候从5*5开始就行，因为2 * 5,3 * 5，4 * 5在前面都已经划掉了

代码

#include<stdio.h>
#include<math.h>
#define N 10000
int a[N];
void shy()
{
	a[1]=0;//去掉1这个非素数
	double s=N;
	for(int i=2;i<=sqrt(s);i++)
	{
		if(a[i])
		{
			for(int j=i*i;j<=N;j+=i)
			{
				a[j]=0;
			}
		}
	}
}
int main()
{
	for(int i=0;i<=N;i++)
		a[i]=1;
	shy();
	for(int i=2;i<=N;i++)
	{
		if(a[i])
			printf("%d\n",i);
	}
	return 0;
}