标签：const int 196 long maxn primes acwing define

#include <bits/stdc++.h>
#define IOS ios_base::sync_with_stdio(0),cin.tie(0),cout.tie(0)
#define _zero(a) memset(a, 0, sizeof(a))
#define endl '\n'
#define int long long
#define mp make_pair
#define PII pair<int, int>
#define x first
#define y second
#define all(a) (a).begin(), (a).end()
#define rall(a) (a).rbegin(), (a).rend()
#define PII pair<int, int>
typedef long long ll;
typedef double dd;
typedef long double ld;

using namespace std;

const int inf = 1e10;
const int M = 998244353;
const ld pi = atan2(0, -1);//arctan(y/x);
const ld eps = 1e-8;

const int maxn = 5e5 + 100, N = 1e6 + 100;
int lp[maxn],primes[maxn],pcnt;
bool st[N];

void sieve()//线性筛
{
    pcnt=0;
    fill_n(lp,maxn,0);
    for(int i=2;i<maxn;i++)
    {
        int& l=lp[i];
        if(l==0)l=i,primes[pcnt++]=i;
        for(int j=0;j<pcnt&&primes[j]<=l;++j)
        {
            int p=primes[j];
            if(i*p>=maxn)break;
            lp[i*p]=p;
        }
    }
}

signed main()
{
    IOS;
    int l, r;
    sieve();
    while (cin >> l >> r) {
        memset(st, false, sizeof st);
        if (l == 1) st[0] = true;
        
        for (int i = 0; i < pcnt; i ++) {
            int p = primes[i];
            //aishai
            int j = ((l + p - 1) / p);
            if (j == 1) ++j;
            j *= p;
            for (; j <= r; j += p) 
                st[j - l] = true;
        }
        
        vector<int> a;
        for (int i = l; i <= r; i ++) if (!st[i - l]) a.push_back(i);
        
        int idx1 = -1, idx2 = -1, mn = INT_MAX, mx = 0;
        if (a.size() < 2) {
            printf("There are no adjacent primes.\n");
            continue;
        }
        for (int i = 0; i < a.size() - 1; i ++) {
            int k = a[i + 1] - a[i];
            if (k > mx) {
                idx2 = i;
                mx = k;
            }
            if (k < mn) {
                mn = k;
                idx1 = i;
            }
        }
        printf("%lld,%lld are closest, %lld,%lld are most distant.\n", a[idx1], a[idx1 + 1], a[idx2], a[idx2 + 1]);
        
        //for (int i = l; i <= r; i ++) cout << st[i - l] << endl;
        //for (int i : a) cout << i << endl;
        //for (int i = 0; i < 100; i ++) cout << primes[i] << endl;
    }
    return 0;
}
 

标签：const,int,196,long,maxn,primes,acwing,define
来源： https://blog.csdn.net/ssd778/article/details/122564413

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。