結果
| 問題 |
No.2531 Coloring Vertices on Namori
|
| コンテスト | |
| ユーザー |
navel_tos
|
| 提出日時 | 2023-11-04 16:12:50 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 602 ms / 2,000 ms |
| コード長 | 2,353 bytes |
| コンパイル時間 | 319 ms |
| コンパイル使用メモリ | 82,532 KB |
| 実行使用メモリ | 160,400 KB |
| 最終ジャッジ日時 | 2024-09-25 22:06:45 |
| 合計ジャッジ時間 | 14,976 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 31 |
ソースコード
#yukicoder411G
#UnionFind
class UnionFind:
def __init__(self,N): self._parent=[-1 for i in[0]*N]
def find(self,v): #頂点vの親を探し、経路圧縮する
vertices=[]
while self._parent[v]>=0: vertices.append(v);v=self._parent[v]
for i in vertices: self._parent[i]=v
return v
def unite(self,x,y): #頂点xとyを併合し、併合の有無を返す
x,y = self.find(x),self.find(y)
if x==y: return 0
if self._parent[x]>self._parent[y]: x,y=y,x #負値で管理
self._parent[x]+=self._parent[y]; self._parent[y]=x; return 1
def same(self,x,y):return self.find(x)==self.find(y) #xとyは同一集合か返す
def size(self,x): return -self._parent[self.find(x)] #xの集合のサイズを求める
#入力受取 UFで頂点を結合するが、「ここを追加すると閉路になる」という辺のみ除外する
N,K=map(int,input().split()); G=[[] for _ in range(N)]
UF=UnionFind(N); MOD=998244353; ExS,ExG=0,0
for _ in range(N):
u,v=map(lambda x: int(x)-1,input().split())
if UF.same(u,v): ExS,ExG=u,v
else: UF.unite(u,v); G[u].append(v); G[v].append(u)
#なもりグラフのサイクルを検知 P[i]: 頂点iの親 とすれば復元可能
P=[-1]*N; visit=[(ExS,-1)]
for now,back in visit:
P[now]=back
for next in G[now]:
if next!=back: visit.append((next,now))
cycle=[ExG]
for now in cycle:
if now==ExS: break
cycle.append(P[now])
#ようやく頂点を追加してグラフを完成させる
G[ExS].append(ExG); G[ExG].append(ExS)
#DP[i]: 頂点iを特定の色で塗ると決めたときの場合の数
CS=set(cycle)
DP=[1 for _ in range(N)]
for i in cycle:
visit=[(i,i)]
for now,back in visit:
for next in G[now]:
if next!=back and next not in CS: visit.append((next,now))
for now,back in visit[::-1]:
if now in CS: continue
DP[back]=DP[back]*DP[now]%MOD*(K-1)%MOD
#sDP[i][x]: サイクルのi番目の要素について、色0で塗るならx=0
sDP=[[1]*2 for _ in range(len(cycle)+1)]
sDP[0][0]=1; sDP[0][1]=0 #cycle[-1]を色0で塗った、と仮定する
for i in range(len(cycle)):
sDP[i+1][0]=sDP[i][1]*DP[cycle[i]]%MOD
sDP[i+1][1]=( sDP[i][0]*(K-1)%MOD+sDP[i][1]*(K-2) )%MOD*DP[cycle[i]]%MOD
#答えを出力
print(sDP[-1][0]*K%MOD)
navel_tos