#592722 (Python3) No.421 しろくろチョコレート

提出ソース

結果

問題	No.421 しろくろチョコレート
ユーザー	tktk_snsn
提出日時	2020-12-13 12:18:50
言語	Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果	AC
実行時間	69 ms / 2,000 ms
コード長	4,113 bytes
コンパイル時間	348 ms
コンパイル使用メモリ	12,444 KB
実行使用メモリ	13,248 KB
最終ジャッジ日時	2023-10-24 15:03:16
合計ジャッジ時間	4,487 ms
ジャッジサーバーID （参考情報）	judge15 / judge12

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	32 ms 10,556 KB
testcase_01	AC	44 ms 11,596 KB
testcase_02	AC	36 ms 10,960 KB
testcase_03	AC	33 ms 10,616 KB
testcase_04	AC	35 ms 10,820 KB
testcase_05	AC	36 ms 10,812 KB
testcase_06	AC	33 ms 10,632 KB
testcase_07	AC	40 ms 11,180 KB
testcase_08	AC	36 ms 10,964 KB
testcase_09	AC	33 ms 10,700 KB
testcase_10	AC	34 ms 10,696 KB
testcase_11	AC	41 ms 11,232 KB
testcase_12	AC	32 ms 10,548 KB
testcase_13	AC	32 ms 10,548 KB
testcase_14	AC	32 ms 10,548 KB
testcase_15	AC	32 ms 10,548 KB
testcase_16	AC	56 ms 12,180 KB
testcase_17	AC	58 ms 12,212 KB
testcase_18	AC	57 ms 12,252 KB
testcase_19	AC	32 ms 10,548 KB
testcase_20	AC	49 ms 11,756 KB
testcase_21	AC	53 ms 11,936 KB
testcase_22	AC	42 ms 11,268 KB
testcase_23	AC	31 ms 10,548 KB
testcase_24	AC	31 ms 10,548 KB
testcase_25	AC	32 ms 10,548 KB
testcase_26	AC	31 ms 10,548 KB
testcase_27	AC	32 ms 10,548 KB
testcase_28	AC	42 ms 11,584 KB
testcase_29	AC	44 ms 11,792 KB
testcase_30	AC	44 ms 11,624 KB
testcase_31	AC	58 ms 13,248 KB
testcase_32	AC	42 ms 11,684 KB
testcase_33	AC	45 ms 11,704 KB
testcase_34	AC	34 ms 10,868 KB
testcase_35	AC	34 ms 10,880 KB
testcase_36	AC	40 ms 11,260 KB
testcase_37	AC	58 ms 12,196 KB
testcase_38	AC	65 ms 12,556 KB
testcase_39	AC	38 ms 11,132 KB
testcase_40	AC	45 ms 11,600 KB
testcase_41	AC	35 ms 10,836 KB
testcase_42	AC	35 ms 10,748 KB
testcase_43	AC	39 ms 11,208 KB
testcase_44	AC	50 ms 12,056 KB
testcase_45	AC	35 ms 10,828 KB
testcase_46	AC	35 ms 10,800 KB
testcase_47	AC	43 ms 11,468 KB
testcase_48	AC	54 ms 12,592 KB
testcase_49	AC	37 ms 11,112 KB
testcase_50	AC	35 ms 10,900 KB
testcase_51	AC	51 ms 11,984 KB
testcase_52	AC	40 ms 11,188 KB
testcase_53	AC	34 ms 10,756 KB
testcase_54	AC	37 ms 11,004 KB
testcase_55	AC	35 ms 10,708 KB
testcase_56	AC	34 ms 10,700 KB
testcase_57	AC	36 ms 10,800 KB
testcase_58	AC	46 ms 11,564 KB
testcase_59	AC	38 ms 10,924 KB
testcase_60	AC	69 ms 12,708 KB
testcase_61	AC	65 ms 12,776 KB
testcase_62	AC	33 ms 10,556 KB
testcase_63	AC	32 ms 10,556 KB
testcase_64	AC	33 ms 10,744 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

