现在的位置: 首页 > 综合 > 正文

//sieve筛选法求一段区间内的素数

2013年12月07日 ⁄ 综合 ⁄ 共 467字 ⁄ 字号 评论关闭

#define maxn 1024000

int L,U;
int p[maxn],sz;//[0,sz)
bool flag[maxn];//U-L<=maxn
void sieve(int L,int U){
     int d=U-L+1;
     for(int i=0;i<d;i++)flag[i]=true;
     for(int i=(L%2!=0);i<d;i+=2)flag[i]=false;
     for(int i=3;i<=sqrt(U+0.0);i+=2){
        if(i>L&&flag[i-L]==false)continue;
        int j=L/i*i;
        if(j<L)j+=i;
        if(j==i)j+=i;
        j=j-L;
        for(;j<d;j+=i)flag[j]=false;
     }
     if(L<=1)flag[1-L]=false;
     if(L<=2)flag[2-L]=true;
     sz=0;
     for(int i=0;i<d;i++)if(flag[i]==true)p[sz++]=i+L;
}

【上篇】
【下篇】

抱歉!评论已关闭.