codeforce 237 C. Primes on Interval

题意:

给出整数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[N];  //第i个素数是prime[i],从0开始记录
int vis[N];  //vis[i]代表当前合数被筛过没有,保证每个合数只被筛一次
bool is_prime[N];  //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[N],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;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值