[大数区间素数]洛谷P1835素数密度

题目

题目

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<cctype>
#include<ctime>
#include<iostream>
#include<string>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<vector>
#include<iomanip>
#include<list>
#include<bitset>
#include<sstream>
#include<fstream>
#include<complex>
#include<algorithm>
#if __cplusplus >= 201103L
#include <unordered_map>
#include <unordered_set>
#endif
#define int long long
using namespace std;
const int INF = 0x3f3f3f3f;
int prime[500010];
bool vis[500010];
int cnt=0;
void getprime(int n){
	for(int i=2;i<=n;i++){
		if(!vis[i]) prime[++cnt]=i;
		for(int j=1;j<=cnt&&prime[j]*i<=n;j++){
			vis[i*prime[j]]=1;
			if(i%prime[j]==0) break;
		}
	}
}
signed main(){
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	getprime(50000);
	int r,l;
	cin>>l>>r;
	unordered_map<int,int> mp;
	int ans=0;
	for(int i=1;i<=cnt;i++){
		for(int j=max(2ll,(l/prime[i]));j*prime[i]<=r;j++){
			if(j*prime[i]<l) continue;
			if(!mp[j*prime[i]]) ans++;
		//	cout<<j*prime[i]<<" "<<ans<<endl;
			mp[j*prime[i]]=1;
		}
	}
	if(l==1) ans++; 
	cout<<r-l+1-ans<<endl;
    return 0;
}