区间素数筛

Posted 2020-10-14

有些题目让你求一个区间的素数

基于普通的埃氏筛法， 衍生出区间筛法，例如求 [ a, b ]内的素数，可以先将 [ 0, 根号b ]内的素数全都找出来，然后在将 [ a, b ]内是其 [ 0, 根号b ] 中素数的倍数的数全部划去，区间剩余的数即为素数。

bool prime[47000];
bool pt[eps]; 

memset(pt, true, sizeof(pt));
        memset(prime, true, sizeof(prime));
        
        for(ll i = 2; i*i <= b; i++){
            if (prime[i]){
                for(ll j = i+i; j*j <= b; j += i){
                    prime[j] = false;
                }
                for(ll j = i*max(2LL, (a+i-1)/i); j <= b; j += i){ // 区间要删除的数至少是素数的2倍，因此再求得倍数的时候，要与2求一个max 
                    pt[j] = false;
                }
            }
        }