数论(筛区间质数) - Help Hanzo LightOJ - 1197
题意:
求 区 间 [ a , b ] 有 多 少 个 素 数 求区间[a,b]有多少个素数 求区间[a,b]有多少个素数
Input
多 组 输 入 , 第 一 行 是 t 表 示 有 t 组 数 据 , 接 下 来 每 行 一 对 a , b 多组输入,第一行是t表示有t组数据,接下来每行一对a,b 多组输入,第一行是t表示有t组数据,接下来每行一对a,b
Output
输 出 " C a s e x : " 后 再 输 出 答 案 输出"Case \ x:"后再输出答案 输出"Case x:"后再输出答案
Sample Input
3
2 36
3 73
3 11
Sample Output
Case 1: 11
Case 2: 20
Case 3: 4
Note
t ≤ 200 , 1 ≤ a ≤ b < 2 31 , b − a ≤ 100000 t ≤ 200 ,1 ≤ a ≤ b < 2^{31},b - a ≤ 100000 t≤200,1≤a≤b<231,b−a≤100000
T i m e l i m i t : 2000 m s , M e m o r y l i m i t : 32768 k B Time\ limit:2000 ms,Memory\ limit:32768 kB Time limit:2000ms,Memory limit:32768kB
分析:
先 用 线 性 筛 法 , 预 处 理 出 b m a x 内 的 质 数 , 先用线性筛法,预处理出\sqrt{b_{max}}内的质数, 先用线性筛法,预处理出bmax内的质数,
接 着 通 过 埃 式 筛 法 , 利 用 预 处 理 出 的 质 数 , 筛 去 区 间 [ a , b ] 内 的 合 数 。 接着通过埃式筛法,利用预处理出的质数,筛去区间[a,b]内的合数。 接着通过埃式筛法,利用预处理出的质数,筛去区间[a,b]内的合数。
标 记 时 , 将 区 间 [ a , b ] 映 射 到 [ 0 , b − a ] 。 标记时,将区间[a,b]映射到[0,b-a]。 标记时,将区间[a,b]映射到[0,b−a]。
参考:线性筛 + 埃式筛 (筛区间质数) - Prime Distance - POJ 2689
代码:
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
const int N=1e5+10;
int primes[N],cnt;
bool st[N];
void get_prime(int n)
{
for(int i=2;i<=n;i++)
{
if(!st[i]) primes[cnt++]=i;
for(int j=0;primes[j]*i<=n;j++)
{
st[primes[j]*i]=true;
if(i%primes[j]==0) break;
}
}
}
int cal(int l,int r)
{
memset(st,false,sizeof st);
int res=0;
for(int i=0;i<cnt;i++)
{
ll p=primes[i];
ll start=max(2*p,(l+p-1)/p*p);
for(ll j=start;j<=r;j+=p)
st[j-l]=true;
}
for(int i=0;i<=r-l;i++)
if(!st[i]&&i+l>=2)
res++;
return res;
}
int main()
{
get_prime(N-1);
int T; cin>>T;
for(int t=1;t<=T;t++)
{
int l,r; cin>>l>>r;
printf("Case %d: %d\n",t,cal(l,r));
}
return 0;
}