結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

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)
0