結果
問題 | No.1611 Minimum Multiple with Double Divisors |
ユーザー |
![]() |
提出日時 | 2021-07-21 23:18:30 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,553 ms / 2,000 ms |
コード長 | 1,871 bytes |
コンパイル時間 | 339 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 77,464 KB |
最終ジャッジ日時 | 2024-10-03 02:37:40 |
合計ジャッジ時間 | 28,501 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 37 |
ソースコード
import sys input = sys.stdin.readline 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(10 ** 2) primes = [] for x in range(100): if sv.isprime(x): primes.append(x) for _ in range(int(input())): x = int(input()) res = 10 ** 15 for p in primes: if x % p: res = min(res, x * p) else: t = x c = 0 while t % p == 0: t //= p c += 1 res = min(res, x * p ** (c + 1)) #print(res, p, c) #print(res) for d in range(2, 33): z = 1. for p in primes: if p > d: break if d % p or x % p: continue cx = 0 cd = 0 t = x while t % p == 0: t //= p cx += 1 t = d while t % p == 0: t //= p cd += 1 z *= (cx + cd + 1) / (cx + 1) if z == 2.: res = min(res, x * d) print(res)