HDU 3549 Flow Problem 最大网络流

本题也是典型的最大流题解。

很好理解最大流建议大家还是去看 Introduction to algorithm这本经典书。

这里题解的特点是： 

1 使用静态邻接表建图

2 建立双向图，以便增加逆向流。

3 合并重复边优化图

使用邻接表解决网络流的关键是第二点。

一开始一直都看不明白为什么需要逆向增加流，很多博客也是一句话概括掉了，故此还是看书比较能弄懂。推荐看书，这里不重复了。

使用邻接矩阵也许方便点。不过这里使用邻接表的速度也不慢， 78 ms，而且感觉自己这样写特别一点的程序有助于更好理解。

毕竟只有理解了算法，才能随意修改实现细节，并保证结果正确。

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <queue>
#include <limits.h>
using std::queue;

const int MAX_EDGE = 1001;
const int MAX_VER = 16;
inline int min(int a, int b) { return a < b? a : b; }

struct Edge
{
	int src, des, wei, next;
};

Edge edge[MAX_EDGE];
int head[MAX_VER], N, M, E;

inline void addEdge(int src, int des, int wei)
{
	for (int e = head[src]; e; e = edge[e].next)
	{
		if (edge[e].des == des)
		{
			edge[e].wei += wei;
			return ;
		}
	}
	edge[E].src = src;
	edge[E].des = des;
	edge[E].wei = wei;
	edge[E].next = head[src];
	head[src] = E++;
}//静态邻接表矩阵构建

int parent[MAX_VER];
bool vis[MAX_VER];

bool bfs(int s, int t)
{
	memset(vis, 0, sizeof(vis));
	queue<int> qu;
	qu.push(s);
	vis[s] = true;

	while (!qu.empty())
	{
		int u = qu.front();
		qu.pop();

		for (int e = head[u]; e ; e = edge[e].next)
		{
			int v = edge[e].des;
			if (!vis[v] && edge[e].wei > 0)
			{
				qu.push(v);
				parent[v] = e;
				vis[v] = true;
			}
		}
	}
	return vis[t];
}

int FordFulkerson(int s, int t)
{
	int maxFlow = 0;
	while (bfs(s, t))
	{		
		int pathFlow = INT_MAX;
		for (int v = t; v != s; v = edge[parent[v]].src)
			pathFlow = min(pathFlow, edge[parent[v]].wei);

		for (int v = t; v != s; v = edge[parent[v]].src)
		{
			int e = parent[v];
			edge[e].wei -= pathFlow;
			addEdge(edge[e].des, edge[e].src, pathFlow);//增加逆向流
		}

		maxFlow += pathFlow;
	}
	return maxFlow;
}

int main()
{
	int T, src, des, wei;
	scanf("%d", &T);
	for (int t = 1; t <= T; t++)
	{
		memset(head, 0, sizeof(head));
		E = 1;
		scanf("%d %d", &N, &M);
		for (int i = 1; i <= M; i++)//edge 下标从1开始
		{
			scanf("%d %d %d", &src, &des, &wei);
			addEdge(src, des, wei);
			addEdge(des, src, 0);
		}
		printf("Case %d: %d\n", t, FordFulkerson(1, N));
	}
	return 0;
}