結果
| 問題 | No.843 Triple Primes | 
| コンテスト | |
| ユーザー |  c-yan | 
| 提出日時 | 2020-12-04 00:25:12 | 
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 150 ms / 2,000 ms | 
| コード長 | 793 bytes | 
| コンパイル時間 | 178 ms | 
| コンパイル使用メモリ | 12,672 KB | 
| 実行使用メモリ | 30,904 KB | 
| 最終ジャッジ日時 | 2024-09-14 09:53:11 | 
| 合計ジャッジ時間 | 5,462 ms | 
| ジャッジサーバーID (参考情報) | judge4 / judge2 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 2 | 
| other | AC * 42 | 
ソースコード
def enumerate_primes(n):
    if n < 2:
        return []
    sieve = list(range(n + 1))
    result = [2]
    for i in range(4, n + 1, 2):
        sieve[i] = 2
    m = int(n ** 0.5)
    if m % 2 == 0:
        m -= 1
    for i in range(3, m + 1, 2):
        if sieve[i] != i:
            continue
        result.append(i)
        for j in range(i * i, n + 1, i * 2):
            if sieve[j] == j:
                sieve[j] = i
    for i in range(m + 2, n + 1, 2):
        if sieve[i] == i:
            result.append(i)
    return result
N = int(input())
if N == 1:
    print(0)
    exit()
elif N == 2:
    print(1)
    exit()
result = 1
ps = enumerate_primes(N)
ops = set(ps[1:])
for r in ps[1:]:
    if r * r - 2 > N:
        break
    if r * r - 2 in ops:
        result += 2
print(result)
            
            
            
        