用C语言实现素数筛法获取一亿（100000000）以内的全部素数

具体筛法是：先把n个自然数按次序排列起来。1不是质数，也不是合数，要划去。第二个数2是质数留下来，而把2后面所有能被2整除的数都划去。2后面第一个没划去的数是3，把3留下，再把3后面所有能被3整除的数都划去。3后面第一个没划去的数是5，把5留下，再把5后面所有能被5整除的数都划去。这样一直做下去，就会把不超过N的全部合数都筛掉，留下的就是不超过N的全部质数。因为希腊人是把数写在涂腊的板上，每要划去一个数，就在上面记以小点，寻求质数的工作完毕后，这许多小点就像一个筛子，所以就把埃拉托斯特尼的方法叫做“埃拉托斯特尼筛法”，简称“筛法”。

下面是代码：

 1 #include <stdio.h>

 2 #include <stdlib.h>

 3 #include <string.h>

 4 

 5 int main(int argc, char * argv[])

 6 {
　　　　//寻找2~num之间的所有素数

 7     if(argc < 2)

 8     {

 9         printf("Usage : %s num\n", argv[0]);

10         return 0;

11     }

12     int iMax = atoi(argv[1]);

13 

14     if(iMax < 2)

15     {

16         printf("num is too little, num >=2");

17         return 0;

18     }

19 

20     char *p = (char *)malloc(sizeof(char) * iMax + 1);

21     memset(p, sizeof(char) * iMax + 1, 0);

22 

23     int i = 0, j = 0, k = 0;

24     for(i = 2; i <= iMax; i++)

25     {

26         for(j = i + i; j <= iMax; j += i)

27         {

28             p[j] = 1;

29         }

30     }

31     FILE * fp = NULL;
　　　　//程序执行完成后，文件 prime-number.txt中就是我们需要的素数

32     if((fp = fopen("prime-number.txt", "w")) == NULL)

33     {

34         return 0;

35     }

36     k = 0;

37     int iAll = 0;

38     for(i = 2; i <= iMax; i++)

39     {

40         if(0 == p[i])

41         {

42             iAll ++;

43             k++;

44             // output to file : fp，把这些素数写入文件

45             fprintf(fp, "%6d ", i);

46             if(10 == k)

47             {

48                 fprintf(fp, "\n");

49                 k = 0;

50             }

51             //printf("%d ", i);

52         }

53     }

54     printf("\n");

55     fclose(fp);

56     free(p);

57     printf("all : %d\n", iAll);

58 

59     return 0;

60 }

 

输出结果放在百度网盘 ：http://pan.baidu.com/s/1pJv58Wb

 

作者：风波

mail : [email protected]