nyoj 187 素数快速筛法

#include <stdio.h>
#include <math.h>
#include <memory.h>
#include <stdlib.h>
long const M=2000010;
long a[M],p;
int main()
{
long i,j,x,y,k,q1,q2;
int t,n;
memset(a,1,sizeof(a));
for(i=2;i<=M;i++)
{
if(a[i])  // 如果已经被筛掉的 它的倍数肯定也已经被筛掉 	  （3筛掉 6 9 12 15。。。） 到了（12 12的倍数肯定是3的倍数 所以肯定也都被筛掉了）
{		            //j为倍数
for(j=2;i*j<=M;j++) //素数的倍数不是素数
{
a[i*j]=0;
}
}
}
while(scanf("%ld",&k)&&k)
{

for(i=2;i<=k;i++)
{
if(a[i])
{
printf("%d ",i);
}
}
printf("\n");

}
return 0;
}