|
扩大抽查的范围,使得结果向真值靠拢。
2022-09-12 19:31:19
1158612828906 2201524576849
用时 3611.7531459331512 秒
-------------------------------------------------
print(dt.datetime.now().strftime("%F %T"))
start = tm.time()
m = 1992249998172004
n, n1, n2 = m + 2, 0, 0
m0 = m // 6469693230 >> 1
m1 = m % 6469693230 >> 2 << 1
for i in range(m0):
j = i * 6469693230 + m1
for k in range(1, 2310, 2):
l = j - k
if isPrime(l):
if isPrime(m - l): n1 += 1
if isPrime(n - l): n2 += 1
print(f'{n1 * 2800733} {n2 * 2800733}')
print("用时 {0} 秒".format(tm.time() - start)) |
|