【NOIP2011模拟9.20】素数密度

Description

求区间[l~r]中素数的个数。

r-l+1<=10^6 l,r<2^31

Solution

呵呵了。刷出maxl−−−−√\sqrt{maxl}中所有的素数，然后再用这一坨素数去刷l~r这个区间。

没有了，大水题。

为什么我开50000过不了！

逼我开100000！

Code

#include<cstdio>
#include<cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define N 1000005
#define M 100000
using namespace std;
bool f[M+5],g
;
int l,r,ans;
int main() {
f[1]=1;
fo(i,2,sqrt(M))
if (!f[i])
fo(j,2,M/i) f[i*j]=1;
scanf("%d%d",&l,&r);
if (r<=M) {
fo(i,l,r) if (!f[i]) ans++;
printf("%d",ans);return 0;
}
if (l<=M) {
fo(i,l,M) if (!f[i]) ans++;
l=M+1;
}
fo(i,1,M)
if (!f[i]) fo(j,l/i,r/i) g[i*j-l+1]=1;
fo(i,1,r-l+1) if (!g[i]) ans++;
printf("%d",ans);
}