POJ2689——区间筛法——Prime Distance

http://poj.org/problem?id=2689

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
/*
区间筛法适用于L, R范围比较大,但是区间长度较小时
套个模板然后取出最小最大就可以
*/
/************************************************
* Author        :Powatr
* Created Time  :2015-8-17 9:38:50
* File Name     :POJ2689.cpp
 ************************************************/
 
#include <cstdio>
#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstring>
#include <cmath>
#include <string>
#include <vector>
#include <queue>
#include <deque>
#include <stack>
#include <list>
#include <map>
#include <set>
#include <bitset>
#include <cstdlib>
#include <ctime>
using namespace std;
 
#define lson l, mid, rt << 1
#define rson mid + 1, r, rt << 1 | 1
typedef long long ll;
const int MAXN = 1e6 + 1e4;
const int INF = 0x3f3f3f3f;
const int MOD = 1e9 + 7;
const int N = 50010;
 
bool prime[MAXN];
bool sieve[N];
int len;
void segent_sieve(ll L, ll R)
{
    len = R - L + 1;
    for(int i = 0; i < len; i++) prime[i] = true;
    if(1 - L >= 0) prime[1-L] = false;//特判
    for(ll i = 2; i*i < R; i++){
        if(seive[i]){
            for(ll j = max(1ll*2, (L-1+i)/i)*i; j <= R; j+=i)//(L-1+i/i)得到最接近L的i的倍数
                prime[j-L] = false;
            //偏移量,因为单纯j会导致数组越界,把范围弄到0到len上
        }
    }
}
int main(){
    //艾氏筛法 O(nlognlogn)
    for(int i = 2; i < N; i++) sieve[i] = true;
    for(int i = 2; i*i < N; i++){
        if(sieve[i]){
            for(int j = i*2; j < N; j+=i){
                sieve[j] = false;
            }
        }
    }
    int L, R;
    while(~scanf("%d%d", &L, &R)){
        segent_seive(L, R);
    int mmax, mmin;
    int lmax, rmax, lmin, rmin;
    mmax = -1, mmin = 1 << 30;
    int t = -1;
    for(int i = 0; i < len; i++){
        if(prime[i]){
            if(t >= 0){
                if(mmax < i - t) mmax = i - t, lmax = L + t, rmax = L + i;
                if(mmin > i - t) mmin = i - t, lmin = L + t, rmin = L + i;
                t = i;
            }
            else t = i;
        }
    }
    if(mmax != -1 || mmin != 1 << 30) printf("%d,%d are closest, %d,%d are most distant.\n", lmin, rmin, lmax, rmax);
    else printf("There are no adjacent primes.\n");
    }
    return 0;
}

  

posted @   Painting、时光  阅读(290)  评论(0编辑  收藏  举报
编辑推荐:
· 智能桌面机器人:用.NET IoT库控制舵机并多方法播放表情
· Linux glibc自带哈希表的用例及性能测试
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
阅读排行:
· 手把手教你在本地部署DeepSeek R1,搭建web-ui ,建议收藏!
· 新年开篇:在本地部署DeepSeek大模型实现联网增强的AI应用
· Janus Pro:DeepSeek 开源革新,多模态 AI 的未来
· 互联网不景气了那就玩玩嵌入式吧,用纯.NET开发并制作一个智能桌面机器人(三):用.NET IoT库
· 【非技术】说说2024年我都干了些啥
点击右上角即可分享
微信分享提示