poj 2689 巧妙地运用素数筛选

poj 2689 巧妙地运用素数筛选

大家好,又见面了,我是全栈君,今天给大家准备了Idea注册码。

称号:

   给出一个区间[L,R]求在该区间内的素数最短,最长距离。 (R < 2 * 10^9 , R – L <= 10 ^ 6)

   由数论知识可得一个数的因子可在开根号内得到。

所以,我们能够打出5*10^4内得素数。然后,在用一次筛法把在[L。R]内得合数找到,则剩下的就是素数了。这里要用到离散化。把一个数 x – L 保存在数组里。由于,直接保存肯定不行。可是我们发现区间特点较小。所以。能够想到离散化。

 

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;

typedef long long LL;
const int MAXN = 50000;
int primes[MAXN];
bool vst[MAXN];
int notPrimes[1000010];
int pos[1000010];
int top,pcnt;

void init(){
    top = 0;
    memset(vst,0,sizeof(vst));
    vst[0] = vst[1] = 1;
    for(int i = 2;i < MAXN;++i)if(!vst[i]){
        primes[top++] = i;
        for(int j = i + i;j < MAXN;j += i) vst[j] = 1;
    }
    //printf("top: %d\n",top);
}

void solve(int L,int R){
    memset(notPrimes,0,sizeof(notPrimes));

    if(L == 1) L = 2;              /// 防止筛掉全部该区间的素数本身!!!!!
    for(int i = 0;i < top&&(LL)primes[i]*primes[i] <= R;++i){  //筛选因子
        int s = L / primes[i] + (L % primes[i] > 0); //当前素数的最小倍数达到L
       
        s = (s == 1 ? 2 : s);    /// 防止筛掉全部该区间的素数本身!!!!!
       
        for(int j = s;(LL)j*primes[i] <= R;++j){
            if((LL)j*primes[i] >= L)   //合数
               notPrimes[j*primes[i] - L] = 1;   // 相当与离散化
        }
    }

    pcnt = 0;
    for(int i = 0;i <= R - L;++i){
        if(!notPrimes[i]){
            pos[pcnt++] = i + L;
            //printf("i -- > %d\n",i + L);
        }
    }
    
    
    if(pcnt < 2){
        puts("There are no adjacent primes.");
    } else {
        int minl,minr,maxl,maxr,minv = 999999,maxv = -1;
        for(int i = 1;i < pcnt;++i){
            if(pos[i] - pos[i-1] > maxv){
                maxv = pos[i] - pos[i-1];
                maxl = pos[i-1];
                maxr = pos[i];
            }
            if(pos[i] - pos[i-1] < minv){
                minv = pos[i] - pos[i-1];
                minl = pos[i-1];
                minr = pos[i];
            }
        }
        printf("%d,%d are closest, %d,%d are most distant.\n",minl,minr,maxl,maxr);
    }
}
int main()
{
    init();
    int L,R;
    while(~scanf("%d%d",&L,&R)){
        solve(L,R);
    }
    return 0;
}




 

版权声明:本文博客原创文章,博客,未经同意,不得转载。

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/117277.html原文链接:https://javaforall.net

(0)
上一篇 2022年1月8日 上午11:00
下一篇 2022年1月8日 下午12:00


相关推荐

  • mysql密码字段类型_MySQL 字段类型[通俗易懂]

    mysql密码字段类型_MySQL 字段类型[通俗易懂]数值MySQL的数值数据类型可以大致划分为两个类别,一个是整数,另一个是浮点数或小数。许多不同的子类型对这些类别中的每一个都是可用的,每个子类型支持不同大小的数据,并且MySQL允许我们指定数值字段中的值是否有正负之分(UNSIGNED)或者用零填补(ZEROFILL)。INT在MySQL中支持的5个主要整数类型是TINYINT,SMALLINT,MEDIUMINT,INT和B…

    2025年7月30日
    7
  • 【TCP/IP】IP地址分类和特殊IP地址

    【TCP/IP】IP地址分类和特殊IP地址IP地址是因特网技术中的一个非常重要的概念,IP地址在IP层实现了底层网络地址的统一,使因特网的网络层地址具有全局唯一性和一致性。IP地址含有位置信息,反映了主机的网络连接,使因特网进行寻址和路由选择的依据。 IP地址概述地址是标识对象所处位置的标识符。传输中的信息带有源地址和目的地址,分别标识通信的源结点和目的结点,即信源和信宿。目的地址是传输设备为信息进行寻址的依据。不同的物理…

    2022年4月29日
    72
  • 2026年AI大模型应用开发保姆级教程:从入门到精通,这一篇开始

    2026年AI大模型应用开发保姆级教程:从入门到精通,这一篇开始

    2026年3月15日
    2
  • 高斯滤波器详解

    高斯滤波器详解高斯滤波器详解本文主要介绍了高斯滤波器的原理及其实现过程高斯滤波器是一种线性滤波器 能够有效的抑制噪声 平滑图像 其作用原理和均值滤波器类似 都是取滤波器窗口内的像素的均值作为输出 其窗口模板的系数和均值滤波器不同 均值滤波器的模板系数都是相同的为 1 而高斯滤波器的模板系数 则随着距离模板中心的增大而系数减小 所以 高斯滤波器相比于均值滤波器对图像个模糊程度较小 什

    2026年3月18日
    1
  • 移位运算用法总结

    移位运算用法总结位运算总结原文一、位运算应用口诀清零取位要用与,某位置一可用或若要取反和交换,轻轻松松用异或二、移位运算它们都是双目运算符,两个运算分量都是整形,结果也是整形。‘&lt;&lt;’左移:右边空出的位置补0,其值相当于乘以2。‘&gt;&gt;’右移:左边空出的位,如果是正数则补0,若为负数则补0或1,取决于所用的计算机系统OSX中补1。其值相当于除以2。…

    2022年7月13日
    19
  • DeepSeek使用指南:从入门到精通的全方位解析

    DeepSeek使用指南:从入门到精通的全方位解析

    2026年3月12日
    2

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号