筛选法求解M以内的素数

void prime(int m)
{
	int *a=(int *)malloc(sizeof(int)*(m+1));
	int i,j=2,k=2;
	for(i=0;i<=m;i++)
	{
		a[i]=i;
	}
	while(j<=m/2)
	{
		if(a[j]!=0)
		{
		   i=2*j;
		   while(i<=m)
		   {
				a[i]=0;
				i=i+j;
			}
		}
		j++;
	}
	for(i=2;i<=m;i++)
	{
		if(a[i]!=0)
		{
			printf("%d ",a[i]);
		}
	}
	free(a);
}

 具体过程为

一个简单的筛素数的过程：

n=30。     1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30        

 第 1 步过后2 4 ... 28 30这15个单元被标成false,其余为true。    

 第 2 步开始：      

i=3;  由于prime[3]=true, 把prime[6], [9], [12], [15], [18], [21], [24], [27], [30]标为false.      

i=4;  由于prime[4]=false,不在继续筛法步骤。     

 i=5;  由于prime[5]=true, 把prime[10],[15],[20],[25],[30]标为false.      i=6>sqrt(30)算法结束。     

第 3 步把prime[]值为true的下标输出来：      

for(i=2; i<=30; i++)      

if(prime) printf("%d ",i);    

 结果是 2 3 5 7 11 13 17 19 23 29