結果

問題 No.1479 Matrix Eraser
ユーザー Today03Today03
提出日時 2024-05-23 15:00:59
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 3,092 bytes
コンパイル時間 3,015 ms
コンパイル使用メモリ 226,740 KB
実行使用メモリ 39,084 KB
最終ジャッジ日時 2024-12-20 18:49:52
合計ジャッジ時間 72,103 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 9 ms
33,512 KB
testcase_01 AC 13 ms
21,896 KB
testcase_02 AC 10 ms
22,052 KB
testcase_03 AC 10 ms
21,740 KB
testcase_04 AC 10 ms
21,748 KB
testcase_05 AC 10 ms
33,960 KB
testcase_06 AC 11 ms
21,724 KB
testcase_07 AC 1,185 ms
22,948 KB
testcase_08 AC 2,149 ms
39,084 KB
testcase_09 TLE -
testcase_10 TLE -
testcase_11 TLE -
testcase_12 AC 1,358 ms
23,264 KB
testcase_13 AC 1,764 ms
22,056 KB
testcase_14 AC 1,405 ms
23,164 KB
testcase_15 AC 176 ms
20,644 KB
testcase_16 AC 1,453 ms
23,304 KB
testcase_17 TLE -
testcase_18 TLE -
testcase_19 TLE -
testcase_20 TLE -
testcase_21 TLE -
testcase_22 TLE -
testcase_23 TLE -
testcase_24 TLE -
testcase_25 TLE -
testcase_26 TLE -
testcase_27 AC 151 ms
18,972 KB
testcase_28 AC 152 ms
19,056 KB
testcase_29 AC 154 ms
18,972 KB
testcase_30 AC 153 ms
19,112 KB
testcase_31 AC 158 ms
19,028 KB
testcase_32 AC 71 ms
24,008 KB
testcase_33 AC 70 ms
23,940 KB
testcase_34 AC 69 ms
23,884 KB
testcase_35 AC 72 ms
24,064 KB
testcase_36 AC 70 ms
24,028 KB
testcase_37 AC 55 ms
18,172 KB
testcase_38 AC 149 ms
18,048 KB
testcase_39 TLE -
testcase_40 AC 10 ms
38,848 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int INF = 1e9 + 10;
const ll INFL = 4e18;

template <typename Cap>
struct max_flow {
    max_flow(int n) {
        this->n = n;
        G = vector<vector<tuple<int, int, Cap>>>(n);
    }
    void add_edge(int u, int v, Cap c) {
        G[u].push_back(make_tuple(v, G[v].size(), c));
        G[v].push_back(make_tuple(u, (int)G[u].size() - 1, 0));
    }
    Cap get_max_flow(int start, int goal) {
        Cap ret = 0;
        while (true) {
            vector<int> dst = calculate_distance(start);
            if (dst[goal] < 0) {
                return ret;
            }
            vector<int> removed(n, 0);
            while (true) {
                Cap add = flowing(start, goal, numeric_limits<Cap>::max(), removed, dst);
                if (add == 0) {
                    break;
                }
                ret += add;
            }
        }
        return ret;
    }

private:
    int n;
    vector<vector<tuple<int, int, Cap>>> G;
    vector<int> calculate_distance(int start) {
        vector<int> dst(n, -1);
        dst[start] = 0;
        queue<int> que;
        que.push(start);
        while (!que.empty()) {
            int now = que.front();
            que.pop();
            for (auto [nxt, _, cap] : G[now]) {
                if (cap > 0 && dst[nxt] == -1) {
                    dst[nxt] = dst[now] + 1;
                    que.push(nxt);
                }
            }
        }
        return dst;
    }
    Cap flowing(int now, int goal, Cap limit, vector<int> &removed, vector<int> &dst) {
        if (now == goal) {
            return limit;
        }
        while (removed[now] < (int)G[now].size()) {
            auto [nxt, rev, cap] = G[now][removed[now]];
            if (cap > 0 && dst[now] < dst[nxt]) {
                Cap flow = flowing(nxt, goal, min(limit, cap), removed, dst);
                if (flow > 0) {
                    get<2>(G[now][removed[now]]) -= flow;
                    get<2>(G[nxt][rev]) += flow;
                    return flow;
                }
            }
            removed[now]++;
        }
        return 0;
    }
};

const int MX = 5e5 + 5;

int main() {
    int H, W;
    cin >> H >> W;
    vector<vector<int>> A(H, vector<int>(W));
    for (int i = 0; i < H; i++) {
        for (int j = 0; j < W; j++) {
            cin >> A[i][j];
        }
    }

    vector<vector<pair<int, int>>> P(MX);
    for (int i = 0; i < H; i++) {
        for (int j = 0; j < W; j++) {
            P[A[i][j]].push_back({i, j});
        }
    }

    int ans = 0;
    for (int i = 1; i < MX; i++) {
        if (!P[i].empty()) {
            max_flow<ll> mf(H + W + 2);
            for (auto [x, y] : P[i]) {
                mf.add_edge(x, H + y, 1);
            }
            for (int j = 0; j < H; j++) {
                mf.add_edge(H + W, j, 1);
            }
            for (int j = 0; j < W; j++) {
                mf.add_edge(H + j, H + W + 1, 1);
            }
            ans += mf.get_max_flow(H + W, H + W + 1);
        }
    }

    cout << ans << endl;
}
0