結果

問題 No.119 旅行のツアーの問題
コンテスト
ユーザー codershifth
提出日時 2015-10-12 01:09:56
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 4,556 bytes
コンパイル時間 1,512 ms
コンパイル使用メモリ 172,188 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-21 10:54:18
合計ジャッジ時間 2,228 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

typedef long long ll;
typedef unsigned long long ull;

#define FOR(i,a,b) for(int (i)=(a);i<(b);i++)
#define REP(i,n) FOR(i,0,n)
#define RANGE(vec) (vec).begin(),(vec).end()


enum EdgeType {
    Directed,
    Undirected
};
template <typename T, EdgeType E=Directed> struct revCap {
    static T value(T cap) { (void)cap; return 0; }
};
template <typename T> struct revCap<T, Undirected> {
    static T value(T cap) { return cap; }
};
template <typename T, EdgeType S=Directed>
class MF {
public:
    struct Edge {
        Edge(int t, T c, int r) : to(t), cap(c), rev(r) {}
        int to;
        T cap;
        int rev;
    };
    std::vector<std::vector<Edge>> edges_;
    std::vector<int> iter_;
    std::vector<int> level_;

    MF(int v) : edges_(v) {}
    MF() {}
    ~MF() {}
    void init(int v) { edges_.resize(v); }
    void bfs(int from) {
        std::fill(level_.begin(), level_.end(), -1);
        std::queue<int> que;
        level_[from] = 0;
        que.push(from);
        while (!que.empty())
        {
            int v = que.front();
            que.pop();
            REP(i, edges_[v].size())
            {
                Edge &e = edges_[v][i];
                if (e.cap > 0 && level_[e.to] < 0)
                {
                    level_[e.to] = level_[v]+1;
                    que.push(e.to);
                }
            }
        }
        return;
    }
    T dfs(int v, int to, T f) {
        if (v == to)
            return f;
        for (int &i = iter_[v]; i < edges_[v].size(); ++i)
        {
            Edge &e = edges_[v][i];
            if (e.cap > 0 && level_[v] < level_[e.to])
            {
                int d = dfs(e.to, to, std::min(f, e.cap));
                if (d > 0)
                {
                    e.cap -= d;
                    edges_[e.to][e.rev].cap += d;
                    return d;
                }
            }
        }
        return 0;
    }
    T max_flow(int from, int to, T inf = std::numeric_limits<T>::max()) {
        int v = edges_.size();
        iter_.resize(v, 0);
        level_.resize(v, -1);
        int flow = 0;
        while (true)
        {
            bfs(from);
            if (level_[to] < 0)
                return flow;
            std::fill(iter_.begin(), iter_.end(), 0);
            int f;
            while ((f = dfs(from, to, inf)))
                flow += f;
        }
        return -1;
    }
    void add(int a, int b, T cap) {
        edges_[a].emplace_back(b,cap,edges_[b].size());
        edges_[b].emplace_back(a,revCap<T,S>::value(cap),edges_[a].size()-1);
    }
};

using namespace std;

class TravelTourProblem {
public:
    // 全ての国に訪れるケースを考えると各国に対してツアーに参加するかどうかの二値分類問題になる。
    // よってその場合は最小カットで解ける。(全流量からカットした分 B or C を取り除く。)
    //
    // 問題は国に訪れないケースも考える必要がある場合。
    // この場合は両方の総和 B+C を取り除くと考える。
    //
    void solve(void) {
            int N,M;
            cin>>N;

            int V = 2*N+2;
            int src = V-2;      // ツアー未参加側
            int sink = src+1;   // ツアー参加側
            int total = 0;
            MF<int> mf(V);
            const auto inf = (1<<30);
            REP(i,N)
            {
                int B,C;
                cin>>B>>C;
                total += (B+C);
                mf.add(src, i, C);
                mf.add(i+N, sink, B);
                mf.add(i,i+N,B+C); // 国に訪れないケースはここを切断すればよい。
            }
            cin>>M;
            REP(i,M)
            {
                int D,E;
                cin>>D>>E;
                //
                //     +-> E1->E2 -/-+
                //     |       v     |
                // src-+   +---+     +-> sink
                //     |   v         |
                //     +-> D1->D2 ---+
                //
                // E 番目のツアー参加への辺を切断しても D->sink が残っているなら
                // flow cut できないようにする。
                //
                mf.add(E+N,D,inf);
            }
            total -= mf.max_flow(src,sink,inf);
            cout<<total<<endl;
    }
};

#if 1
int main(int argc, char *argv[])
{
        ios::sync_with_stdio(false);
        auto obj = new TravelTourProblem();
        obj->solve();
        delete obj;
        return 0;
}
#endif
0