埃拉托色尼筛法(来自:算法:C语言实现)

这个程序的功能是: 如果自然数i是素数,则设a[i]为1,否则设为0.首先把数组中的所有元素设为1,否则设为0.首先把数组中的所有元素设为1,已表明没有任何数以被证明是非素数.然后,把数组中所对应索引处已证明是非素数(已知素数的倍数)的元素设为0.如果所有更小素数的倍数都已设为0,a[i]仍为1,则可知它是素数.

因为程序中所用的数组由最简单的元素类型(0-1值)组成的数组,所以直接使用由位组成的数组比整数组成的更省空间.另外,如果N值过大,某些编程环境可能要求把数组定义为全局变量,或者我们可以动态地为它分配空间

#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <malloc.h>
using namespace std;
#define OK 1
#define ERROR -1
#define TRUE 1
#define FALSE 0
#define N 1000
typedef int Status;

int main()
{
    //定义两个变量i,j,一个数组a[N]
    int i, j, a[N];
    //将数组初始化为1,a[i] = 1表示i为素数
    //初始时默认从2到N-1均为素数
    for(i = 2; i < N; i++ )
    {
        a[i] = 1;
    }
    //遍历数组找出下标为素数的,并将所有下标为该素数的倍数的值改为0
    for(i = 2; i < N; i++)
    {
        if(a[i])
        {
            for(j = i; i*j < N; j++)
            {
                a[i*j] = 0;
            }
        }
    }
    //输出2~N范围内的素数
    for(i = 2; i < N; i++)
    {
        if(a[i])
        {
            printf("%4d", i);
        }
    }
    return 0;
}