結果
| 問題 | No.1301 Strange Graph Shortest Path | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2021-02-15 19:55:23 | 
| 言語 | PyPy3 (7.3.15) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 1,335 ms / 3,000 ms | 
| コード長 | 1,564 bytes | 
| コンパイル時間 | 325 ms | 
| コンパイル使用メモリ | 82,304 KB | 
| 実行使用メモリ | 218,884 KB | 
| 最終ジャッジ日時 | 2024-07-23 11:55:14 | 
| 合計ジャッジ時間 | 40,053 ms | 
| ジャッジサーバーID (参考情報) | judge4 / judge2 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 2 | 
| other | AC * 33 | 
ソースコード
# 最小費用流
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
n,m=map(int,input().split())
abcd=[list(map(int,input().split())) for _ in range(m)]
mcf=MinCostFlow(n)
for a,b,c,d in abcd:
  a,b=a-1,b-1
  mcf.add_edge(a,b,1,c)
  mcf.add_edge(b,a,1,c)
  mcf.add_edge(a,b,1,d)
  mcf.add_edge(b,a,1,d)
print(mcf.flow(0,n-1,2))
            
            
            
        