[线筛] [洛谷] P2640

最新推荐文章于 2024-07-09 01:32:36 发布

原创最新推荐文章于 2024-07-09 01:32:36 发布 · 217 阅读

0 ·

CC 4.0 BY-SA版权

数论专栏收录该内容

25 篇文章

订阅专栏

本文介绍了一种高效的素数筛选算法，并利用该算法查找一定范围内的孪生素数对。通过初始化一个布尔型数组，标记所有可能的合数，从而得到素数列表。然后，遍历这个列表，找出所有相差为特定值M的孪生素数对。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

//#pragma GCC optimize(2)
#include <cstdio>
#include <iostream>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <string>
#include <cstring>
#include <algorithm>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <ctime>
#include <vector>
#include <fstream>
#include <list>
#include <iomanip>
#include <numeric>
using namespace std;
typedef long long ll;

const int MAXN = 1e6 + 10;

bool isprime[MAXN + 100];

int prime[MAXN + 100], tf = 0;

void init()
{
    for(int i = 2; i < MAXN; i++)
    {
        if(!isprime[i])
            prime[tf++] = i;

        for(int j = 0; j < MAXN && i * prime[j] < MAXN; j++)
        {
            isprime[i * prime[j]] = 1;
            if(i * prime[j] == 0)
                break;
        }
    }

    /*for(int i = 0; i < 100; i++)
        cout<<prime[i]<<endl;*/
}

int main()
{
    //ios::sync_with_stdio(false);

    //cin.tie(0);     cout.tie(0);

    init();

    int N, M, ans = 0;

    cin>>N>>M;

    bool flag = false;

    for(int i = 2; i <= N - M; i++)
    {
        if(!isprime[i] && !isprime[i + M])
        {
            flag = true;

            cout<<i<<' '<<i + M<<endl;
        }
    }

    flag ? cout<<endl : cout<<"empty"<<endl;

    return 0;
}