POJ2689 Prime Distance（素数、线性筛、试除法）

最新推荐文章于 2021-05-22 22:05:17 发布

shawwwOvO

最新推荐文章于 2021-05-22 22:05:17 发布

阅读量163

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/weixin_44986601/article/details/105290590

版权

数论专栏收录该内容

13 篇文章 0 订阅

订阅专栏

POJ2689 Prime Distance
坑死我了，除0导致re到si。还一直以为数组越界自闭了一早上。

题意：给定两个整数L、U。求闭区间中相邻两素数的差的最大值和最小值，输出这两个质数。LU的范围2e9。U-L范围1e6。

思路：
LU的范围2e9，直接线性筛子会 MLE&&TLE。但是U-L范围是1e6。
任何合数n必定包含一个不超过 $\sqrt{n}$ 的质因子。
所以只要用线性筛找出 $\sqrt{U}$ 范围内的所有素数，然后对每一个素数在[L，U]范围内的合数打标记mp，最后遍历标记数组mp，找到minn和maxx即可。

LL v[maxn],prime[maxn];//v存质数，vis判断是不是质数
bool mp[maxn];
int primes(LL n){
    int m=0;
    for(LL i=2;i<=n;i++){
        if(v[i]==0){//i是质数
            v[i]=i;prime[++m]=i;
        }
        for(int j=1;j<=m;j++){
            if(prime[j]>v[i]||prime[j]>n/i)break;
            v[i*prime[j]]=prime[j];
        }
    }
    return m;
}
int main(){
    LL l,u,minn,maxx,c1,c2,d1,d2,pre,sum;
    int co=primes(sqrt(2147483647)),tt;
    while(scanf("%lld%lld",&l,&u)!=EOF){
        minn=1e9+10;maxx=-1;
        sum=u-l+1;
        memset(mp,0,sizeof mp);
        if(l==1){
            mp[1]=1;sum--;
        }
        for(LL i=1;prime[i]<=sqrt(u);i++){
            if(!prime[i])break;
            LL j=max(2ll,l/prime[i]+(l%prime[i]!=0));
            while(j*prime[i]<=u){
                if(!mp[j*prime[i]-l+1])sum--;
                mp[j*prime[i]-l+1]=true;
                j++;
            }
        }
        if(sum<2){
            printf("There are no adjacent primes.\n");
            continue;
        }
        pre=0;
        for(LL i=l;i<=u;i++){
            if(!mp[i-l+1]){
                if(pre!=0){
                    if(minn>i-pre){minn=i-pre;c1=pre;c2=i;}
                    if(maxx<i-pre){maxx=i-pre;d1=pre;d2=i;}
                }
                pre=i;
            }
        }
        printf("%lld,%lld are closest, %lld,%lld are most distant.\n",c1,c2,d1,d2);
    }
	return 0;
}