結果
問題 |
No.2604 Initial Motion
|
ユーザー |
|
提出日時 | 2024-01-13 01:34:39 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,030 bytes |
コンパイル時間 | 190 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 83,860 KB |
最終ジャッジ日時 | 2024-09-28 01:01:14 |
合計ジャッジ時間 | 15,126 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 11 WA * 23 TLE * 5 |
ソースコード
# https://tjkendev.github.io/procon-library/python/min_cost_flow/primal-dual.html # 最小費用流(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())) MCF = MinCostFlow(N+2) INF = 10**16 for _ in range(M): u, v, d = map(int, input().split()) MCF.addEdge(u, v, INF, d) MCF.addEdge(v, u, INF, d) for a in A: MCF.addEdge(0, a, 1, 0) for i in range(N): MCF.addEdge(i+1, N+1, B[i], 0) ans = MCF.minCostFlow(0, N+1, K) print(ans)