AcWing 1292. 哥德巴赫猜想

$AcWing$ $1292$. 哥德巴赫猜想

一、题目描述

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

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

例如：

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

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

输入格式
输入包含多组数据。

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

读入以 $0$ 结束。

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

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

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

数据范围
$6≤n<10^6$

输入样例

8
20
42
0

输出样例：

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

二、一些常识

• 欧拉筛时间复杂度

$O(n)$

• 自然对数$e$

自然对数$e$是一个数学常数，约等于$2.71828$。它是一个无理数，也是一个重要的数学常数。$e$是自然对数的底数，它在许多数学和科学领域中都有广泛的应用。

自然对数$e$可以通过以下级数展开来定义：

$\large e = 1 + \frac{1}{1}! + \frac{1}{2}! + \frac{1}{3}! + \frac{1}{4}! + ...$

$n!$表示$n$的阶乘。

• 质数个数定理

$1\sim n$之间质数的个数 约 是$\large \frac{n}{lgn}=\frac{n}{log_en}$

举个例子：$n=100$,按如下的$C++$代码+质数个数定理进行计算,结果是$21$,实际上是$25$个，两者在数量级上是一致的，但不是非常准确，只能用于估算。

$n \div lg(n) =100 \div 4.60517 = 21$

		
#include <iostream>
#include <cmath>

bool is_prime(int num) {
    if (num < 2) {
        return false;
    }

    for (int i = 2; i <= sqrt(num); i++) {
        if (num % i == 0) {
            return false;
        }
    }

    return true;
}

int main() {
    int n = 100; // n的值为100

    int prime_count = 0; // 质数个数的计数器

    for (int i = 2; i <= n; i++) {
        if (is_prime(i)) {
            prime_count++;
        }
    }

    std::cout << "质数个数: " << prime_count << std::endl;

    double estimated_count = (n / log(n)); // 使用质数个数定理估计质数的个数

    std::cout << "质数个数定理估计值: " << estimated_count << std::endl;

    return 0;
}

    
	

• 算数基本定理

任何一个大于$1$的自然数 $N$，如果$N$不为质数，那么$N$可以唯一分解成有限个质数的乘积$N = P_1^{a_1}P_2^{a_2}P_3^{a_3}…P_n^{a_n}$，这里$P_1<P_2<P_3…<P_n$均为质数，其中指数$a_i$是正整数,这样的分解称为 $N$ 的标准分解式。

三、实现思路

• ① 欧拉筛求出区间内所有质数
• ② 枚举每个质数，（注意 要放过数字$2$,因为$2$是偶数，本题要求是奇数）
• ③ 判断$n-a$是不是也是质数，此时$st[i]$数组发挥巨大作用
• ④ 找到的第一个就是结果

四、实现代码

#include <bits/stdc++.h>
using namespace std;

// 欧拉筛
const int N = 1e6 + 10;
int primes[N], cnt; // primes[]存储所有素数
bool st[N];         // st[x]存储x是否被筛掉
void get_primes(int n) {
    for (int i = 2; i <= n; i++) {
        if (!st[i]) primes[cnt++] = i;
        for (int j = 0; primes[j] * i <= n; j++) {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0) break;
        }
    }
}

int main() {
    // 加快读入
    ios::sync_with_stdio(false), cin.tie(0);
    // 欧拉筛求出区间内所有质数
    get_primes(N - 1); // 不能到N，因为数组下标从0开始！要不会越界~

    int n;

    while (cin >> n, n) {
        for (int i = 1;; i++) { // 枚举每个奇数质数,不用上界，因为认为哥德巴赫猜想是正确的，肯定有解
            int a = primes[i];  // primes[0]=2,因为2特殊，是质数，并且是奇数，本题要求是奇数，放过数字2
            int b = n - a;      // 差值
            if (!st[b]) {       // 差值也是质数,这个st数组用的漂亮啊，我还以为要在primes中二分呢
                          // 第一个找的就是最小奇数质数因子，所以相对就是最大的奇数质数因子，两者的差就是最大
                printf("%d = %d + %d\n", n, a, b);
                break;
            }
        }
    }
    return 0;
}