【洛谷】P1835 素数密度

洛谷P1835 素数密度

题目描述
给定区间[ L，R ] (L≤R≤2147483647,R-L≤1000000)，请计算区间中素数的个数。

输入格式
两个数L和R。

输出格式
一行，区间中素数的个数。

输入输出样例
输入#1:

2 11

输出#1:

5

这里是要统计一段范围内的素数个数。而这段范围可以在int的任意范围内，只不过区间的大小小于1000000。显然之前线性筛法已经无法直接实现，而朴素的$O(n\sqrt{n})$判断应该是也会超时的，所以一时没想到什么思路。
但其实仔细想想，所有int范围内的合数都一定存在一个$\sqrt{IN{T}_{M}AX}$以下的质因数。所以我们可以先用线性筛把这个范围内的质因数求出来，然后再用这些质因数标记[L,R]范围内的合数即可。

• prime是素数列表，isprime标记是否是素数
• build()线性筛先预处理素数
• 对于素数prime[i]，要得到第一个大于等于L的倍数，令j=(L-1)/prime[i]+1的技巧挺不错的
• 但是还要注意j至少为2才有意义，否则j=0或者1的时候会出错
• 另外要注意(long)j*prime[i]<=(long)r，这个long很关键。因为如果都是int会溢出

代码如下：

#include <bits/stdc++.h>
using namespace std;
int l,r,now=1;
int prime[50000],isprime[50010];
void build(){
    isprime[1]=0;
    for(int i=2;i<50000;i++){
        if(isprime[i]){
            prime[now++]=i;
        }
        for(int j=1;j<=now&&i*prime[j]<50000;j++){
            isprime[i*prime[j]]=0;
            if(i%prime[j]==0) break;
        }
    }
}
int res[1000010],ans=0;
int main(){
    cin>>l>>r;
    memset(isprime,1,sizeof(isprime));
    build();
    memset(res,1,sizeof(res));
    for(int i=1;i<now&&prime[i]<=r;i++){
        for(long j=max(2,(l-1)/prime[i]+1);j*prime[i]<=(long)r;j++){
            res[j*prime[i]-l]=0;
        }
    }
    if(l<=1) res[1-l]=0;
    for(int i=0;i<=r-l;i++){
        if(res[i]) ans++;
    }
    cout<<ans;
    return 0;
}