結果
| 問題 |
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)
wolgnik