from collections import deque
import sys
input = sys.stdin.readline
sys.setrecursionlimit(10 ** 7)
DX = (-1, 0, 1, 0, -1, -1, 1, 1)
DY = (0, 1, 0, -1, -1, 1, -1, 1)
DX = DX[:4]
DY = DY[:4]


class MF_graph(object):
    def __init__(self, n):
        self.n = n
        self.g = [[] for _ in range(n)]  # to, rev, cap
        self.pos = []

    def add_edge(self, frm, to, cap):
        """ frmからtoへ容量cap, 流量0の辺を追加, 何番目の辺かを返す"""
        m = len(self.pos)
        self.pos.append((frm, len(self.g[frm])))
        self.g[frm].append([to, len(self.g[to]), cap])
        self.g[to].append([frm, len(self.g[frm]) - 1, 0])
        return m

    def get_edge(self, i):
        e_to, e_rev, e_cap = self.g[self.pos[i][0]][self.pos[i][1]]
        re_to, _, re_cap = self.g[e_to][e_rev]
        # from, to, cap, flow
        return (re_to, e_to, e_cap + re_cap, re_cap)

    def edges(self):
        m = len(self.pos)
        for i in range(m):
            yield self.get_edge(i)

    def change_edge(self, i, new_cap, new_flow):
        """ i番目に追加された辺の容量,流量をnew_cap, new_flowに変更する """
        f, s = self.pos[i]
        rf, rs, _ = self.g[f][s]
        self.g[f][s][2] = new_cap - new_flow
        self.g[rf][rs][2] = new_flow
        return

    def __dfs(self, s, v, up):
        if v == s:
            return up
        res = 0
        level_v = self.level[v]
        for i in range(self.iter[v], len(self.g[v])):
            u_to, u_rev, _ = self.g[v][i]
            if level_v <= self.level[u_to] or self.g[u_to][u_rev][2] == 0:
                continue
            d = self.__dfs(s, u_to, min(up - res, self.g[u_to][u_rev][2]))
            if d <= 0:
                continue
            self.g[v][i][2] += d
            self.g[u_to][u_rev][2] -= d
            res += d
            if res == up:
                break
        return res

    def flow(self, s, t, flow_limit=10**18):
        """ sからtへflow_limitだけ流す。流せた量を返す """
        self.iter = [0] * self.n
        flow = 0
        while flow < flow_limit:
            self.level = [-1] * self.n
            self.level[s] = 0
            que = deque([s])
            while que:
                v = que.popleft()
                for u_to, _, u_cap in self.g[v]:
                    if u_cap == 0 or self.level[u_to] >= 0:
                        continue
                    self.level[u_to] = self.level[v] + 1
                    if u_to == t:
                        break
                    que.append(u_to)

            if self.level[t] == -1:
                break
            self.iter = [0] * self.n
            while flow < flow_limit:
                f = self.__dfs(s, t, flow_limit - flow)
                if not f:
                    break
                flow += f
        return flow

    def min_cut(self, s):
        """ 最小カットを返す。sから到達可能だとTrue """
        visited = [False] * self.n
        que = deque([s])
        while que:
            v = que.popleft()
            visited[v] = True
            for u_to, _, u_cap in self.g[v]:
                if u_cap and (not visited[u_to]):
                    visited[u_to] = True
                    que.append(u_to)
        return visited


N, M = map(int, input().split())
S = [input().rstrip() for _ in range(N)]
source = N * M
sink = source + 1

G = MF_graph(N * M + 2)
black = 0
white = 0
for x in range(N):
    for y in range(M):
        if S[x][y] == ".":
            continue
        if (x + y) % 2 == 0:
            white += 1
            G.add_edge(source, x * M + y, 1)
            for dx, dy in zip(DX, DY):
                nx = x + dx
                ny = y + dy
                if 0 <= nx < N and 0 <= ny < M and S[nx][ny] != ".":
                    G.add_edge(x * M + y, nx * M + ny, 1)
        else:
            black += 1
            G.add_edge(x * M + y, sink, 1)

ans = 0
flow = G.flow(source, sink)
ans += 100 * flow

black -= flow
white -= flow
d = min(black, white)
ans += 10 * d

black -= d
white -= d
ans += black + white

print(ans)

yukicoder

結果

テストケース

ソースコード