poj2689(Prime Distance两次筛)

两次筛,这里第一次用的是欧拉筛,第二次用的是埃氏筛(似乎第二次只能用埃氏筛)。
for(int i=0;i<=R-L;i++)
if(!isnPres[i]) res[l++]=i+L;
/*
for(int i=L;i<=R;i++)//如果这样写的话就runtime error
if(!isnPres[i-L]) res[l++]=i;
*/
那个好心人能告诉我这两者的差别!

#include <iostream>
#include <algorithm>
#include <string.h>
#include <math.h>
#include <stdio.h>
#define Inf 0x7fffffff
#define ms(X) memset(X,0,sizeof(X))
bool isnPres[1008020];
int t=0,mPrime[101000],res[1008020];
void create(void)
{
    ms(isnPres);
    for(int i=2;i<101000;i++)
    {
        if(!isnPres[i]) mPrime[t++]=i;
        for(int j=0;mPrime[j]*i<101000;j++)
        {
            isnPres[mPrime[j]*i]=true;
            if(i%mPrime[j]==0) break;
        }
    }
}
void createres(int st,int ed)
{
    ms(isnPres);
    for(int i=0;i<t&&(long long )mPrime[i]*mPrime[i]<=ed;i++)
    {
        int s=st/mPrime[i]+(st%mPrime[i]>0);
        if(s==1) s=2;
        for(int j=s;(long long)j*mPrime[i]<=ed;j++)
            if((long long)j*mPrime[i]>=st) isnPres[j*mPrime[i]-st]=true;
    }
}
int main(int argc, char const *argv[])
{
    create();
    int L,R,l;
    while(scanf("%d %d",&L,&R)!=EOF)
    {
        l=0;
        ms(res);
        if(L>2) L|=1;
        if(L<2) L=2;
        createres(L,R);
        for(int i=0;i<=R-L;i++)//这里想不通
            if(!isnPres[i]) res[l++]=i+L;
            /*
            for(int i=L;i<=R;i++)//如果这样写的话就runtime error
            if(!isnPres[i-L]) res[l++]=i;
            */
        if(l<2) {printf("There are no adjacent primes.\n");continue;}
        int s=res[1]-res[0],m=res[1]-res[0],s1=res[0],s2=res[1],m1=res[0],m2=res[1];
        for(int j=1;j<l-1;j++)
        {
            if(res[j+1]-res[j]>m) m=res[j+1]-res[j],m1=res[j],m2=res[j+1];
            if(res[j+1]-res[j]<s) s=res[j+1]-res[j],s1=res[j],s2=res[j+1];
        }
        printf("%d,%d are closest, %d,%d are most distant.\n",s1,s2,m1,m2);
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值