結果
問題 | No.421 しろくろチョコレート |
ユーザー |
|
提出日時 | 2021-06-27 14:40:37 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 238 ms / 2,000 ms |
コード長 | 2,790 bytes |
コンパイル時間 | 261 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 91,792 KB |
最終ジャッジ日時 | 2024-09-23 07:31:54 |
合計ジャッジ時間 | 13,339 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 65 |
ソースコード
#!/usr/bin/env python3# from typing import *import sysimport ioimport mathimport collectionsimport decimalimport itertoolsimport bisectimport heapqdef input():return sys.stdin.readline()[:-1]# sys.setrecursionlimit(1000000)# _INPUT = """# paste here...# """# sys.stdin = io.StringIO(_INPUT)class Dinic:def __init__(self, n):self.n = nself.links = [[] for _ in range(n)]self.depth = Noneself.progress = Nonedef add_link(self, _from, to, cap):self.links[_from].append([cap, to, len(self.links[to])])self.links[to].append([0, _from, len(self.links[_from]) - 1])def bfs(self, s):depth = [-1] * self.ndepth[s] = 0q = collections.deque([s])while q:v = q.popleft()for cap, to, rev in self.links[v]:if cap > 0 and depth[to] < 0:depth[to] = depth[v] + 1q.append(to)self.depth = depthdef dfs(self, v, t, flow):if v == t:return flowlinks_v = self.links[v]for i in range(self.progress[v], len(links_v)):self.progress[v] = icap, to, rev = link = links_v[i]if cap == 0 or self.depth[v] >= self.depth[to]:continued = self.dfs(to, t, min(flow, cap))if d == 0:continuelink[0] -= dself.links[to][rev][0] += dreturn dreturn 0def max_flow(self, s, t):flow = 0while True:self.bfs(s)if self.depth[t] < 0:return flowself.progress = [0] * self.ncurrent_flow = self.dfs(s, t, INF)while current_flow > 0:flow += current_flowcurrent_flow = self.dfs(s, t, INF)INF = 10**10H, W = map(int, input().split())S = [input() for _ in range(H)]# 2分グラフ: S -> (W) -> (B) -> Tmf = Dinic(H*W+2)n_w, n_b = 0, 0for h in range(H):for w in range(W):if S[h][w] == 'w':n_w += 1mf.add_link(0, h*W+w+1, 1)if 0 < w and S[h][w-1] != '.':mf.add_link(h*W+w+1, h*W+(w-1)+1, 1)if w < W-1 and S[h][w+1] != '.':mf.add_link(h*W+w+1, h*W+(w+1)+1, 1)if 0 < h and S[h-1][w] != '.':mf.add_link(h*W+w+1, (h-1)*W+w+1, 1)if h < H-1 and S[h+1][w] != '.':mf.add_link(h*W+w+1, (h+1)*W+w+1, 1)elif S[h][w] == 'b':n_b += 1mf.add_link(h*W+w+1, H*W+1, 1)flow = mf.max_flow(0, H*W+1)n_w -= flown_b -= flowm = min(n_w, n_b)score = flow*100 + m*10 + (n_w-m) + (n_b-m)print(score)