哥德巴赫猜想之一是指一个偶数(2除外)可以拆分为两个素数之和。请验证这个猜想。
因为同一个偶数可能可以拆分为不同的素数对之和,这里要求结果素数对彼此最接近。
输入格式:
首先输入一个正整数T,表示测试数据的组数,然后是T组测试数据。每组测试输入1个偶数n(6≤n≤10000)。
输出格式:
对于每组测试,输出两个彼此最接近的素数a、b(a≤b),两个素数之间留1个空格。
输入样例:
2
30
40
输出样例:
13 17
17 23
输入格式:
首先输入一个正整数T,表示测试数据的组数,然后是T组测试数据。每组测试输入1个偶数n(6≤n≤10000)。
输出格式:
对于每组测试,输出两个彼此最接近的素数a、b(a≤b),两个素数之间留1个空格。
输入样例:
2
30
40
输出样例:
13 17
17 23
import time
def getPrime(START=2, END = 10000):
data = [1] * END
data[:2] = [0, 0]
for idx in range(int(END ** 0.5) + 1):
if idx > 1:
if data[idx] == 1:
for idy, _ in enumerate(data):
if idy > idx and idy % idx == 0:
data[idy] = 0
return [idx for idx, v in enumerate(data) if END + 1 > idx > START - 1 and v == 1]
n = int(input(">>>"))
l = []
for _ in range(n):
l.append(int(input(">>>")))
t = time.time()
primes = getPrime()
for i in l:
p = [j for j in primes if j < i]
p_t = {tuple(sorted([j, i - j])) for j in p if (i - j) in p}
res = sorted([[abs(k - m),k,m] for k, m in p_t])[0][1:]
print(res[0], res[1])
print(time.time() - t)