您的位置:首页 > 其它

codeforce 237 C. Primes on Interval

2018-01-21 12:09 501 查看
题意:

给出整数a,b 和数值k

a,b代表[a,b]区间,

求出最小的长度l,使得[a,b]区间内任何一个长度为l的连续区间中素数的数量>=k.

题目链接:链接传送门

思路:素数筛+二分区间长度

#include<cstdio>
#include<cmath>
#include<cstring>
#include<string>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<queue>
#include<stack>

using namespace std;

#define FOU(i,x,y) for(int i=x;i<=y;i++)
#define FOD(i,x,y) for(int i=x;i>=y;i--)
#define MEM(a,val) memset(a,val,sizeof(a))
#define PI acos(-1.0)

const double EXP = 1e-9;
typedef long long ll;
typedef unsigned long long ull;
const int INF = 0x3f3f3f3f;
const ll MINF = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9+7;
const int N = 1e6+5;

int prime
; //第i个素数是prime[i],从0开始记录
int vis
; //vis[i]代表当前合数被筛过没有,保证每个合数只被筛一次
bool is_prime
; //is_prime[i]为true代表i是素数

int euler(int n) //筛出小于等于n的所有素数并返回个数
{
int cnt=0; //记录个数
for(int i=0;i<=n;i++)
{
vis[i]=0;
is_prime[i]=false;
}
for(int i=2;i<=n;i++)
{
if(vis[i]==0)
{
prime[cnt++]=i;
is_prime[i]=true;
}
for(int j=0;j<cnt&&i*prime[j]<=n;j++)
{
vis[i*prime[j]]=1;
if(i%prime[j]==0)
break;
}
}
return cnt;
}

int sum
,a,b,k;

bool check(int x)
{
for(int i=a;i<=b-x+1;i++)
{
if(sum[i+x-1]-sum[i-1]<k)
{
//cout<<sum[i+x-1]<<" "<<sum[i-1]<<endl;
//cout<<"false"<<endl;
return true;
}
}
//cout<<"true"<<endl;
return false;
}

int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
std::ios::sync_with_stdio(false);
scanf("%d%d%d",&a,&b,&k);
int tmp=euler(b);
sum[0]=0;
sum[1]=0;
for(int i=2;i<=b;i++)
{
if(is_prime[i])
sum[i]=sum[i-1]+1;
else
sum[i]=sum[i-1];
//cout<<"sum["<<i<<"]="<<sum[i]<<endl;
}
int l=1;
int r=b-a+1;
if(sum[b]-sum[a-1]<k)
{
printf("-1\n");
return 0;
}
while(l<=r)
{
int mid = (l+r)/2;
if(check(mid))
l=mid+1;
else
r=mid-1;
//cout<<"l="<<l<<"r="<<r<<endl;
}
printf("%d\n",l);
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签:  codeforce 二分 素数筛