20/07/19 19:37:07.01 G6h2CBYL.net
>>277
>>294
10億まで調べたが（GPUで10分ぐらい）それ以降のが見つからんな
import math
import sys
from tqdm import tqdm
MAX=int(sys.argv[1])
merk=int(MAX/20)
record=dict()
sosu=[2,3,]
kinko_sosu=[2,3,]
def isPrime(n):
m = math.floor(math.sqrt(n)) + 1
for p in sosu:
if n % p == 0:
return False
if p >= m:
return True
cnt_1mod6=cnt_5mod6=0
cnt_p=2
for i in tqdm(range(5,MAX)):
if isPrime(i):
cnt_p+=1
sosu.append(i)
if i % 6 == 1:
cnt_1mod6 += 1
elif i % 6 == 5:
cnt_5mod6 += 1
if cnt_1mod6 == cnt_5mod6:
kinko_sosu.append(i)
if i % merk == 0:
record[i] = str(round(100*(cnt_5mod6 - cnt_1mod6)/cnt_1mod6,3))+'%'
print(kinko_sosu)
print(record)
100%|██████████| 99999995/99999995 [09:58<00:00, 167059.18it/s]
[2, 3, 7, 13, 19, 37, 43, 79, 163, 223, 229]
{5000000: '0.075%', 10000000: '0.057%', 15000000: '0.063%', 20000000: '0.056%',
25000000: '0.027%', 30000000: '0.041%', 35000000: '0.054%', 40000000: '0.038%',
45000000: '0.039%', 50000000: '0.011%', 55000000: '0.028%', 60000000: '0.02%',
65000000: '0.033%', 70000000: '0.033%', 75000000: '0.029%', 80000000: '0.028%',
85000000: '0.02