【NOIP2011模拟9.20】素数密度

Description

求区间[l~r]中素数的个数。
r-l+1<=10^6 l,r<2^31

Solution

呵呵了。刷出$\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[N];
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);
}