描述
给定区间[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;
}