区间素数筛模版

区间素数筛模版

筛出区间[a,b]的素数。(b-a<=10000,1<=a<=b<=2^31)

存在P中,素数个数即为P的size()。

ll a,b;
bool isprime[maxn];
vector<ll> prime;
bool isP[maxn];
vector<ll> P;

void play_prime()
{
    memset(isprime,1,sizeof(isprime));
    isprime[1]=0;
    for(int i=2;i<maxn;i++){
        if(!isprime[i]) continue;
        for(int j=i*2;j<maxn;j+=i){
            isprime[j]=0;
        }
    }
    for(int i=1;i<maxn;i++){
        if(isprime[i]) prime.push_back(i);
    }
}

void get_prime(ll a,ll b)
{
    if(a<=1) a=2;
    memset(isP,1,sizeof(isP));
    for(int i=0;i<prime.size(),prime[i]*prime[i]<=b;i++){
        ll p=prime[i];
        for(ll j=a+(p-a%p)%p;j<=b;j+=p){
            if(j!=p) isP[j-a]=0;
        }
    }
    P.clear();
    for(ll i=a;i<=b;i++){
        if(isP[i-a]) P.push_back(i);
    }
}
View Code

 

posted @ 2015-06-10 13:35  __560  阅读(207)  评论(0编辑  收藏  举报