結果
問題 |
No.2604 Initial Motion
|
ユーザー |
![]() |
提出日時 | 2024-01-12 22:55:16 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,978 bytes |
コンパイル時間 | 336 ms |
コンパイル使用メモリ | 82,208 KB |
実行使用メモリ | 89,028 KB |
最終ジャッジ日時 | 2024-09-27 23:42:42 |
合計ジャッジ時間 | 14,908 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 10 WA * 3 TLE * 1 -- * 25 |
ソースコード
# 最小費用流(minimum cost flow) class MinCostFlow: def __init__(self, n): self.n = n self.G = [[] for i in range(n)] def addEdge(self, f, t, cap, cost): # [to, cap, cost, rev] self.G[f].append([t, cap, cost, len(self.G[t])]) self.G[t].append([f, 0, -cost, len(self.G[f])-1]) def minCostFlow(self, s, t, f): n = self.n G = self.G prevv = [0]*n; preve = [0]*n INF = 10**9+7 res = 0 while f: dist = [INF]*n dist[s] = 0 update = 1 while update: update = 0 for v in range(n): if dist[v] == INF: continue gv = G[v] for i in range(len(gv)): to, cap, cost, rev = gv[i] if cap > 0 and dist[v] + cost < dist[to]: dist[to] = dist[v] + cost prevv[to] = v; preve[to] = i update = 1 if dist[t] == INF: return -1 d = f; v = t while v != s: d = min(d, G[prevv[v]][preve[v]][1]) v = prevv[v] f -= d res += d * dist[t] v = t while v != s: e = G[prevv[v]][preve[v]] e[1] -= d G[v][e[3]][1] += d v = prevv[v] return res import sys input=sys.stdin.readline K,N,M=map(int,input().split()) A=list(map(int,input().split())) B=list(map(int,input().split())) #初期化 mf=MinCostFlow(N+2) for i in range(K): mf.addEdge(0,A[i],1,0) for i in range(N): mf.addEdge(i+1,N+1,B[i],0) for i in range(M): u,v,d=map(int,input().split()) mf.addEdge(u,v,4000,d) mf.addEdge(v,u,4000,d) #flowを流す。返り値は(流れの量,費用)。O(F(N+M)log(N+M)) print(mf.minCostFlow(0,N+1,K))