具体筛法是：先把n个自然数按次序排列起来。1不是质数，也不是合数，要划去。第二个数2是质数留下来，而把2后面所有能被2整除的数都划去。2后面第一个没划去的数是3，把3留下，再把3后面所有能被3整除的数都划去。3后面第一个没划去的数是5，把5留下，再把5后面所有能被5整除的数都划去。这样一直做下去，就会把不超过N的全部合数都筛掉，留下的就是不超过N的全部质数。因为希腊人是把数写在涂腊的板上，每要划去一个数，就在上面记以小点，寻求质数的工作完毕后，这许多小点就像一个筛子，所以就把埃拉托斯特尼的方法叫做“埃拉托斯特尼筛法”，简称“筛法”。
下面是代码：
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <strings.h>
int main(int argc, char * argv[])
{
//寻找2~num之间的所有素数
if(argc < )
{
printf("Usage : %s num\n", argv[]);
return ;
}
int iMax = atoi(argv[]); if(iMax < )
{
printf("num is too little, num >=2");
return ;
} char *p = (char *)malloc(sizeof(char) * iMax + );
bzero(p, sizeof(char) * iMax + ); int i = , j = , k = ;
for(i = ; i <= iMax; i++)
{
for(j = i + i; j <= iMax; j += i)
{
p[j] = ;
}
}
FILE * fp = NULL;
//程序执行完成后，文件 prime-number.txt中就是我们需要的素数
if((fp = fopen("prime-number.txt", "w")) == NULL)
{
return ;
}
k = ;
int iAll = ;
for(i = ; i <= iMax; i++)
{
if( == p[i])
{
iAll ++;
k++;
// output to file : fp，把这些素数写入文件
fprintf(fp, "%6d ", i);
if( == k)
{
fprintf(fp, "\n");
k = ;
}
//printf("%d ", i);
}
}
printf("\n");
fclose(fp);
free(p);
printf("all : %d\n", iAll); return ;
}
输出结果放在百度网盘 ：http://pan.baidu.com/s/1pJv58Wb
作者：风波
mail : [email protected]