用C語言實現素數篩法獲取一億（100000000）以內的全部素數

具體篩法是：先把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 < 2)
    {
        printf("Usage : %s num\n", argv[0]);
        return 0;
    }
    int iMax = atoi(argv[1]);

    if(iMax < 2)
    {
        printf("num is too little, num >=2");
        return 0;
    }

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

    int i = 0, j = 0, k = 0;
    for(i = 2; i <= iMax; i++)
    {
        for(j = i + i; j <= iMax; j += i)
        {
            p[j] = 1;
        }
    }
    FILE * fp = NULL;
　　　　//程序執行完成后，文件 prime-number.txt中就是我們需要的素數
    if((fp = fopen("prime-number.txt", "w")) == NULL)
    {
        return 0;
    }
    k = 0;
    int iAll = 0;
    for(i = 2; i <= iMax; i++)
    {
        if(0 == p[i])
        {
            iAll ++;
            k++;
            // output to file : fp，把這些素數寫入文件
            fprintf(fp, "%6d ", i);
            if(10 == k)
            {
                fprintf(fp, "\n");
                k = 0;
            }
            //printf("%d ", i);
        }
    }
    printf("\n");
    fclose(fp);
    free(p);
    printf("all : %d\n", iAll);

    return 0;
}

 

輸出結果放在百度網盤 ：http://pan.baidu.com/s/1pJv58Wb

 

作者：風波

mail : [e-mail]