Control HDU-4289 （最小割）

题目连接

题意：

N个城市， M条路 

一群恐怖份子想要从一个城市 S 将 一批武器送到 E 城市，

你现在需要在城市设防，不同的城市需要的花费不同，

然后是 M行，每一行是 一个城市到另一个城市有路

思路：

1）读题画图，题意就是把找出权值和最小的点集，把这些点消掉，使 S  和 T 不连通

2）怎么消点呢？看了别人的思路还

3）把每个城市都拆成两个点，然后两点连线，权值为这个城市的设防值。

4）然后城市之间的路设为INF而且是双向的

5）这样问题就变成了求最小割的问题

6）最小割等于最大流

AC代码：

#include<iostream>
#include<cstring>
#include<math.h>
#include<algorithm>
#include<queue>
#include<stdio.h>
using namespace std;
const int maxn = 1010;
const int maxm = 1e5 + 10;
const int inf = 0x3f3f3f3f;
int M, N;
int s, e;
struct ac{
    int v, c, nex;
}edge[maxm];

int head[maxn], dis[maxn], curedge[maxn], cnt;
void init() {
    cnt = 0;
    memset(head, -1, sizeof(head));
}
void addedge(int u, int v, int c) {
	// 正向建边
    edge[cnt] = {v, c, head[u]};
    head[u] = cnt++;
    // 反向建边, 流量为0
    edge[cnt] = {u, 0, head[v]};
    head[v] = cnt++;
}
bool bfs() {
    queue<int> que;
    que.push(s);
    memset(dis, 0, sizeof(dis)); // 对图进行分层
    dis[s] = 1;
    while (!que.empty()) {
        int u = que.front();
        que.pop();
        for (int i = head[u]; i != -1; i = edge[i].nex) {
        	int v = edge[i].v;
        	int c = edge[i].c;
        	// 如果节点v已经分过层或者u->v流量为0, continue
            if (dis[v] || c == 0) continue;
            dis[v] = dis[u] + 1; // 对v进行标记并加入队列
            que.push(v);
        }
    }
    return dis[e] > 0;  // 判断是否存在增广路，s是否能到达e
}
 
int dfs(int u, int flow) { // 增广路走到u点的最小流量为flow
    if (u == e || flow == 0) return flow;
    // 遍历u的所有出边
    for (int &i = curedge[u]; i != -1; i = edge[i].nex) { // 当前弧优化
    	int v = edge[i].v;
    	int c = edge[i].c;
    	// 判断能否u->v增广
        if (dis[v] != dis[u] + 1 || c == 0) continue;
        int d = dfs(v, min(flow, c));
        if (d > 0) { // 找到一条增广路，修改增广路上的正反向边
            edge[i].c -= d;
            edge[i^1].c += d;
            return d;
        }            
    }
    dis[u] = -1; // // 炸点优化
    return 0;
}
 
int Dinic() {
    int sum = 0, d;
    while (bfs()) { // 判读是否存在增广路
        for (int i = 0; i <= N * 2; ++i) curedge[i] = head[i]; // copy head数组，在dfs中可以直接得到下一条没有被增广过的边
        while ((d = dfs(s, inf)) > 0) sum += d; // 多次dfs找增广路
	}
    return sum;
}
 
int main () {
	int Cost;
	int u, v;
	while (scanf("%d%d", &N, &M) != EOF) {
		init();
		scanf("%d%d", &s, &e);
		e = e + N;
		for (int i = 1; i <= N; ++i) {
			scanf("%d", &Cost);
			addedge(i, i + N, Cost);
		}
		for (int i = 1; i <= M; ++i) {
			scanf("%d%d", &u, &v);
			addedge(u + N, v, inf);
			addedge(v + N, u, inf);
		}
		printf("%d\n", Dinic());
	}
}