結果
問題 | 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 sysinput = sys.stdin.readlineN = int(input())mod = 998244353class sieve:def __init__(self, n):self.n = nself.sv = [1] * (n + 1)self.sv[0] = 0self.sv[1] = 0for i in range(2, n + 1):if self.sv[i]:for j in range(i * 2, n + 1, i):self.sv[j] = 0def isprime(self, x):if x > self.n:return Falsereturn self.sv[x] == 1def factorize(self, x):res = []for i in range(2, int(x ** 0.5) + 1):if self.sv[i]:while x % i == 0:x //= ires.append(i)if x != 1:res.append(x)return resdef 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)continueres[j] = max(f.count(j), res[j])rres = 1for i in range(self.n + 1):if res[i] != 0:rres *= pow(i, res[i], mod)rres %= modfor i in ex:rres *= irres %= modreturn rressv = sieve(N)res = 1for p in range(2, N + 1):if sv.isprime(p) == False: continuec = 0for x in range(p, N, p):y = x * (N - x)t = 0while y % p == 0:y //= pt += 1c = max(c, t)res *= pow(p, c, mod)res %= modprint(res)