素数筛法模板

代码:

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <math.h>
#include <queue>
#define MAX 1000001
using namespace std;
bool b[MAX];
int main()
{
    b[0]=b[1]=false;
    b[2]=true;
    for(int i=3; i<MAX; i++)
        if(i%2==0) b[i]=false;
        else b[i]=true;
    double t=sqrt(1000000*1.0);
    for(int i=3; i<=t; i++)
    {
        if(b[i])
        {
            for(int j=i*i; j<MAX; j=j+i)//j=i*i;如:j=5*5(因为2*5,3*5,4*5在之前已经被筛去,节约时间),因为i<=sqrt(MAX-1),所以i*i不会超出int范围。
            {
                b[j]=false;
            }
        }
    }
    for(int i=1; i<=100; i++)
    {
        if(b[i]) printf("%d ",i);
    }
}

 

posted @ 2014-11-26 17:18  人艰不拆_zmc  阅读(187)  评论(0编辑  收藏  举报