nyoj 187 素数快速筛法

http://acm.nyist.net/JudgeOnline/problem.php?pid=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;
}