結果
| 問題 | 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 | 
ソースコード
# 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)
            
            
            
        