結果

問題 No.1479 Matrix Eraser
ユーザー tamatotamato
提出日時 2021-04-16 20:31:25
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 3,520 bytes
コンパイル時間 265 ms
コンパイル使用メモリ 87,488 KB
実行使用メモリ 204,584 KB
最終ジャッジ日時 2023-09-15 21:21:33
合計ジャッジ時間 6,827 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 105 ms
85,192 KB
testcase_01 AC 107 ms
80,820 KB
testcase_02 AC 108 ms
80,644 KB
testcase_03 AC 107 ms
81,276 KB
testcase_04 AC 110 ms
80,832 KB
testcase_05 AC 108 ms
81,084 KB
testcase_06 AC 108 ms
80,900 KB
testcase_07 TLE -
testcase_08 TLE -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

mod = 1000000007
eps = 10**-9


def main():
    import sys
    input = sys.stdin.readline
    from collections import deque

    class Dinic:
        def __init__(self, N):
            self.N = N
            self.adj = [[] for _ in range(N + 1)]
            self.lv = None
            self.progress = None
            self.inf = 1 << 60

        def add_edge(self, fr, to, cap):
            # [to, cap, rev]
            forward = [to, cap, None]
            backward = forward[2] = [fr, 0, forward]
            self.adj[fr].append(forward)
            self.adj[to].append(backward)

        def bfs(self, s, t):
            que = deque([s])
            lv = [-1] * (self.N + 1)
            lv[s] = 0
            while que:
                v = que.popleft()
                lvv = lv[v]
                for u, cap, _ in self.adj[v]:
                    if cap and lv[u] == -1:
                        que.append(u)
                        lv[u] = lvv + 1
                        if u == t:
                            break
            self.lv = lv

        def dfs(self, s, t, flow):
            st = deque([s])
            st_f = deque([flow])
            while st:
                v = st[-1]
                f = st_f[-1]
                if v == t:
                    for i in range(len(st) - 1, 0, -1):
                        p = st[i - 1]
                        edge = self.adj[p][self.progress[p] - 1]
                        edge[1] -= f
                        edge[2][1] += f
                    return f
                if self.progress[v] == len(self.adj[v]):
                    st.pop()
                    st_f.pop()
                    continue
                else:
                    adj_v = self.adj[v]
                    for i in range(self.progress[v], len(adj_v)):
                        i = self.progress[v]
                        self.progress[v] = i + 1
                        u, cap, rev = adj_v[i]
                        if not cap or self.lv[u] != self.lv[v] + 1:
                            continue
                        else:
                            st.append(u)
                            st_f.append(min(st_f[-1], cap))
                            break

        def max_flow(self, s, t):
            flow = 0
            while True:
                self.bfs(s, t)
                if self.lv[t] < 0:
                    return flow
                self.progress = [0] * (self.N + 1)
                flow_new = self.dfs(s, t, self.inf)
                while flow_new:
                    flow += flow_new
                    flow_new = self.dfs(s, t, self.inf)

    H, W = map(int, input().split())
    A = []
    for _ in range(H):
        A.append(list(map(int, input().split())))

    D = Dinic(H + W + 2)
    s = H + W + 1
    t = s + 1
    for h in range(H):
        D.add_edge(s, h+1, 1)
    for w in range(W):
        D.add_edge(H + w + 1, t, 1)
    edge = [[] for _ in range(5 * 10 ** 5 + 1)]
    for h in range(H):
        for w in range(W):
            a = A[h][w]
            edge[a].append((h, w))
    ans = 0
    for a in range(5 * 10 ** 5, 0, -1):
        if edge[a]:
            for h, w in edge[a]:
                D.add_edge(h+1, H+w+1, 1)
            ans += D.max_flow(s, t)

            D = Dinic(H + W + 2)
            s = H + W + 1
            t = s + 1
            for h in range(H):
                D.add_edge(s, h + 1, 1)
            for w in range(W):
                D.add_edge(H + w + 1, t, 1)
    print(ans)


if __name__ == '__main__':
    main()
0