素数密度(求区间内素数个数)

题目传送门

题意:
给你一个区间,求这个区间内素数的个数。
数据范围
(L≤R≤2147483647,R-L≤1000000)

思路
区间是在int范围内的,所以一个数的最大质因子不会超过5e4。但是注意区间长度,如果我们遍历区间,然后一个一个判断,喜提TLE。但是一般的筛法又处理不到这个数据范围,所以我们就另辟蹊径。
前面说到区间内任何一个数的最大质因子不超过50000,那么我们可以先用筛法求出50000以内的素数,然后直接用这些素数对区间内的数进行筛,这样就不用处理无关数据了。

#include
using namespace std;
const int N=1e6+5;
const int inf=0x7fffffff;
const int mod=1e9+7;
const int eps=1e-6;
typedef long long ll;
typedef unsigned long long ull;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define pii pair
#define int long long
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define endl '\n'
bool ok[N];
int p[N];
int cnt=0;
bool isp[N];
void f()
{
    for(int i=2;i<=50000;i++)
    {
        if(!ok[i])
        {
            p[++cnt]=i;
            for(int j=i*i;j<=50000;j+=i)
            {
                ok[j]=1;
            }
        }
    }
}
signed main()
{
    IOS;
    int l,r;
    cin>>l>>r;
    f();
    for(int i=1;i<=cnt&&p[i]<=r;i++)
    {
        for(int j=l/p[i]*p[i];j<=r;j+=p[i])
        {
            if(j==p[i]||j<l)
            {
                continue;
            }
            isp[j-l]=1;
        }
    }
    int crt=0;
    for(int i=l;i<=r;i++)
    {
        if(isp[i-l]==0)
        {
            crt++;
        }
    }
    cout<<crt<<endl;
}

你可能感兴趣的:(素数筛)