結果

問題 No.904 サメトロ
ユーザー MisterMister
提出日時 2020-04-21 23:24:27
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 35 ms / 1,000 ms
コード長 3,118 bytes
コンパイル時間 891 ms
コンパイル使用メモリ 85,180 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-10-09 14:16:24
合計ジャッジ時間 2,187 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 3 ms
6,820 KB
testcase_03 AC 4 ms
6,816 KB
testcase_04 AC 3 ms
6,816 KB
testcase_05 AC 2 ms
6,816 KB
testcase_06 AC 5 ms
6,816 KB
testcase_07 AC 3 ms
6,820 KB
testcase_08 AC 14 ms
6,820 KB
testcase_09 AC 2 ms
6,816 KB
testcase_10 AC 21 ms
6,816 KB
testcase_11 AC 19 ms
6,824 KB
testcase_12 AC 3 ms
6,816 KB
testcase_13 AC 2 ms
6,816 KB
testcase_14 AC 2 ms
6,820 KB
testcase_15 AC 6 ms
6,816 KB
testcase_16 AC 35 ms
6,820 KB
testcase_17 AC 2 ms
6,816 KB
testcase_18 AC 2 ms
6,820 KB
testcase_19 AC 2 ms
6,820 KB
testcase_20 AC 2 ms
6,816 KB
testcase_21 AC 4 ms
6,820 KB
testcase_22 AC 9 ms
6,816 KB
testcase_23 AC 14 ms
6,816 KB
testcase_24 AC 2 ms
6,820 KB
testcase_25 AC 2 ms
6,820 KB
testcase_26 AC 29 ms
6,816 KB
testcase_27 AC 14 ms
6,816 KB
testcase_28 AC 2 ms
6,816 KB
testcase_29 AC 3 ms
6,816 KB
testcase_30 AC 2 ms
6,820 KB
testcase_31 AC 2 ms
6,820 KB
testcase_32 AC 2 ms
6,820 KB
testcase_33 AC 2 ms
6,816 KB
testcase_34 AC 2 ms
6,816 KB
testcase_35 AC 2 ms
6,816 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <limits>

template <class Cap, class Cost>
struct MinCostFlow {
    struct Edge {
        int src, dst;
        Cap cap;
        Cost cost;
        Edge(int src, int dst, Cap cap, Cost cost)
            : src(src), dst(dst), cap(cap), cost(cost){};
    };

    using Edges = std::vector<Edge>;
    using Graph = std::vector<std::vector<int>>;

    Edges edges;
    Graph graph;
    std::vector<Cost> dist;
    std::vector<int> rev;

    const Cost INF = std::numeric_limits<Cost>::max() / 2;

    explicit MinCostFlow(int n) : graph(n), dist(n), rev(n) {}

    void span(int u, int v, Cap cap, Cost cost) {
        graph[u].push_back(edges.size());
        edges.emplace_back(u, v, cap, cost);

        graph[v].push_back(edges.size());
        edges.emplace_back(v, u, 0, -cost);
    }

    void bellman_ford(int s) {
        std::fill(dist.begin(), dist.end(), INF);
        dist[s] = 0;

        for (int i = 0; i < (int)graph.size(); ++i) {
            for (int eidx = 0; eidx < (int)edges.size(); ++eidx) {
                const auto& edge = edges[eidx];
                int u = edge.src, v = edge.dst;

                if (edge.cap > 0 &&
                    dist[u] < INF &&
                    dist[v] > dist[u] + edge.cost) {
                    dist[v] = dist[u] + edge.cost;
                    rev[v] = eidx;
                }
            }
        }
    }

    Cost exec(int s, int g, Cap flow) {
        Cost ret = 0;

        while (flow > 0) {
            bellman_ford(s);
            if (dist[g] == INF) break;

            Cap f = flow;
            int v = g;
            while (v != s) {
                const auto& edge = edges[rev[v]];
                f = std::min(f, edge.cap);
                v = edge.src;
            }

            flow -= f;
            ret += f * dist[g];

            v = g;
            while (v != s) {
                auto& edge = edges[rev[v]];
                auto& redge = edges[rev[v] ^ 1];
                edge.cap -= f;
                redge.cap += f;
                v = edge.src;
            }
        }
        return (flow > 0 ? -1 : ret);
    }
};

using lint = long long;
constexpr lint INF = 1LL << 30;

void solve() {
    int n;
    std::cin >> n;

    lint amax = 0, absum = 0;
    int s = n * 2 + 0, g = n * 2 + 1;
    MinCostFlow<lint, lint> mcf(n * 2 + 2);

    for (int v = 0; v < n; ++v) {
        if (v == 0) {
            mcf.span(s, v, INF, 1);
            mcf.span(v + n, g, INF, 0);
        } else {
            lint a, b;
            std::cin >> a >> b;

            mcf.span(s, v, a, -INF);
            mcf.span(v + n, g, b, -INF);
            amax += b;
            absum += a + b;
        }
    }

    for (int u = 0; u < n; ++u) {
        for (int v = 0; v < n; ++v) {
            if (u != v) mcf.span(u, v + n, INF, 0);
        }
    }

    mcf.span(s, g, INF, 0);

    lint c = mcf.exec(s, g, INF);
    lint amin = c + INF * absum;

    std::cout << amax - amin + 1 << std::endl;
}

int main() {
    std::cin.tie(nullptr);
    std::ios::sync_with_stdio(false);

    solve();

    return 0;
}
0