【浙大版《Python 程序设计》题目集（解）】第4章-4 验证“哥德巴赫猜想”（20分）

数学领域著名的“哥德巴赫猜想”的大致意思是：任何一个大于2的偶数总能表示为两个素数之和。比如：24=5+19，其中5和19都是素数。本实验的任务是设计一个程序，验证20亿以内的偶数都可以分解成两个素数之和。

输入格式：

输入在一行中给出一个(2, 2 000 000 000]范围内的偶数N。

输出格式：

在一行中按照格式“N = p + q”输出N的素数分解，其中p ≤ q均为素数。又因为这样的分解不唯一（例如24还可以分解为7+17），要求必须输出所有解中p最小的解。

输入样例：

24

输出样例：

24 = 5 + 19

解析:

本题要验证哥德巴赫猜想，证明任何一个大于2的偶数总能表示为两个素数之和。

方案 1 思想（不推荐）：

        建立素数库，使用素数库里的素数去验证。不推荐，原因超时！！！

方案 2 思想：

        找到第一个素数，计算差值，验证差值是否也是素数。

方案 3 思想：（新手会有有点难受）

        方案 3 采用了引入‘库’以及‘函数‘，如果是学一个单元做相应单元的练习的学习者，这一步就比较难受了。回到正题：建立 t_f 函数，用于判断参数是否为素数，返回值只有 'True' 和 'False'，通过差值来验证是否都为素数，是素数则输出其结果并结束循环，原理和方案 2 同理。

方案 4 思想（不推荐）：

        使用列表推导式获取一定范围内的素数库，和方案 1 同理 ，不推荐，放在这主要原因是’短‘！！！可以减少代码量

解题:

# 方案1：（不推荐）
# 输入大于2的偶数n
n = int(input())

# 创建素数库
sushu_list = list()
for i in range(2, n):
    for j in range(2, int(i ** 0.5) + 1):
        if i % j == 0:
            break
    else:
        sushu_list.append(i)

# 验证‘哥德巴赫猜想’
flag = True
for i in sushu_list[:-1]:
    for j in sushu_list[1:]:
        if n == i + j:
            print("{0:} = {1:} + {2:}".format(n, i, j))
            flag = False
            break
    if flag is False:
        break

#####################################################################

# 方案2：
# 输入任何一个大于2的偶数
n = int(input())

# 验证‘哥德巴赫猜想’
flag = True  # 控制循环的结束
for i in range(2, n // 2 + 1):  # 单个素数小于n/2
    for j in range(2, int(i ** 0.5) + 1):  # 寻找素数
        if i % j == 0:
            break
    else:  # 找到素数去验证’哥德巴赫猜想‘
        for k in range(2, int((n - i) ** 0.5 + 1)):  # 查看另一个是否也是素数
            if (n - i) % k == 0:
                break
        else:  # 验证成功，打印结果
            print("{0:} = {1:} + {2:}".format(n, i, n - i))
            flag = False
    if flag is False:
        break  # 结束外循环

########################################################################

# 方案3：
# 引用数学库
import math

# 创建用于判别素数函数
def t_f(n):
  if n <= 1:
    return False
  for i in range(2, int(math.sqrt(n)) + 1):
    if n % i == 0:
        return False
  return True

# 输入任何一个大于2的偶数
c = int(input())

for a in range(2, c // 2 + 1):
    b = c - a
    if (t_f(a) is True and t_f(b) is True):
        print('{} = {} + {}'.format(c, a, b))
        break

############################################################################

# 方案4（不推荐）
n = eval(input())
a = [i for i in range(2, n) if 0 not in [i % j for j in range(2, i)]]
for i in a:
    if n - i in a:
        print("{} = {} + {}".format(n, i, n - i))
        break

如有更好的建议请评论！！！