結果

問題 No.957 植林
ユーザー lam6er
提出日時 2025-04-16 16:11:22
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 3,220 bytes
コンパイル時間 173 ms
コンパイル使用メモリ 81,324 KB
実行使用メモリ 166,260 KB
最終ジャッジ日時 2025-04-16 16:16:36
合計ジャッジ時間 17,104 ms
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 15 TLE * 1 -- * 29
権限があれば一括ダウンロードができます

ソースコード

diff #

from collections import deque

class Edge:
    def __init__(self, to, rev, residual):
        self.to = to
        self.rev = rev
        self.residual = residual

class Dinic:
    def __init__(self, n):
        self.size = n
        self.graph = [[] for _ in range(n)]
    
    def add_edge(self, fr, to, cap):
        forward = Edge(to, len(self.graph[to]), cap)
        backward = Edge(fr, len(self.graph[fr]), 0)
        self.graph[fr].append(forward)
        self.graph[to].append(backward)
    
    def bfs_level(self, s, t, level):
        q = deque()
        level[:] = [-1] * self.size
        level[s] = 0
        q.append(s)
        while q:
            v = q.popleft()
            for edge in self.graph[v]:
                if edge.residual > 0 and level[edge.to] == -1:
                    level[edge.to] = level[v] + 1
                    q.append(edge.to)
                    if edge.to == t:
                        return
    
    def dfs_flow(self, v, t, upTo, iter_, level):
        if v == t:
            return upTo
        for i in range(iter_[v], len(self.graph[v])):
            edge = self.graph[v][i]
            if edge.residual > 0 and level[v] < level[edge.to]:
                d = self.dfs_flow(edge.to, t, min(upTo, edge.residual), iter_, level)
                if d > 0:
                    edge.residual -= d
                    self.graph[edge.to][edge.rev].residual += d
                    return d
            iter_[v] += 1
        return 0
    
    def max_flow(self, s, t):
        flow = 0
        level = [-1] * self.size
        while True:
            self.bfs_level(s, t, level)
            if level[t] == -1:
                return flow
            iter_ = [0] * self.size
            while True:
                f = self.dfs_flow(s, t, float('inf'), iter_, level)
                if f == 0:
                    break
                flow += f
            level = [-1] * self.size
        return flow

def main():
    import sys
    input = sys.stdin.read().split()
    idx = 0
    H = int(input[idx]); idx +=1
    W = int(input[idx]); idx +=1
    
    G = []
    for _ in range(H):
        row = list(map(int, input[idx:idx+W]))
        idx += W
        G.append(row)
    
    R = list(map(int, input[idx:idx+H]))
    idx += H
    C = list(map(int, input[idx:idx+W]))
    idx += W
    
    total_R = sum(R)
    total_C = sum(C)
    INF = 1 << 60
    
    num_nodes = 2 + H + W + H * W
    dinic = Dinic(num_nodes)
    
    # Add edges from S to R_i
    for i in range(H):
        dinic.add_edge(0, 2 + i, R[i])
    
    # Add edges from S to C_j
    for j in range(W):
        dinic.add_edge(0, 2 + H + j, C[j])
    
    # Add edges for X_{i,j}
    for i in range(H):
        for j in range(W):
            x_node = 2 + H + W + i * W + j
            dinic.add_edge(x_node, 1, G[i][j])
            # Add edges from R_i to X_{i,j}
            r_node = 2 + i
            dinic.add_edge(r_node, x_node, INF)
            # Add edges from C_j to X_{i,j}
            c_node = 2 + H + j
            dinic.add_edge(c_node, x_node, INF)
    
    max_flow = dinic.max_flow(0, 1)
    profit = total_R + total_C - max_flow
    print(max(profit, 0))

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