結果
| 問題 |
No.2604 Initial Motion
|
| コンテスト | |
| ユーザー |
ニックネーム
|
| 提出日時 | 2024-01-12 22:09:21 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,966 bytes |
| コンパイル時間 | 193 ms |
| コンパイル使用メモリ | 82,596 KB |
| 実行使用メモリ | 89,780 KB |
| 最終ジャッジ日時 | 2024-09-27 22:33:16 |
| 合計ジャッジ時間 | 8,288 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 13 TLE * 2 -- * 24 |
ソースコード
from heapq import heappush, heappop
class MinCostFlow:
INF = 10**18
def __init__(self, N):
self.N = N
self.G = [[] for i in range(N)]
def add_edge(self, fr, to, cap, cost):
forward = [to, cap, cost, None]
backward = forward[3] = [fr, 0, -cost, forward]
self.G[fr].append(forward)
self.G[to].append(backward)
def flow(self, s, t, f):
N = self.N; G = self.G
INF = MinCostFlow.INF
res = 0
H = [0]*N
prv_v = [0]*N
prv_e = [None]*N
d0 = [INF]*N
dist = [INF]*N
while f:
dist[:] = d0
dist[s] = 0
que = [(0, s)]
while que:
c, v = heappop(que)
if dist[v] < c:
continue
r0 = dist[v] + H[v]
for e in G[v]:
w, cap, cost, _ = e
if cap > 0 and r0 + cost - H[w] < dist[w]:
dist[w] = r = r0 + cost - H[w]
prv_v[w] = v; prv_e[w] = e
heappush(que, (r, w))
if dist[t] == INF:
return None
for i in range(N):
H[i] += dist[i]
d = f; v = t
while v != s:
d = min(d, prv_e[v][1])
v = prv_v[v]
f -= d
res += d * H[t]
v = t
while v != s:
e = prv_e[v]
e[1] -= d
e[3][1] += d
v = prv_v[v]
return res
from collections import Counter
k,n,m = map(int,input().split())
a = Counter(map(int,input().split()))
b = list(map(int,input().split()))
g = MinCostFlow(n+2)
for i,v in a.items(): g.add_edge(0,i,v,0)
for i,v in enumerate(b,1): g.add_edge(i,n+1,v,0)
for _ in range(m):
u,v,d = map(int,input().split())
g.add_edge(u,v,k,d); g.add_edge(v,u,k,d)
print(g.flow(0,n+1,k))
ニックネーム