結果
| 問題 |
No.371 ぼく悪いプライムじゃないよ
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2016-05-17 00:02:03 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
AC
|
| 実行時間 | 68 ms / 1,000 ms |
| コード長 | 824 bytes |
| コンパイル時間 | 329 ms |
| コンパイル使用メモリ | 12,544 KB |
| 実行使用メモリ | 12,160 KB |
| 最終ジャッジ日時 | 2024-12-30 18:19:05 |
| 合計ジャッジ時間 | 4,655 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 42 |
ソースコード
#!/usr/bin/python3
def sieve(n):
is_prime = [0, 0] + [1] * (n-1)
sq = int(n ** .5)
for i in range(2, sq+1):
if is_prime[i]:
for j in range(i*i, n+1, i):
is_prime[j] = 0
return [i for i, p in enumerate(is_prime) if p]
from collections import defaultdict
def factorize(n):
res = defaultdict(int)
d = 2
while d * d <= n:
while n % d == 0:
res[d] += 1
n //= d
d += 1
if n > 1:
res[n] += 1
return res
L, H = map(int, input().split())
primes = sieve(10**5 + 100)
for p in primes[::-1]:
l = (L+p-1) // p
h = H // p
for d in reversed(range(max(p, l), h+1)):
fac = factorize(d)
if min(list(fac.keys()) + [p]) == p:
res = p * d
print(res)
exit(0)