AcWing哥德巴赫猜想(C++)筛素数模板题

题目

哥德巴赫猜想的内容如下：

任意一个大于 44 的偶数都可以拆成两个奇素数之和。

例如：

8=3+58=3+5
20=3+17=7+1320=3+17=7+13
42=5+37=11+31=13+29=19+2342=5+37=11+31=13+29=19+23

现在，你的任务是验证所有小于一百万的偶数能否满足哥德巴赫猜想。

输入格式

输入包含多组数据。

每组数据占一行，包含一个偶数 n。

读入以 00 结束。

输出格式

对于每组数据，输出形如 n = a + b，其中 a,b是奇素数。

若有多组满足条件的 a,b，输出 b−a最大的一组。

若无解，输出 Goldbach's conjecture is wrong.。

数据范围

输入最多包含 5000650006 组数据。
6≤n<1066≤<106

输入样例：

8
20
42
0

输出样例：

8 = 3 + 5
20 = 3 + 17
42 = 5 + 37

代码

#include<iostream>

using namespace std;

const int N = 1e6 + 10;
int primes[N], cnt;
bool st[N];

//线性筛质数模板 
void init()
{
	for(int i = 2; i <= N; i ++)
	{
		if(!st[i]) primes[cnt ++] = i;
			
		for(int j = 0; j < cnt && (primes[j] * i < N); j ++)
		{
			st[primes[j] * i] = true;
			if(i % primes[j] == 0) break;
		}
	}
}

int main()
{
	//初始化，先预处理出1e6之内的质数 
	init();
	
	int x;
	while(cin >> x && x)
	{
		for(int i = 1; ; i ++) //这里从下标1开始遍历 ，是因为题目中要求奇质数，而第一个素数为2，因此从3开始遍历 
		{
			int b = x - primes[i];
			if(!st[b]) //如果这个数也是素数 
			{
				cout << x << " = " << primes[i] << " + " << b << endl;//哥德巴赫猜想目前没被证伪，因此不用考虑不存在情况 
				break;
			}	
		} 
	}
	
	return 0;
}