hdu 3879(最小割模型求解最大权闭合图)

题意：

公司得到了一共N个可以作为通讯信号中转站的地址，而由于这些地址的地理位置差异，在不同的地方建造通讯中转站需要投入的成本也是不一样的，所幸在前期调查之后这些都是已知数据：建立第i个通讯中转站需要的成本为Pi（1≤i≤N）。

 

•另外公司调查得出了所有期望中的用户群，一共M个。关于第i个用户群的信息概括为Ai, Bi和Ci：这些用户会使用中转站Ai和中转站Bi进行通讯，公司可以获益Ci。（1≤i≤M, 1≤Ai, Bi≤N）

•THU集团的CS&T公司可以有选择的建立一些中转站（投入成本），为一些用户提供服务并获得收益（获益之和）。那么如何选择最终建立的中转站才能让公司的净获利最大呢？（净获利 = 获益之和 - 投入成本之和）

解题思路：这道题是《最小割模型在信息学竞赛中的应用》介绍到的“最大获利问题”，详细的证明过程要参看论文。这里只讲建图的思路。

把每个用户和每个站点看成是一个顶点，建立网络，从源点s向每个用户连一条容量为利润的边，每个用户向相关站点连一条容量为无穷大的边，每个站点向汇点连一条容量为成本的边。求出的最小割就是maxflow = （未被选的用户收益之和+被选的站点成本之和），设sum为总收益，我们要求的是（被选的用户收益之和-被选的站点成本之和），刚好等于sum-maxflow。至于原因参看论文。

这里要使用非递归版本的dinic，我的是递归版本的超时了。

TLE:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;

const int maxn = 60000;
const int inf = 0x3f3f3f3f;
struct Edge
{
	int to,next,flow;
}edge[maxn<<2];
int n,m,cnt,pre[maxn],layer[maxn];

void addedge(int u,int v,int flow)
{
	edge[cnt].to = v;
	edge[cnt].flow = flow;
	edge[cnt].next = pre[u];
	pre[u] = cnt++;
	swap(u,v);
	edge[cnt].to = v;
	edge[cnt].flow = 0;
	edge[cnt].next = pre[u];
	pre[u] = cnt++;
}

bool bfs(int s,int t)
{
	queue<int> q;
	memset(layer,0,sizeof(layer));
	layer[s] = 0;
	q.push(s);
	while(!q.empty())
	{
		int u = q.front();
		q.pop();
		if(u == t) return true;
		for(int i = pre[u]; i != -1; i = edge[i].next)
		{
			int v = edge[i].to;
			if(edge[i].flow > 0 && layer[v] == 0)
			{
				layer[v] = layer[u] + 1;
				q.push(v);
			}
		}
	}
	return false;
}

int dfs(int u,int t,int maxflow)
{
	if(u == t) return maxflow;
	int uflow = 0;
	for(int i = pre[u]; i != -1; i = edge[i].next)
	{
		int v = edge[i].to;
		if(layer[v] == layer[u] + 1 && edge[i].flow > 0)
		{
			int flow = min(maxflow - uflow,edge[i].flow);
			flow = dfs(v,t,flow);
			edge[i].flow -= flow;
			edge[i^1].flow += flow;
			uflow += flow;
			if(uflow == maxflow) break;
		}
	}
	if(uflow == 0)
		layer[u] = 0;
	return uflow;
}

int dinic(int s,int t)
{
	int maxflow = 0;
	while(bfs(s,t) == true)
		maxflow += dfs(s,t,inf);
	return maxflow;
}

int main()
{
	int s,t,u,v,w;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		s = 0, t = n + m + 1;
		cnt = 0;
		memset(pre,-1,sizeof(pre));
		for(int i = 1; i <= n; i++)
		{
			scanf("%d",&w);
			addedge(i,t,w);
		}
		int sum = 0;
		for(int i = 1; i <= m; i++)
		{
			scanf("%d%d%d",&u,&v,&w);
			sum += w;
			addedge(s,n+i,w);
			addedge(n+i,u,inf);
			addedge(n+i,v,inf);
		}
		printf("%d\n",sum - dinic(s,t));
	}
	return 0;
}