結果
問題 | No.1529 Constant Lcm |
ユーザー |
![]() |
提出日時 | 2021-06-04 21:07:39 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 258 ms / 3,000 ms |
コード長 | 1,344 bytes |
コンパイル時間 | 315 ms |
コンパイル使用メモリ | 82,112 KB |
実行使用メモリ | 83,872 KB |
最終ジャッジ日時 | 2024-11-19 11:20:25 |
合計ジャッジ時間 | 4,529 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 24 |
ソースコード
import sys input = sys.stdin.readline N = int(input()) mod = 998244353 class sieve: def __init__(self, n): self.n = n self.sv = [1] * (n + 1) self.sv[0] = 0 self.sv[1] = 0 for i in range(2, n + 1): if self.sv[i]: for j in range(i * 2, n + 1, i): self.sv[j] = 0 def isprime(self, x): if x > self.n: return False return self.sv[x] == 1 def factorize(self, x): res = [] for i in range(2, int(x ** 0.5) + 1): if self.sv[i]: while x % i == 0: x //= i res.append(i) if x != 1: res.append(x) 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]) fs = set(f) for j in fs: 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(N) res = 1 for p in range(2, N + 1): if sv.isprime(p) == False: continue c = 0 for x in range(p, N, p): y = x * (N - x) t = 0 while y % p == 0: y //= p t += 1 c = max(c, t) res *= pow(p, c, mod) res %= mod print(res)