素数筛法C++代码

这种素数筛法时间复杂度大约为O(n*sqrt(n)),略差于最优线性筛法的O(n)，但是面对一般的OJ题目足够用了，而且简单易记。

使用素数筛法，输出0~Max之间的所有素数

#include <iostream>

using namespace std;

#define Max 100         //0~max_num之间的素数
int prime[Max+1];       //存储筛选出来的素数

//生成0~Max之间的素数，存放在prime数组中
int generate_prime()
{
    int p[Max+1];       //标记数组
    int cnt=0;          //素数的个数
    for(int i=0;i<=Max;i++)
        p[i]=1;
    p[0]=p[1]=0;        //0,1不是素数
    for(int i=2;i<=Max;i++)
    {
        if(p[i]==1)     //如果p[i]=1，则i是素数
        {
            prime[cnt++]=i;             //将素数i存入prime数组中，cnt加1
            for(int j=i;j<=Max;j=j+i)   //将素数i的倍数都标记为非素数
                p[j]=0;
        }
    }
    return cnt;//返回素数的个数
}

int main()
{
    int cnt=generate_prime();   //生成素数
    for(int i=0;i<cnt;i++)      //输出素数
        cout<<"  " << prime[i];
    cout<<endl;
    return 0;
}

将Max设置为100，生成0~100之间的所有素数，运行结果如下：