結果

問題 No.843 Triple Primes
ユーザー masatosismasatosis
提出日時 2023-05-16 13:25:27
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 331 ms / 2,000 ms
コード長 894 bytes
コンパイル時間 1,702 ms
コンパイル使用メモリ 87,092 KB
実行使用メモリ 100,500 KB
最終ジャッジ日時 2023-08-21 05:18:34
合計ジャッジ時間 11,415 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 74 ms
71,864 KB
testcase_01 AC 331 ms
100,500 KB
testcase_02 AC 83 ms
76,788 KB
testcase_03 AC 79 ms
76,728 KB
testcase_04 AC 81 ms
76,800 KB
testcase_05 AC 79 ms
76,896 KB
testcase_06 AC 84 ms
76,824 KB
testcase_07 AC 291 ms
95,328 KB
testcase_08 AC 276 ms
97,724 KB
testcase_09 AC 321 ms
100,436 KB
testcase_10 AC 265 ms
97,212 KB
testcase_11 AC 293 ms
99,668 KB
testcase_12 AC 314 ms
100,340 KB
testcase_13 AC 301 ms
100,052 KB
testcase_14 AC 305 ms
99,908 KB
testcase_15 AC 295 ms
96,960 KB
testcase_16 AC 259 ms
97,332 KB
testcase_17 AC 68 ms
71,584 KB
testcase_18 AC 69 ms
71,636 KB
testcase_19 AC 70 ms
71,480 KB
testcase_20 AC 121 ms
81,872 KB
testcase_21 AC 100 ms
78,776 KB
testcase_22 AC 204 ms
90,748 KB
testcase_23 AC 212 ms
90,680 KB
testcase_24 AC 123 ms
82,216 KB
testcase_25 AC 108 ms
80,320 KB
testcase_26 AC 321 ms
100,312 KB
testcase_27 AC 83 ms
77,232 KB
testcase_28 AC 304 ms
100,072 KB
testcase_29 AC 126 ms
82,868 KB
testcase_30 AC 315 ms
100,248 KB
testcase_31 AC 94 ms
77,672 KB
testcase_32 AC 82 ms
77,084 KB
testcase_33 AC 150 ms
83,824 KB
testcase_34 AC 226 ms
88,028 KB
testcase_35 AC 321 ms
100,384 KB
testcase_36 AC 102 ms
79,712 KB
testcase_37 AC 279 ms
95,320 KB
testcase_38 AC 208 ms
92,524 KB
testcase_39 AC 307 ms
100,196 KB
testcase_40 AC 69 ms
71,720 KB
testcase_41 AC 68 ms
71,900 KB
testcase_42 AC 267 ms
97,428 KB
testcase_43 AC 288 ms
99,916 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