CQUPT Python 哥德巴赫猜想

【问题描述】

数学领域著名的“哥德巴赫猜想”的大致意思是：任何一个大于2的偶数总能表示为两个素数之和。比如：24=5+19，其中5和19都是素数。本实验的任务是设计一个程序，验证20亿以内的偶数都可以分解成两个素数之和。输入一个大于2的正整数，当输入为偶数时，在一行中按照格式“N = p + q”输出N的素数分解，其中p 、 q均为素数且p ≤ q。因为这样的分解可能不唯一（例如24还可以分解为7+17），要求必须输出所有解中p最小的解。当输入为奇数时，输出'Data error!' 。

【输入形式】

输入一个不小于2的正整数

【输出形式】

当输入为偶数时，按照格式“N = p + q”输出N的素数分解；当输入为奇数时，输出'Data error!' 。‪‬‪‬‪‬‪‬‪‬‮‬‪‬‭‬‪‬‪‬‪‬‪‬‪‬‮‬‭‬‪‬‪‬‪‬‪‬‪‬‪‬‮‬‫‬‫‬‪‬‪‬‪‬‪‬‪‬‮‬‭‬‫‬‪‬‪‬‪‬‪‬‪‬‮‬‭‬‫‬‪‬‪‬‪‬‪‬‪‬‮‬‪‬‮‬

【样例输入】

36
【样例输出】

N = 5 + 31

【完整代码】

# 判断n是否为素数
def isprime(n):
    for i in range(2, int(n ** 0.5) + 1):
        if n % i == 0:
            return False
    return True


n = int(input())
if n % 2 == 1:
    print("Data error!")
else:
    for i in range(2, n):
        if isprime(i):  # p是素数
            if isprime(n - i):  # q也是素数
                print(f"N = {i} + {n - i}")
                break

【代码讲解】

ispreime函数判断n是否为素数，遍历n存不存在除了1和它本身的因素。

本来以为20亿会造成超时，结果只找一个p最小的解，没花多久挺快的。