埃氏筛法的原理
从最小的素数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;
}