信息学奥赛一本通 1622:Goldbach’s Conjecture | 洛谷 UVA543 Goldbach‘s Conjecture

【题目链接】

ybt 1622:Goldbach's Conjecture
洛谷 UVA543 Goldbach's Conjecture

【题目考点】

1. 筛法求质数表

【解题思路】

首先使用埃筛或线性筛求出质数表。

包括isPrime数组,isPrime[i]表示数值i是否是质数。以及prime数组,prime[i]保存第i个质数,pn是保存在prime数组中的质数的个数。

判断整数n是否可以写成两个奇素数的加和,枚举第一个较小的奇素数。prime[1]是2,是偶数,略过。i从2循环到pn,第一个奇素数为prime[i],而该数是相加的两个素数中的较小的数,因此该数需要满足不超过n的一半,即需要满足prime[i] <= n/2。第一个奇素数是prime[i],要想使两个数加和为n,则第二个数为n-prime[i],判断第二个数是否是素数,可以使用前面求出的isPrime数组,isPrime[n-prime[i]]表示n-prime[i]是否为素数。如果第二个数也是素数,则输出n等于两个数相加的公式,并跳出循环。

虽然哥的巴赫猜想还没有被证明,但在题目给定的范围找到一个偶数不满足哥的巴赫猜想是不可能的,如果你找到了哥的巴赫猜想的反例,都可以得菲尔兹奖了。所以不需要考虑无解的情况。

【题解代码】

解法1:埃筛求质数表

cpp 复制代码
#include <bits/stdc++.h>
using namespace std;
#define N 1000005
bool isPrime[N];
int prime[N], pn;
void initPrime(int n)
{
	memset(isPrime, 1, sizeof(isPrime));
	for(int i = 2; i*i <= n; ++i) if(isPrime[i])
		for(int j = i*i; j <= n; j += i)
			isPrime[j] = false;
	for(int i = 1; i <= n; ++i) if(isPrime[i])
		prime[++pn] = i;
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(nullptr); 
	initPrime(1e6);
	int n;
	while(cin >> n && n != 0)
	{
		for(int i = 2; i <= pn && prime[i] <= n/2; ++i) if(isPrime[n-prime[i]])
		{	
			cout << n << " = " << prime[i] << " + " << n-prime[i] << '\n';
			break;
		}
	}
	return 0;
}

解法2:线性筛求质数表

cpp 复制代码
#include <bits/stdc++.h>
using namespace std;
#define N 1000005
bool isPrime[N];
int prime[N], pn;
void initPrime(int n)
{
	memset(isPrime, 1, sizeof(isPrime));
	for(int i = 2; i <= n; ++i)
	{
		if(isPrime[i])
			prime[++pn] = i;
		for(int j = 1; j <= pn && i*prime[j] <= n; ++j)
		{
			isPrime[i*prime[j]] = false;
			if(i%prime[j] == 0)
				break;
		}
	}
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(nullptr); 
	initPrime(1e6);
	int n;
	while(cin >> n && n != 0)
	{
		for(int i = 2; i <= pn && prime[i] <= n/2; ++i) if(isPrime[n-prime[i]])
		{	
			cout << n << " = " << prime[i] << " + " << n-prime[i] << '\n';
			break;
		}
	}
	return 0;
}
相关推荐
CoovallyAIHub2 分钟前
超越传统3D生成:OccScene实现感知与生成的跨任务共赢
深度学习·算法·计算机视觉
Mr.H012711 分钟前
克鲁斯卡尔(Kruskal)算法
数据结构·算法·图论
Tisfy13 分钟前
LeetCode 3346.执行操作后元素的最高频率 I:滑动窗口(正好适合本题数据,II再另某他法)
算法·leetcode·题解·滑动窗口·哈希表
黑菜钟23 分钟前
代码随想录第53天 | 图论二三题
c++·图论
CoovallyAIHub26 分钟前
华为世界模型来了!30分钟生成272㎡室内场景,虚拟人导航不迷路
深度学习·算法·计算机视觉
西哥写代码27 分钟前
基于dcmtk的dicom工具 第十二章 响应鼠标消息实现图像的调窗、缩放、移动
c++·mfc·dicom·dcmtk·vs2017
熬了夜的程序员29 分钟前
【LeetCode】94. 二叉树的中序遍历
数据结构·算法·leetcode·职场和发展·深度优先
放羊郎30 分钟前
SLAM各类算法特点对比
人工智能·算法·slam·视觉slam·建图·激光slam
熬了夜的程序员34 分钟前
【LeetCode】92. 反转链表 II
数据结构·算法·leetcode·链表·职场和发展·排序算法
Aurorar0rua41 分钟前
C Primer Plus Notes 10
c语言·开发语言·算法