結果

問題 No.843 Triple Primes
ユーザー convexineq
提出日時 2019-06-28 21:31:57
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 189 ms / 2,000 ms
コード長 625 bytes
コンパイル時間 84 ms
コンパイル使用メモリ 12,416 KB
実行使用メモリ 20,572 KB
最終ジャッジ日時 2024-09-19 13:49:43
合計ジャッジ時間 5,726 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff #

# coding: utf-8
# Your code here!

n=int(input())
if n==1: print(0); exit()



def Eratosthenes(N):
    N+=1
    is_prime_list = [True]*N
    is_prime_list[0], is_prime_list[1] = False, False
    for j in range(4,N,2): is_prime_list[j] = False
    m = min(int(N**0.5)+10,N)
    for i in range(3,m,2):
        if is_prime_list[i]:
            for j in range(i*i,N,i):
                is_prime_list[j] = False
    return [number for number, TF in enumerate(is_prime_list) if TF]

heiho = {i*i for i in Eratosthenes(n)}

ans = 0
for i in Eratosthenes(n+1):
    if i+2 in heiho:
        ans += 1
#        print(i)

print(ans*2-1)
0