洛谷 P1835 素数密度

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
int prime[N],cnt,ans;
bool v[N];
void Prime(int x){
	for(int i=2;i<x;i++){
		if(v[i]){
			continue;
		}
		prime[++cnt]=i;
		for(int j=i;j<x/i;j++){
			v[i*j]=1;
		}
	}
}
int main(){
	long long l,r;
	scanf("%11d %11d",&l,&r);
	if(l==1){
		l=2;
	}
	Prime(50000);
	memset(v,0,sizeof(v));
	for(int i=1;i<=cnt;++i){
		long long p=prime[i],start;
		if((l+p-1)/p*p>2*p){
			start=(l+p-1)/p*p;
		}
		else{
			start=2*p;
		}
		for(long long j=start;j<=r;j+=p){
			v[j-l+1]=1;
		}
	}
	for(int i=1;i<=r-l+1;i++){
		if(!v[i]){
			ans++;
		}
	}
	printf("%d",ans);
	return 0;
}