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

题目传送门

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

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

#include<bits/stdc++.h>
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<int,int>
#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;
}

发布了93 篇原创文章 · 获赞 9 · 访问量 4203

猜你喜欢

转载自blog.csdn.net/Joker_He/article/details/104884002