POJ 1459 Power Network （多源点多汇点求最大流问题）

题目链接

POJ-1459 Power Network

题目分析

 题意：

简单来说，由多个源点向多个汇点传输数据，每个源点传出的数据是有限的，每个汇点接受的数据也是有限的，某些点之间可以传输一定量的数据，求最大流。

思路：

相比于普通的最大流问题，这个地方的源点和汇点不唯一，而且源点的最大流出容量也有所限制，其实我们只要简单处理了一下，就可以将这类问题转化为普通的最大流问题。

我们在图中新增一个超级源点s和一个超级汇点t，由超级源点s向每个源点连一条容量为对应最大流出容量的边，再从每个汇点向超级汇点t连一条容量为对应最大流入容量的边，这样处理之后，问题就转化为了求超级源点到超级汇点之间的最大流了

代码区

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#include<string>
#include<fstream>
#include<vector>
#include<stack>
#include <map>
#include <iomanip>
#define bug cout << "**********" << endl
#define show(x,y) cout<<"["<<x<<","<<y<<"] "
//#define LOCAL = 1;
using namespace std;
typedef long long ll;
const int inf = 0x3f3f3f3f;
const ll mod = 1e6 + 3;
const int Max = 1e4 + 10;

struct Edge
{
	int to, next;
	int flow;
}edge[Max << 2];

int n, out, in, m;
int head[110], tot;
int dis[110];			//表示顶点到源点的距离标号

void init()
{
	memset(head, -1, sizeof(head));tot = 0;
}

void add(int u, int v, int flow)
{
	edge[tot].to = v;
	edge[tot].flow = flow;
	edge[tot].next = head[u];
	head[u] = tot++;
}

bool bfs()												//判断残留网络是否连通，并求出dis
{
	queue<int>q;
	memset(dis, -1, sizeof(dis));
	dis[0] = 0;
	q.push(0);
	while (!q.empty())
	{
		int u = q.front();q.pop();
		for (int i = head[u]; i != -1; i = edge[i].next)
		{
			int v = edge[i].to;
			if (dis[v] == -1 && edge[i].flow > 0)		//可以借助边i到达新的结点
			{
				dis[v] = dis[u] + 1;					//求顶点到源点的距离编号
				q.push(v);
			}
		}
	}
	return dis[n] != -1;								//确认是否连通
}

int dfs(int u, int flow_in)								//求增广路
{
	if (u == n) return flow_in;
	int flow_out = 0;									//记录这一点实际流出的流量
	for (int i = head[u]; i != -1;i = edge[i].next)
	{
		int v = edge[i].to;
		if(dis[v] == dis[u] + 1 && edge[i].flow > 0 )
		{
			int flow_part = dfs(v, min(flow_in, edge[i].flow));

			if (flow_part == 0) continue;

			flow_in -= flow_part;						//流出了一部分，剩余可分配流入就减少了
			flow_out += flow_part;						//记录这一点最大的流出

			edge[i].flow -= flow_part;
			edge[i ^ 1].flow += flow_part;				//常规操作了...
			if (flow_in == 0)
				break;
		}
	}
	return flow_out;
}


int max_flow()
{
	int sum = 0;
	while (bfs())
	{
		sum += dfs(0, inf);
	}
	return sum;
}

int main()
{
#ifdef LOCAL
	freopen("input.txt", "r", stdin);
	freopen("output.txt", "w", stdout);
#endif
	while (scanf("%d%d%d%d", &n, &out, &in, &m) != EOF)
	{
		init();
		n++;										//新增超级源点0和超超级汇点n
		for (int i = 1, u, v, flow;i <= m;i++)
		{
			scanf(" (%d,%d)%d", &u, &v, &flow);
			u++, v++;								//把0当作超级源点
			add(u, v, flow);add(v, u, 0);
		}

		for (int i = 1, v, flow;i <= out;i++)		//由超级源点向每个源点连一条容量为对应最大流出容量的边
		{
			scanf(" (%d)%d", &v, &flow);v++;
			add(0, v, flow);add(v, 0, 0);
		}

		for (int i = 1, u, flow;i <= in; i++)		//从每个汇点向超级汇点连一条容量为对应最大流入容量的边
		{
			scanf(" (%d)%d", &u, &flow);u++;
			add(u, n, flow);add(n, u, 0);
		}
		printf("%d\n", max_flow());
	}
	return 0;
}