結果

問題 No.2604 Initial Motion
コンテスト
ユーザー rlangevin
提出日時 2024-01-19 08:57:30
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 2,261 bytes
コンパイル時間 222 ms
コンパイル使用メモリ 82,576 KB
実行使用メモリ 89,920 KB
最終ジャッジ日時 2024-09-28 03:20:20
合計ジャッジ時間 7,514 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10 TLE * 1 -- * 28
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import sys
input = sys.stdin.readline

class MinCostFlow:
    def __init__(self, N):
        self.N = N
        self.inf = 10 ** 18
        self.G = [[] for i in range(N)]

    def add_edge(self, u, v, cap, cost):
        self.G[u].append((v, cap, cost, len(self.G[v])))
        self.G[v].append((u, 0, -cost, len(self.G[u]) - 1))

    def bellman_ford(self, s):
        dist = [self.inf] * self.N
        dist[s] = 0
        pv = [0] * self.N
        pe = [0] * self.N
        while True:
            update = False
            for v in range(self.N):
                if dist[v] == self.inf:
                    continue
                for i in range(len(self.G[v])):
                    next, cap, cost, _ = self.G[v][i]
                    if cap > 0 and dist[next] > dist[v] + cost:
                        dist[next] = dist[v] + cost
                        update = True
                        pv[next] = v
                        pe[next] = i
            if not update:
                break
        return dist, pv, pe

    def calc_min_cost_flow(self, s, t, f):
        result = 0
        while f > 0:
            dist, pv, pe = self.bellman_ford(s)
            if dist[t] == self.inf:
                return self.inf
            flow = f
            v = t
            while v != s:
                flow = min(flow, self.G[pv[v]][pe[v]][1])
                v = pv[v]
            result += flow * dist[t]
            f -= flow
            v = t
            while v != s:
                d, cap, cost, r = self.G[pv[v]][pe[v]]
                cap -= flow
                self.G[pv[v]][pe[v]] = (d, cap, cost, r)
                rev = self.G[pv[v]][pe[v]][3]
                d, cap, cost, r = self.G[v][rev]
                cap += flow
                self.G[v][rev] = (d, cap, cost, r)
                v = pv[v]
        return result
    
    
K, N, M = map(int, input().split())
A = list(map(int, input().split()))
B = list(map(int, input().split()))
G = MinCostFlow(N + 2)
s, g = 0, N + 1
for i in range(M):
    u, v, d = map(int, input().split())
    G.add_edge(u, v, K, d)
    G.add_edge(v, u, K, d)
    
for i in range(K):
    G.add_edge(s, A[i], 1, 0)
for i in range(N):
    G.add_edge(i + 1, g, B[i], 0)
    
print(G.calc_min_cost_flow(s, g, K))
0