結果
問題 | No.1657 Sum is Prime (Easy Version) |
ユーザー | wolgnik |
提出日時 | 2021-08-27 21:32:40 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 265 ms / 2,000 ms |
コード長 | 1,392 bytes |
コンパイル時間 | 212 ms |
コンパイル使用メモリ | 82,404 KB |
実行使用メモリ | 106,864 KB |
最終ジャッジ日時 | 2024-11-21 01:01:49 |
合計ジャッジ時間 | 4,235 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
import sys input = sys.stdin.readline L, R = map(int, input().split()) class sieve: def __init__(self, n): self.n = n self.sv = [1] * (n + 1) self.sv[0] = 0 self.sv[1] = 0 self.minp = [n] * (n + 1) for i in range(2, n + 1): if self.sv[i]: self.minp[i] = i for j in range(i * 2, n + 1, i): self.sv[j] = 0 self.minp[j] = min(self.minp[j], i) def isprime(self, x): if x > self.n: return False return self.sv[x] == 1 def factorize(self, x): if self.sv[x]: return [(x, 1)] res = [] while x > 1: p = self.minp[x] c = 0 while x % p == 0: x //= p c += 1 res.append((p, c)) return res def modlcm(self, a, mod): res = [0] * (self.n + 1) ex = set() for i in range(len(a)): f = self.factorize(a[i]) for j in f: if j > self.n: ex.add(j) continue res[j] = max(f.count(j), res[j]) rres = 1 for i in range(self.n + 1): if res[i] != 0: rres *= pow(i, res[i], mod) rres %= mod for i in ex: rres *= i rres %= mod return rres sv = sieve(R * 2) res = 0 for p in range(R * 2 + 1): if sv.isprime(p): if p in range(L, R + 1): res += 1 x = (p - 1) // 2 y = -((1 - p) // 2) if x in range(L, R + 1) and y in range(L, R + 1): res += 1 print(res)