P1304 哥德巴赫猜想 ---素数筛

最新推荐文章于 2025-05-21 12:04:25 发布

chlhoa

最新推荐文章于 2025-05-21 12:04:25 发布

阅读量256

点赞数

文章标签：算法 c++ 数据结构

本文链接：https://blog.csdn.net/weixin_73739312/article/details/129464487

版权

该程序使用C++实现了一个线性筛法求解素数，并应用此方法验证哥德巴赫猜想，即检查所有大于等于4的偶数是否可以表示为两个质数之和。程序首先找到给定范围内的所有素数，然后组合这些素数来寻找符合猜想的配对。

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

**P1304 哥德巴赫猜想 ---素数筛
**来源：洛谷
#include<bits/stdc++.h>
using namespace std;
const int N=1e4+10;
int n;
bool f[N],f1[N];
struct node
{
    int x,y,z;
}a[N];
bool cmp(node x1,node x2)
{
    return x1.z<x2.z;
}
int len=0;
int primes[N],cnt=0;
int main()
{
    cin>>n;
    f[0]=1;
    f[1]=1;
    // 线性筛，求出所有素数
    for(int i=2;i<=n;i++)
    {
        if(!f[i]) primes[++cnt]=i;
        for(int j=1;i*primes[j]<=n;j++)
        {
            f[i*primes[j]]=1;
            if(i%primes[j]==0) break;
        }
    }
    int now;
    for(int i=1;2*primes[i]<=n;i++)
    {
        for(int j=i;j<=cnt;j++)
        {
            now=primes[i]+primes[j];
            if(now>=4&&now<=n&&now%2==0&&f1[now]==0)
            {
                f1[now]=1;
                a[++len].z=now;
                a[len].x=primes[i];
                a[len].y=primes[j];
            }
        }
    }
    sort(a+1,a+1+len,cmp);
    for(int i=1;i<=len;i++)
    {
        cout<< a[i].z<< "="<< a[i].x<< "+"<<a[i].y<<endl;
    }
    return 0;
}