【质数距离-线性筛+埃氏筛】

Kent_J_Truman

于 2024-08-17 20:41:36 发布

阅读量136

点赞数 5

文章标签：算法

本文链接：https://blog.csdn.net/m0_73669127/article/details/141286436

版权

题目

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;

const int N = 1e6+10;
int primes[N], cnt;
bool st[N];
void get_primes(int n)
{
    memset(st, 0, sizeof st);
    cnt = 0;
    for(int i = 2; i <= n; i++)
    {
        if(!st[i]) primes[++cnt] = i;
        for(int j = 1; primes[j] * i <= n; j++)
        {
            st[primes[j] * i] = true;
            if(i % primes[j] == 0) break;
        }
    }
}
int main()
{
    int l, r;
    while(~scanf("%d%d", &l, &r))
    {
        get_primes(50000);
        memset(st, 0, sizeof st);
        for(int i = 1; i <= cnt; i++)
        {
            LL p = primes[i];
            for(LL j = max(2 * p, (l + p - 1) / p * p); j <= r; j += p)
            {
                st[j-l] = true;
            }
        }
        
        cnt = 0;
        for(int i = 0; i <= r - l; i++)
        {
            if(!st[i] && i+l > 1) primes[++cnt] = i + l;
        }
        
        if(cnt < 2)
        {
            cout << "There are no adjacent primes." << '\n';
        }
        else
        {
            int minp = 1, maxp = 1;
            for(int i = 1; i + 1 <= cnt; i++)
            {
                int d = primes[i+1] - primes[i];
                if(d < primes[minp+1] - primes[minp]) minp = i;
                if(d > primes[maxp+1] - primes[maxp]) maxp = i;
            }
            
            printf("%d,%d are closest, %d,%d are most distant.\n"
            , primes[minp], primes[minp+1]
            , primes[maxp], primes[maxp+1]);
        }
    }
}