結果

問題 No.1301 Strange Graph Shortest Path
ユーザー persimmon-persimmonpersimmon-persimmon
提出日時 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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
53,120 KB
testcase_01 AC 38 ms
52,992 KB
testcase_02 AC 1,098 ms
209,564 KB
testcase_03 AC 1,041 ms
193,208 KB
testcase_04 AC 1,262 ms
214,408 KB
testcase_05 AC 1,100 ms
211,484 KB
testcase_06 AC 1,091 ms
203,828 KB
testcase_07 AC 1,107 ms
205,648 KB
testcase_08 AC 1,079 ms
195,380 KB
testcase_09 AC 1,044 ms
196,712 KB
testcase_10 AC 990 ms
194,728 KB
testcase_11 AC 1,335 ms
207,800 KB
testcase_12 AC 1,233 ms
208,032 KB
testcase_13 AC 1,130 ms
210,512 KB
testcase_14 AC 1,089 ms
196,724 KB
testcase_15 AC 1,099 ms
198,820 KB
testcase_16 AC 1,210 ms
214,236 KB
testcase_17 AC 1,142 ms
213,516 KB
testcase_18 AC 1,216 ms
200,432 KB
testcase_19 AC 1,135 ms
205,008 KB
testcase_20 AC 1,091 ms
202,240 KB
testcase_21 AC 1,113 ms
210,500 KB
testcase_22 AC 1,229 ms
206,024 KB
testcase_23 AC 1,105 ms
211,524 KB
testcase_24 AC 1,125 ms
203,836 KB
testcase_25 AC 1,263 ms
214,872 KB
testcase_26 AC 1,102 ms
205,184 KB
testcase_27 AC 1,161 ms
207,196 KB
testcase_28 AC 981 ms
205,232 KB
testcase_29 AC 1,146 ms
213,116 KB
testcase_30 AC 1,282 ms
213,796 KB
testcase_31 AC 1,197 ms
212,420 KB
testcase_32 AC 39 ms
53,144 KB
testcase_33 AC 546 ms
204,416 KB
testcase_34 AC 1,251 ms
218,884 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

# 最小費用流
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))
0