关闭

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

标签: 网络流
143人阅读 评论(0) 收藏 举报
分类:

题意:

公司得到了一共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;
}


0
0

猜你在找
【直播】计算机视觉原理及实战—屈教授
【套餐】深度学习入门视频课程—唐宇迪
【套餐】Hadoop生态系统零基础入门--侯勇蛟
【套餐】嵌入式Linux C编程基础--朱有鹏
【套餐】2017软考系统集成项目——任铄
【套餐】Android 5.x顶级视频课程——李宁
【直播】广义线性模型及其应用——李科
【直播】从0到1 区块链的概念到实践
【直播】机器学习之凸优化——马博士
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:65370次
    • 积分:4499
    • 等级:
    • 排名:第6543名
    • 原创:395篇
    • 转载:54篇
    • 译文:0篇
    • 评论:3条
    最新评论