Article Directory
Elbow Sieve
Find the number of prime numbers less than the integer n.
int prime[MAX_N];// 存储素数
bool is_prime[MAX_N];//表示第i个数是否为素数
int sieve(int n){
int p=0;
for(int i=0;i<=n;i++) is_prime[i]=true;
is_prime[0]=is_prime[1]=false;
for(int i=2;i<=n;i++){
if(is_prime[i]){
prime[p++]=i;
for(int j=2*i;j<=n;j+=i) is_prime[j]=false;
}
}
return p;
}
Interval sieve
Find the number of prime numbers in the left open and right closed interval [a, b).
Idea: First find [2, b \sqrt bb) In the prime numbers, and then draw these numbers from the target interval.
typedef long long ll;
bool is_prime[MAX_L];
bool is_prime_small[MAX_SQRT_B];
void segment_sieve(ll a,ll b){
//初始化
for(int i=0;(ll)i*i<b;i++) is_prime_small[i]=true;
for(int i=0;i<b-a;i++) is_prime[i]=true;
//枚举小区间里的数
for(int i=2;(ll)i*i<b;i++){
if(is_prime_small[i]){
//筛小区间
for(int j=2*i;(ll) j*j<b;j+=i) is_prime_small[j]=false;
//筛大区间
for(ll j=max(2LL,(a+i-1)/i)*i;j<b;j+=i) is_prime[j-a]=false;
}
}
}