結果

問題 No.421 しろくろチョコレート
ユーザー mkawa2mkawa2
提出日時 2020-04-21 18:44:50
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 217 ms / 2,000 ms
コード長 3,962 bytes
コンパイル時間 257 ms
コンパイル使用メモリ 12,452 KB
実行使用メモリ 12,668 KB
最終ジャッジ日時 2023-10-24 15:00:33
合計ジャッジ時間 4,560 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
10,664 KB
testcase_01 AC 42 ms
11,400 KB
testcase_02 AC 34 ms
10,900 KB
testcase_03 AC 33 ms
10,704 KB
testcase_04 AC 34 ms
10,796 KB
testcase_05 AC 35 ms
10,860 KB
testcase_06 AC 33 ms
10,720 KB
testcase_07 AC 40 ms
11,148 KB
testcase_08 AC 34 ms
10,884 KB
testcase_09 AC 33 ms
10,768 KB
testcase_10 AC 33 ms
10,756 KB
testcase_11 AC 38 ms
11,180 KB
testcase_12 AC 31 ms
10,660 KB
testcase_13 AC 31 ms
10,660 KB
testcase_14 AC 31 ms
10,660 KB
testcase_15 AC 31 ms
10,660 KB
testcase_16 AC 48 ms
11,916 KB
testcase_17 AC 55 ms
11,940 KB
testcase_18 AC 58 ms
11,968 KB
testcase_19 AC 31 ms
10,660 KB
testcase_20 AC 45 ms
11,540 KB
testcase_21 AC 53 ms
11,728 KB
testcase_22 AC 41 ms
11,204 KB
testcase_23 AC 32 ms
10,660 KB
testcase_24 AC 32 ms
10,660 KB
testcase_25 AC 31 ms
10,660 KB
testcase_26 AC 31 ms
10,660 KB
testcase_27 AC 31 ms
10,660 KB
testcase_28 AC 36 ms
11,356 KB
testcase_29 AC 38 ms
11,548 KB
testcase_30 AC 39 ms
11,400 KB
testcase_31 AC 48 ms
12,668 KB
testcase_32 AC 37 ms
11,456 KB
testcase_33 AC 43 ms
11,472 KB
testcase_34 AC 33 ms
10,764 KB
testcase_35 AC 33 ms
10,776 KB
testcase_36 AC 38 ms
11,180 KB
testcase_37 AC 65 ms
11,928 KB
testcase_38 AC 72 ms
12,164 KB
testcase_39 AC 35 ms
11,056 KB
testcase_40 AC 74 ms
11,452 KB
testcase_41 AC 34 ms
10,872 KB
testcase_42 AC 34 ms
10,804 KB
testcase_43 AC 36 ms
11,132 KB
testcase_44 AC 217 ms
11,840 KB
testcase_45 AC 34 ms
10,828 KB
testcase_46 AC 35 ms
10,844 KB
testcase_47 AC 42 ms
11,308 KB
testcase_48 AC 50 ms
12,228 KB
testcase_49 AC 35 ms
11,016 KB
testcase_50 AC 33 ms
10,884 KB
testcase_51 AC 67 ms
11,784 KB
testcase_52 AC 38 ms
11,124 KB
testcase_53 AC 33 ms
10,780 KB
testcase_54 AC 34 ms
10,920 KB
testcase_55 AC 33 ms
10,768 KB
testcase_56 AC 34 ms
10,760 KB
testcase_57 AC 34 ms
10,828 KB
testcase_58 AC 42 ms
11,408 KB
testcase_59 AC 37 ms
10,932 KB
testcase_60 AC 93 ms
12,284 KB
testcase_61 AC 66 ms
12,324 KB
testcase_62 AC 32 ms
10,660 KB
testcase_63 AC 32 ms
10,660 KB
testcase_64 AC 31 ms
10,664 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from itertools import *
from bisect import *
from collections import *
from heapq import *
import sys

sys.setrecursionlimit(10 ** 6)

