基础数学——素数密度

素数密度

题目描述

给定区间 $[L,R]$（$1\le L\le R<2^{31}$，$R-L\le 10^6$），请计算区间中素数的个数。

输入格式

第一行，两个正整数 $L$ 和 $R$。

输出格式

一行，一个整数，表示区间中素数的个数。

样例 #1

样例输入 #1

2 11

样例输出 #1

5

思路分析

本道题为区间素数问题，又因为时间复杂度$R-L<=10^6$，所以我们得考虑二次筛法，细节：二次筛法适用于区间求素数的个数，相邻或相近问题，比如质数距离，而且这里的区间范围给定的最大值是$2^{31}-1$,而用线性筛法求的是[1,n]中的所有质数，因此直接用线性筛法求肯定会直接TLE，这里先说几个性质：

• 性质1：若一个数n是一个合数，必然存在2个因子d,n/d，假设d <= nd，则d<= √n，因此必然存在一个小于等于√n 的因子
• 性质2：若x∈[L,R],且x是合数，则一定存在P <=√2^31-1(< 50000)，使得P能整除x，其中P < x

性质来源，这里讲的更详细

代码演示

//因为只要求求素数个数，个人感觉可以

#include<iostream>
#include<algorithm>
#include<cstring>

using namespace std;

const int N = 1e6+10;

typedef long long LL;

int primes[N],cnt;
bool st[N];
int l,r;

void init(int x){
    for(int i=2;i<=x;i++){
        if(!st[i])primes[cnt++]=i;
        for(int j=0;i*primes[j]<=x;j++){
            st[i*primes[j]]=true;
            if(i%primes[j]==0)break;
        }
    }
}

int main(){
    cin>>l>>r;
    init(50000);
    l+=(l==1);
    // int res=0;
    // for(int i=l;i<=r;i++){
    //     if(!st[i])res++;
    // }
    
    // cout<<res;
    
    memset(st,0,sizeof st);
    LL res=0;
    for(int i=0;i<cnt;i++){
        LL p=primes[i];
        
        for(LL j=max(p*2,(p+l-1)/p*p);j<=r;j+=p){
            st[j-l]=true;
        }
    }
    
    for(int i=0;i<=r-l;i++){
        if(!st[i])res++;
    }
    
    
    cout<<res;
    
    return 0;
}