HDU-4289-コントロール(最大流量、スプリットポイント)

リンク:

https://vjudge.net/problem/HDU-4289

質問の意味:

あなたは、セキュリティの部長は、最近、テロリストのグループが、別の1(宛先)に1つの都市(ソース)からいくつかのWMD 1を輸送することを計画しているという極秘情報を受け取りました。あなたはそれらの日付、送信元と送信先を知っており、彼らは高速道路網を使用しています。
  高速道路網は、2つの異なる都市を結ぶ、双方向の高速道路で構成されています。車両は唯一の都市で高速道路網を終了/のみ入力することができます。
  テロリストが(それはSAは、この街である、である)観察下に都市を入力したとき、彼らはすぐに捕まえられるようにするには、いくつかの選択された都市で、いくつかのSA(特別エージェント)を見つけることがあります。
  すべての都市でSAを検索することも可能ですが、SAで街を制御することは、あなたの部門、都市から異なる場合があります一定の金額を、コスト、およびかもしれないので、あなたの予算は、すべての制御の全費用を負担することができない場合があります:都市は、あなたがいることを、都市のセットを識別しなければならない
  *テロリストのすべてのトラフィックは、セットの少なくとも1つの都市を渡す必要があります。
  *セット内のすべての都市を制御するコストの合計が最小になります。
  あなたは彼らの目的地へのテロリストのソースから取得することは常に可能であると仮定してもよいです。
-------------------------------------------------- ----------
大量破壊兵器1

アイデア:

オリジナルタイトルの最小値を求めると、sおよびtは今二点それぞれ、一つの入口側に接続された出口、各2つの点間も無向エッジに分割することができる。うカットポイント
になりました最小カット、最大流量に応じて、最大流量を実行することができます。

コード:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
//#include <memory.h>
#include <queue>
#include <set>
#include <map>
#include <algorithm>
#include <math.h>
#include <stack>
#include <string>
#define MINF 0x3f3f3f3f
using namespace std;
typedef long long LL;

const int MAXN = 200+10;
const int INF = 1e9;

struct Edge
{
    int from, to, cap;
};
vector<Edge> edges;
vector<int> G[MAXN*4];
int Dis[MAXN*4];
int n, m, s, t;

void AddEdge(int from, int to, int cap)
{
    edges.push_back(Edge{from, to, cap});
    edges.push_back(Edge{to, from, 0});
    G[from].push_back(edges.size()-2);
    G[to].push_back(edges.size()-1);
}

bool Bfs()
{
    memset(Dis, -1, sizeof(Dis));
    queue<int> que;
    que.push(s);
    Dis[s] = 0;
    while (!que.empty())
    {
        int u = que.front();
        que.pop();
//        cout << u << endl;
        for (int i = 0;i < G[u].size();i++)
        {
            Edge &e = edges[G[u][i]];
            if (e.cap > 0 && Dis[e.to] == -1)
            {
                Dis[e.to] = Dis[u]+1;
                que.push(e.to);
            }
        }
    }
    return Dis[t] != -1;
}

int Dfs(int u, int flow)
{
    if (u == t)
        return flow;
    int res = 0;
    for (int i = 0;i < G[u].size();i++)
    {
        Edge &e = edges[G[u][i]];
        if (e.cap > 0 && Dis[u]+1 == Dis[e.to])
        {
            int tmp = Dfs(e.to, min(flow, e.cap));
//            cout << "flow:" << e.from << ' ' << e.to << ' ' << tmp << endl;
            e.cap -= tmp;
            flow -= tmp;
            edges[G[u][i]^1].cap += tmp;
            res += tmp;
            if (flow == 0)
                break;
        }
    }
    if (res == 0)
        Dis[u] = -1;
    return res;
}

int MaxFlow()
{
    int res = 0;
    while (Bfs())
    {
        res += Dfs(s, INF);
//        cout << res << endl;
    }
    return res;
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    while (cin >> n >> m)
    {
        cin >> s >> t;
        for (int i = 0;i <= n*2;i++)
            G[i].clear();
        edges.clear();
        s = s*2-1;
        t = t*2;
        int w;
        for (int i = 1;i <= n;i++)
        {
            cin >> w;
            AddEdge(i*2-1, i*2, w);
        }
        int u, v;
        for (int i = 1;i <= m;i++)
        {
            cin >> u >> v;
            AddEdge(u*2, v*2-1, INF);
            AddEdge(v*2, u*2-1, INF);
        }
        LL res = MaxFlow();
        cout << res << endl;
    }

    return 0;
}

おすすめ

転載: www.cnblogs.com/YDDDD/p/11331155.html