def II(): return int(sys.stdin.readline())
def MI(): return map(int, sys.stdin.readline().split())
def LI(): return list(map(int, sys.stdin.readline().split()))
def SI(): return sys.stdin.readline()[:-1]
def LLI(rows_number): return [LI() for _ in range(rows_number)]
def LLI1(rows_number): return [LI1() for _ in range(rows_number)]
int1 = lambda x: int(x) - 1
def MI1(): return map(int1, sys.stdin.readline().split())
def LI1(): return list(map(int1, sys.stdin.readline().split()))
p2D = lambda x: print(*x, sep="\n")
dij = [(1, 0), (0, 1), (-1, 0), (0, -1)]

class Dinic():
    def __init__(self, n, s, t):
        self.n, self.s, self.t = n, s, t
        self.to = defaultdict(list)
        self.level = [-1]
        self.max_flow = -1

    def add_edge(self, u, v, cap):
        u_index_in_to_v = len(self.to[v])
        v_index_in_to_u = len(self.to[u])
        self.to[u].append([v, cap, u_index_in_to_v])
        self.to[v].append([u, 0, v_index_in_to_u])

    # 無向辺の追加
    def add_undir_edge(self, u, v, cap):
        u_index_in_to_v = len(self.to[v])
        v_index_in_to_u = len(self.to[u])
        self.to[u].append([v, cap, u_index_in_to_v])
        self.to[v].append([u, cap, v_index_in_to_u])

    def set_level(self):
        s = self.s
        level = [-1] * self.n
        level[s] = 0
        q = deque()
        q.append([s, 0])
        while q:
            u, u_level = q.popleft()
            for v, cap, _ in self.to[u]:
                if cap == 0: continue
                if level[v] != -1: continue
                level[v] = u_level + 1
                if v == self.t:
                    self.level = level
                    return True
                q.append([v, u_level + 1])
        return False

    def dfs(self, u=-1, flow_to_me=10 ** 16):
        if u == -1: u = self.s
        if u == self.t: return flow_to_me
        flow_from_me = 0
        u_level = self.level[u]
        for utov_i, (v, cap, vtou_i) in enumerate(self.to[u]):
            if self.level[v] != u_level + 1: continue
            if cap == 0: continue
            flow_to_v = self.dfs(v, min(cap, flow_to_me - flow_from_me))
            if not flow_to_v: continue
            flow_from_me += flow_to_v
            self.to[u][utov_i][1] -= flow_to_v
            self.to[v][vtou_i][1] += flow_to_v
        return flow_from_me

    def calculation(self):
        res = 0
        while self.set_level():
            res += self.dfs()
        return res

    # これが出力用
    def get_max_flow(self):
        if self.max_flow == -1:
            self.max_flow = self.calculation()
        return self.max_flow

def main():
    h,w=MI()
    ss=[SI() for _ in range(h)]
    aa=[[0]*w for _ in range(h)]
    # 始点を頂点0、終点を頂点1とするので、W,Bに2から頂点番号を振っていく
    cnt=2
    for i,j in product(range(h),range(w)):
        if ss[i][j]==".":continue
        aa[i][j]=cnt
        cnt+=1
    # 最大流のライブラリ
    mf=Dinic(cnt,0,1)
    # それぞれのマスから右と下に辺を張るとともにw,bの数を数える
    wc=bc=0
    for i,j in product(range(h),range(w)):
        a=aa[i][j]
        if a==0:continue
        # 黒のとき（2部グラフの右側）
        if (i+j)%2:
            bc+=1
            mf.add_edge(a,1,1)
            if j+1<w and aa[i][j+1]:mf.add_edge(aa[i][j+1],a,1)
            if i+1<h and aa[i+1][j]:mf.add_edge(aa[i+1][j],a,1)
        # 白のとき（2部グラフの左側）
        else:
            wc+=1
            mf.add_edge(0,a,1)
            if j+1<w and aa[i][j+1]:mf.add_edge(a,aa[i][j+1],1)
            if i+1<h and aa[i+1][j]:mf.add_edge(a,aa[i+1][j],1)
    # ペアの組数がp
    p=mf.get_max_flow()
    wc,bc=wc-p,bc-p
    mn=min(wc,bc)
    ans=p*100+mn*10+wc+bc-2*mn
    print(ans)

main()
0