結果

問題 No.843 Triple Primes
ユーザー masatosismasatosis
提出日時 2023-05-16 13:25:27
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 309 ms / 2,000 ms
コード長 894 bytes
コンパイル時間 183 ms
コンパイル使用メモリ 82,156 KB
実行使用メモリ 95,040 KB
最終ジャッジ日時 2024-12-14 14:10:24
合計ジャッジ時間 8,690 ms
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
53,832 KB
testcase_01 AC 309 ms
95,040 KB
testcase_02 AC 49 ms
62,632 KB
testcase_03 AC 49 ms
62,400 KB
testcase_04 AC 50 ms
62,436 KB
testcase_05 AC 50 ms
63,320 KB
testcase_06 AC 50 ms
63,588 KB
testcase_07 AC 265 ms
87,832 KB
testcase_08 AC 263 ms
90,428 KB
testcase_09 AC 290 ms
94,372 KB
testcase_10 AC 234 ms
90,844 KB
testcase_11 AC 270 ms
94,588 KB
testcase_12 AC 284 ms
93,420 KB
testcase_13 AC 280 ms
93,616 KB
testcase_14 AC 278 ms
93,208 KB
testcase_15 AC 272 ms
90,332 KB
testcase_16 AC 233 ms
89,664 KB
testcase_17 AC 39 ms
52,816 KB
testcase_18 AC 38 ms
53,480 KB
testcase_19 AC 38 ms
53,348 KB
testcase_20 AC 90 ms
71,044 KB
testcase_21 AC 70 ms
67,556 KB
testcase_22 AC 183 ms
82,008 KB
testcase_23 AC 193 ms
82,984 KB
testcase_24 AC 95 ms
71,408 KB
testcase_25 AC 79 ms
69,616 KB
testcase_26 AC 295 ms
94,832 KB
testcase_27 AC 54 ms
63,600 KB
testcase_28 AC 285 ms
93,440 KB
testcase_29 AC 102 ms
73,352 KB
testcase_30 AC 289 ms
93,592 KB
testcase_31 AC 63 ms
65,084 KB
testcase_32 AC 54 ms
64,544 KB
testcase_33 AC 123 ms
74,748 KB
testcase_34 AC 210 ms
79,948 KB
testcase_35 AC 294 ms
94,480 KB
testcase_36 AC 74 ms
68,884 KB
testcase_37 AC 251 ms
88,068 KB
testcase_38 AC 179 ms
84,376 KB
testcase_39 AC 281 ms
93,488 KB
testcase_40 AC 39 ms
52,868 KB
testcase_41 AC 38 ms
53,104 KB
testcase_42 AC 238 ms
89,964 KB
testcase_43 AC 266 ms
93,252 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

N = int(input())


def get_primes(n):
    """nまでの素数リスト"""
    candidates = [0, 1] * (n // 2 + 1)
    if n % 2 == 0:
        candidates.pop()
    candidates[1] = 0
    candidates[2] = 1
    for p in range(3, int(n**0.5) + 1, 2):
        if candidates[p]:
            for m in range(p * p, n + 1, 2 * p):
                candidates[m] = 0
    return [p for p, b in enumerate(candidates) if b]


primes = get_primes(2 * N)
setprimes = set(primes)
# print(len(primes))


# R2 = dict(x: x * x for x in primes)
R2 = dict()
# R2 = set()
for x in primes:
    if x * x > 2 * N:
        break
    R2[x * x] = x
    # R2.add(x*x)
# print(len(R2))

ans = set()
for p in primes:
    if p > N:
        break
    for r2 in R2:
        if R2[r2] > N:
            break
        q = r2 - p
        if q <= N and q in setprimes:
            ans.add((p, q, R2[p + q]))
print(len(ans))
# print(ans)
0