結果
| 問題 |
No.2497 GCD of LCMs
|
| コンテスト | |
| ユーザー |
mymelochan
|
| 提出日時 | 2023-10-06 23:50:25 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 402 ms / 2,000 ms |
| コード長 | 1,598 bytes |
| コンパイル時間 | 174 ms |
| コンパイル使用メモリ | 81,856 KB |
| 実行使用メモリ | 84,564 KB |
| 最終ジャッジ日時 | 2024-07-26 17:23:19 |
| 合計ジャッジ時間 | 3,919 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 14 |
ソースコード
#############################################################
import sys
sys.setrecursionlimit(10**7)
from heapq import heappop,heappush
from collections import deque,defaultdict,Counter
from bisect import bisect_left, bisect_right
from itertools import product,combinations,permutations
ipt = sys.stdin.readline
def iin():
return int(ipt())
def lmin():
return list(map(int,ipt().split()))
def prime_factorize(n):
d = defaultdict(int)
while n % 2 == 0:
d[2] += 1
n //= 2
f = 3
while f * f <= n:
if n % f == 0:
d[f] += 1
n //= f
else:
f += 2
if n != 1:
d[n] += 1
return d
MOD = 998244353
#############################################################
N,M = lmin()
A = lmin()
D = [prime_factorize(a) for a in A]
S = set()
for d in D:
for k in d.keys():
S.add(k)
G = [[] for _ in range(N)]
for _ in range(M):
u,v = lmin()
u,v = u-1,v-1
G[u].append(v)
G[v].append(u)
cost = defaultdict(list)
for s in S:
tmp_cost = [MOD]*N
tmp_cost[0] = D[0][s]
hq = [(tmp_cost[0],0)]
while hq:
c1,cur = heappop(hq)
if tmp_cost[cur] < c1:
continue
tmp_cost[cur] = c1
for nxt in G[cur]:
c2 = D[nxt][s]
v = max(c1,c2)
if v < tmp_cost[nxt]:
tmp_cost[nxt] = v
heappush(hq,(v,nxt))
cost[s] = tmp_cost
#print(cost)
for i in range(N):
ans = 1
for s in S:
ans *= pow(s,cost[s][i],MOD)
ans %= MOD
print(ans)
mymelochan