《找素数》 题解

描述

给定区间[L, R] ， 请计算区间中素数的个数。

输入描述

两个数L和R。

输出描述

一行，区间中素数的个数。

用例输入 1 

2 11

用例输出 1 

5

提示

2 <= L <= R <= 2147483647 R-L <= 1000000

解题思路
本题用传统方法怎么都超时，看了题解才知道，需要用到埃氏筛法：这个方法用于寻找自然数n以内的所有素数，具体执行是剔除n以内所有的2到sqrt(n)之间的所有素数的倍数（不包含2、3、5等），剩下的即为素数。

易错点
注意变量是否需要long int，以防发生运行错误；
为了使hash散列的数组长度不超限，需要把[L,R]对应到[0,R-L]；
题目要求从L到R的素数，因此，只需要将在L和R之间的2到sqrt®的倍数删去即可；否则，从2开始查找会时间超限。

#include <bits/stdc++.h>
using namespace std;
const int maxn=1e6+7;
int m=(int)sqrt(2147483647+0.5)+10;
bool isp[maxn];
int p[maxn/10];
int len;
void init(){
    isp[0]=isp[1]=true;
    for(int i=2;i<m;i++){
        if(!isp[i])p[++len]=i;
        for(int j=1;j<=len&&p[j]*i<m;j++){
            isp[i*p[j]]=true;
            if(i%p[j]==0)break;
        }
    }
}

int main() {
    init();
    memset(isp,false,sizeof isp);
    int l,r;
    cin>>l>>r;
    int Len=r-l+1;
    for(int i=1;i<=len;i++){
        int L;
        if(p[i]>=l)L=p[i]-l;
        else{
            if(l%p[i]==0)L=0;
            else L=p[i]-l%p[i];
        }
        if(l+L==p[i])L+=p[i];
        for(int j=L;j<=Len;j+=p[i]){
            isp[j]=true;
        }
    }
    int cnt=0;
    for(int i=0;i<Len;i++){
        if(!isp[i])cnt++;
    }
    if(l==1)cnt--;
    cout<<cnt<<endl;
    return 0;
}