nbu 2429 Transfer stations

最新推荐文章于 2020-12-24 07:08:43 发布
最新推荐文章于 2020-12-24 07:08:43 发布
阅读量568 收藏
点赞数
分类专栏： 网络流
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/gotoac/article/details/8293271
版权

题目链接:http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2429


题目大意:

现在有一项工程,有n(n<=5000,编号1~n)个站点可以建立,建立i站点需要花费pi.

而假设a站点和b站点均建立了,那么可以获利C<a,b>,这样的关系有m个(m<=50000).

问建立那些站点可获得最大纯利润.


题目思路:

裸的最大权闭合图.

这种题型的建图有两种方法.

一:(建完后的图的原图完全不同)

(1)将所有边看成点(成为利润点),并各向构成这条边的两个点建立一条容量为正无穷大的边.

(3)用源点(建立一个超级源点)向所有利润点建立一条容量为其获利的边.

(2)将原来所有点(成为花费点),向汇点(建立一个超级汇点)建立一条容量为其花费的边.

这样之后就是要求出该图的最小割,而且最小割必然是割在和源点或者汇点相连的边上,因为其余的边容量均为无穷大.

通俗的说最小割的割集就是我们不需要的建立的站点的花费和无法获利的边.

这种方法的正确性,其实可以这样考虑,从源点到汇点必然是这样的路径: 源点--->利润边--->花费边--->汇点,

而最小割肯定是割在这三条边上容量最小的一条.

将问题转为最小割之后,由于最小割容量等于最大流,所以用所有利润边的获利和减去求得最大流即可.


二:(只需要在原图的基础上改进)

对于第一种我们可以分析,

点数=n+m+2,n最大是5000,m最大是50000,所以构造出来的图会有55002个点.

边数=n+m+m,大概11万条边.

显然这是一个 很庞大的图,效率并不好.

所以我们需要一个更好的建模方式.

(1)同样需要一个超级源点S,一个超级汇点T.

(2)源点连向所有点,容量为U (U为所有获利的和),称为零边.

(3)所有点连向汇点,容量为U+2*pi-di (pi为i的点权di为所有所有和i相关联的边的权值和),称为差边.

ps:U的作用是为了使差边的容量非负.

这样我们可以看一下从源点到汇点的路径是什么样的.

假设原图中存在一条利润边a-->b

那么必然存在这样几条路径使得S和T连通.

S ---> a ---> T

S ---> b ---> T

S ---> a ---> b ---> T

S ---> b ---> a ---> T

为了使得S和T不连通,那么最小割可能为以下两种.

S ---> a 和 S ---> b ,即都割零边,割的容量为2U

a ---> T 和 b ---> T ,即都割差边,割的容量为2U+2*(pa+pb)-(va+vb)

S ---> a 和 b ---> a 和 b ---> T , 即割零边和差边,割的容量为2U+2*pb-vb+weight(a,b).

S ---> b 和 a ---> b 和 a ---> T , 即割零边和差边,割的容量为2U+2*pa-va+weight(a,b).

如果割零边,表示a站点和b站点都不建立,显然利润为0.

反之,a和b 都建立,利润应该是 weigt(a,b)-(pa+pb).

如果与a,与b关联的边均只有一条,即(a,b),那么va = vb = weigt(a,b).

那么割差边时,割的容量2U+2*(pa+pb)-(va+vb) =2U+2*(pa+pb)-2*weigt(a,b).

假设还有一割点c与a关联,且除a没有其他点和c关联.

1).假设c站点建立,那么多了一个割为c ---> T.

那么割的容量为 3U+2*(pa+pb+pc)-(va+vb+vc) = 3U+2*(pa+pb+pc)-2*(weigt(a,b)+weight(a,c)).

2).假设c站点不建立,那么多两个割为S ---> c 和 a ---> c

那么割的容量为 3U+2*(pa+pb) - (va+vb) + weight(a,c) 

  = 3U+2*(pa+pb)-2*weigt(a,b)-weight(a,c)+weight(a,c)

          = 3U+2*(pa+pb)-2*weigt(a,b)

同理可证当c或b还有其他关联边时也满足该式子.

又所以利润 = (nU-割[S,T])/2.



代码(一):

#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;

#define ll long long
#define ls rt<<1
#define rs ls|1
#define lson l,mid,ls
#define rson mid+1,r,rs
#define middle (l+r)>>1
#define eps (1e-8)
#define clr_all(x,c) memset(x,c,sizeof(x))
#define clr(x,c,n) memset(x,c,sizeof(x[0])*(n+1))
#define MOD 1000000009
#define INF 0x3f3f3f3f
#define PI acos(-1.0)
#define _max(x,y) (((x)>(y))? (x):(y))
#define _min(x,y) (((x)<(y))? (x):(y))
#define _abs(x) ((x)<0? (-(x)):(x))
#define getmin(x,y) (x= ((x)<0 || (y)<(x))? (y):(x))
#define getmax(x,y) (x= ((y)>(x))? (y):(x))
template <class T> void _swap(T &x,T &y){T t=x;x=y;y=t;}
int TS,cas=1;
const int M=60000+5;
int n,m;

struct sap{
	typedef int type;
	struct edge{
		int to,next;
		type flow;
	}edg[999999];
	int n,m;
	int g[M],cur[M],pre[M];
	int dis[M],gap[M];
	int q[M],head,tail;
	void init(int _n){n=_n,m=0,clr_all(g,-1);}
	void insert(int u,int v,type f,type c=0){
		edg[m].to=v,edg[m].flow=f,edg[m].next=g[u],g[u]=m++;
		edg[m].to=u,edg[m].flow=0,edg[m].next=g[v],g[v]=m++;
	}
	bool bfs(int s,int t){
		clr(gap,0,n),clr(dis,-1,n);
		gap[dis[t]=0]++,dis[s]=-1;
		for(q[head=tail=0]=t;head<=tail;){
			int u=q[head++];
			for(int i=g[u];i!=-1;i=edg[i].next){
				edge& e=edg[i],ee=edg[i^1];
				if(dis[e.to]==-1 && ee.flow>0){
					gap[dis[e.to]=dis[u]+1]++;
					q[++tail]=e.to;
				}
			}
		}
		return dis[s]!=-1;
	}
	type maxFlow(int s,int t){
		if(!bfs(s,t)) return 0;
		type res=0,a;
		int i;
		for(i=0;i<n;i++) cur[i]=g[i];
		pre[s]=s,cur[s]=g[s],cur[t]=g[t];
		for(int u=s;dis[s]<n;){
			if(u==t){
				for(a=-1;u!=s;u=pre[u])
					getmin(a,edg[cur[pre[u]]].flow);
				for(u=t;u!=s;u=pre[u]){
					edg[cur[pre[u]]].flow-=a;
					edg[cur[pre[u]]^1].flow+=a;
				}
				res+=a;
			}
			bool ok=0;
			for(i=cur[u];i!=-1;i=edg[i].next){
				edge& e=edg[i];
				if(dis[u]==dis[e.to]+1 && e.flow>0){
					cur[u]=i,pre[e.to]=u,u=e.to;
					ok=1;break;
				}
			}
			if(ok) continue;
			int mindis=n-1;
			for(i=g[u];i!=-1;i=edg[i].next){
				edge& e=edg[i];
				if(mindis>dis[e.to] && e.flow>0)
					mindis=dis[e.to],cur[u]=i;
			}
			if(--gap[dis[u]]==0) break;
			gap[dis[u]=mindis+1]++,u=pre[u];
		}
		return res;
	}
}p;

void run(){
	int i,j;
	p.init(n+m+2);
	int tot=0;
	for(i=1;i<=n;i++){
		scanf("%d",&j);
		p.insert(i,n+m+1,j);
	}
	for(i=1;i<=m;i++){
		int a,b,c;
		scanf("%d%d%d",&a,&b,&c);
		p.insert(0,n+i,c);
		p.insert(n+i,a,INF);
		p.insert(n+i,b,INF);
		tot+=c;
	}
	printf("%d\n",tot-p.maxFlow(0,n+m+1));
}

void preSof(){
}

int main(){
    //freopen("input.txt","r",stdin);
    //freopen("output.txt","w",stdout);
    preSof();
    //run();
    while(~scanf("%d%d",&n,&m)) run();
    //for(scanf("%d",&TS);cas<=TS;cas++) run();
    return 0;
}

代码(二): 

#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;

#define ll long long
#define ls rt<<1
#define rs ls|1
#define lson l,mid,ls
#define rson mid+1,r,rs
#define middle (l+r)>>1
#define eps (1e-8)
#define clr_all(x,c) memset(x,c,sizeof(x))
#define clr(x,c,n) memset(x,c,sizeof(x[0])*(n+1))
#define MOD 1000000009
#define INF 0x3f3f3f3f
#define PI acos(-1.0)
#define _max(x,y) (((x)>(y))? (x):(y))
#define _min(x,y) (((x)<(y))? (x):(y))
#define _abs(x) ((x)<0? (-(x)):(x))
#define getmin(x,y) (x= ((x)<0 || (y)<(x))? (y):(x))
#define getmax(x,y) (x= ((y)>(x))? (y):(x))
template <class T> void _swap(T &x,T &y){T t=x;x=y;y=t;}
int TS,cas=1;
const int M=5000+5;
int n,m,U;
int pr[M],d[M];
int a[M*10],b[M*10],c[M*10];

struct sap{
	typedef int type;
	struct edge{
		int to,next;
		type flow;
	}edg[999999];
	int n,m;
	int g[M],cur[M],pre[M];
	int dis[M],gap[M];
	int q[M],head,tail;
	void init(int _n){n=_n,m=0,clr_all(g,-1);}
	void insert(int u,int v,type uf,type vf=0){
		edg[m].to=v,edg[m].flow=uf,edg[m].next=g[u],g[u]=m++;
		edg[m].to=u,edg[m].flow=vf,edg[m].next=g[v],g[v]=m++;
	}
	bool bfs(int s,int t){
		clr(gap,0,n),clr(dis,-1,n);
		gap[dis[t]=0]++,dis[s]=-1;
		for(q[head=tail=0]=t;head<=tail;){
			int u=q[head++];
			for(int i=g[u];i!=-1;i=edg[i].next){
				edge& e=edg[i],ee=edg[i^1];
				if(dis[e.to]==-1 && ee.flow>0){
					gap[dis[e.to]=dis[u]+1]++;
					q[++tail]=e.to;
				}
			}
		}
		return dis[s]!=-1;
	}
	type maxFlow(int s,int t){
		if(!bfs(s,t)) return 0;
		type res=0,a;
		int i;
		for(i=0;i<n;i++) cur[i]=g[i];
		pre[s]=s,cur[s]=g[s],cur[t]=g[t];
		for(int u=s;dis[s]<n;){
			if(u==t){
				for(a=-1;u!=s;u=pre[u])
					getmin(a,edg[cur[pre[u]]].flow);
				for(u=t;u!=s;u=pre[u]){
					edg[cur[pre[u]]].flow-=a;
					edg[cur[pre[u]]^1].flow+=a;
				}
				res+=a;
			}
			bool ok=0;
			for(i=cur[u];i!=-1;i=edg[i].next){
				edge& e=edg[i];
				if(dis[u]==dis[e.to]+1 && e.flow>0){
					cur[u]=i,pre[e.to]=u,u=e.to;
					ok=1;break;
				}
			}
			if(ok) continue;
			int mindis=n-1;
			for(i=g[u];i!=-1;i=edg[i].next){
				edge& e=edg[i];
				if(mindis>dis[e.to] && e.flow>0)
					mindis=dis[e.to],cur[u]=i;
			}
			if(--gap[dis[u]]==0) break;
			gap[dis[u]=mindis+1]++,u=pre[u];
		}
		return res;
	}
}p;

void run(){
	int i,j;
	p.init(n+1);
	int tot=0;
	for(i=1;i<=n;i++)
		scanf("%d",&pr[i]);
	clr(d,0,n);
	for(U=0,i=1;i<=m;i++){
		scanf("%d%d%d",&a[i],&b[i],&c[i]);
		U+=c[i],d[a[i]]+=c[i],d[b[i]]+=c[i];
		p.insert(a[i],b[i],c[i],c[i]);
	}
	for(i=1;i<=n;i++){
		p.insert(0,i,U);
		p.insert(i,n+1,U+2*pr[i]-d[i]);
	}
	printf("%d\n",(U*n-p.maxFlow(0,n+1))/2);
}

void preSof(){
}

int main(){
    //freopen("input.txt","r",stdin);
    //freopen("output.txt","w",stdout);
    preSof();
    //run();
    while(~scanf("%d%d",&n,&m)) run();
    //for(scanf("%d",&TS);cas<=TS;cas++) run();
    return 0;
}


gotoac
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
nbu15_NBU15
weixin_29876887的博客
12-24 125
beigeTexturehomogeneous,short-fibredDensityat20°C,g/cm30.99Servicetemperaturerange*,°C,approx.–40to130Droppoint,DINISO2176,°C,approx.>220Workedpenetration,DINISO2137,(ASTMD217...
笔记：NBU客户端安装
Lizi_TT的博客
05-07 3325
NBU 7.6 安装 echo "kernel.sem=300 307200 32 1024">> /etc/sysctl.conf 解压后安装： [root@nbu76 NetBackup_7.6.0.1_LinuxR_x86_64]# ./install Symantec Installation Script Copyright 1993 - 2013 Symante...
nbu15_NBU常用的命令
weixin_33087827的博客
12-24 800
1.在命令提示符下输入available_media，查看当前磁带使用情况2.在命令提示符下输入bpmedialist命令查看当前磁带使用情况及是否冻结3.在命令提示符号下输入bpexpdate –m a00001 –d 0来使磁带过期4.在命令提示符号下输入bpmedia –unfreeze –m a00001来使磁带解除冻结状态5.在命令提示符号下输入vmquery –m a00001来查看磁...
nbu15_克鲁勃NBU 15
weixin_29088655的博客
12-24 472
ISOFLEXNBU15â用于滚动轴承的特殊润滑脂高速度和负载描述:包装尺寸:ISOFLEXNBU1550克PE罐400克墨盒1公斤罐ISOFLEXNBU15是一款高速具有良好压力的油脂ab-吸附能力。它由一个酯油/合成油的组合烃油/矿物油和钡复合肥皂。它提供良好的耐磨和防腐保护，并且耐水，耐环境媒体和氧化。•良好的压力吸收容量25公斤桶•高速•良好的磨损保护•良好的防腐蚀性•抗老化...
NBU10.3.0.1介质
最新发布
03-21
vxupdate_nb_10.3.0.1_windows_x64.sja vxupdate_nb_10.3.0.1_suse_x64.sja vxupdate_nb_10.3.0.1_redhat_x64.sja NetBackup_10.3.0.1_Win.zip NetBackup_10.3.0.1_LinuxS_x86_64.tar.gz NetBackup_10.3.0.1_LinuxR_...
NBU.rar_NBU
09-22
VERITA NetBackup（NBU）是一款广泛应用于企业级数据保护的备份和恢复软件。它以其高效、稳定和全面的数据管理能力而闻名。标题中的“NBU.rar_NBU”可能指的是NetBackup相关的压缩文件，其中包含了关于该软件的信息...
NBU8.3下载
08-11
veritas netbackup 8.3 安装包
最新NBU售后培训PPT
01-06
NBU 备份软件 售后培训PPT 全是硬货干活
nbu&mysql备份恢复.docx
04-24
NBU&MySQL备份恢复 NBU（NetBackup）是一款功能强大的备份恢复软件，而MySQL是一款流行的关系型数据库管理系统。本文将详细介绍如何使用NBU来备份和恢复MySQL数据库。 NBU备份mysql的优势 使用NBU备份mysql数据库...
nbu 2431 星际转移问题
Go.To.AC的专栏
12-13 763
题目链接：http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2431   题目大意： 中文题意，不解释。   题目思路： 源自王晓东《线性规划和网络流24题》。 原题的解法，一种枚举答案+动态建图，一种二分答案+多次建图，因为原题最多30天。 但是这题最多可以是1000多天，所以二分答案+多次建图会TLE。 不过两种解法的思路
nbu 2428 Risk
Go.To.AC的专栏
12-14 575
题目链接:http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2428 题目大意: 有n个区域(1 现在要通过一次调兵(设从a调兵,那么只能调到和a相邻的自己的区域)使得和敌人接壤的区域的最小士兵数最大,且自己的所有区域的兵的数量为正. 题目思路: 假设答案为x. (1)建立一个超级源点S,并把所有的敌人的区域看成一个
poj 1459 Power Network （测模板）
Go.To.AC的专栏
11-01 532
题目链接:http://poj.org/problem?id=1459 题目大意: 多源多汇网络流. 题目思路: 学习网络流的第一题,也是拿来测模板. 定一个虚源点指向所有源点,容量为每个原点的Pmax. 定一个虚汇点使所有汇点指向它,容量为每个汇点的Cmax. 流一下得到最大流即可. 该题用邻接表写速度可能会变慢,因为邻接表里,会出现重边,但是不影响最终结果.
nbu 2427 Pigs
Go.To.AC的专栏
12-14 532
题目链接:http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2427 题目大意: 有m(1)个猪圈,n(1 第i人第i天去买猪,有ki把钥匙,最多买bi头猪. 第i天的时候,钥匙对应的猪圈会打开,这个时候可以对猪圈任意分配猪的数量(和不会超过打开的猪圈的猪的数量之和.) 问最多能卖多少猪. 题目思路: (1)建立一
nbu 2430 war of energy
Go.To.AC的专栏
12-13 523
题目链接：http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2430 题目大意： 有n个城市（编号1~n），m条管道，每条管道连接两个城市，管道i最多可容纳c[i]的能量通过，求从1最多能向n传送多少能量。 0 题目思路： 裸的最大流，按照题意建图，1为源点，n为汇点，求最大流。 代码： #include
NK 2122 太空飞行计划问题
Go.To.AC的专栏
12-01 507
题目链接：http://acm.nankai.edu.cn/p2122.html 题目大意： 中文题意，不解释。 题目思路： 最大权闭合图。 详解：http://blog.csdn.net/gotoac/article/details/8247023 由于要输出非选择的实验和仪器，喜欢用sap的兄弟最后从源点搜一遍，能遍历到的点就是选择了的。 代码： #
nbu 2432 羊群之防
Go.To.AC的专栏
12-13 496
题目链接：http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2432   题目大意： 中文题意，不解释。   题目思路： 裸的费用流。 羊进洞所需的时间为费用。 建图： 每只羊对应每个洞的建立一条容量为1的边。 设一个源点，连向每只羊，容量为1。 设一个汇点，每个洞连向汇点，容量为1。 建完图后发现，还是一个二分图，所以
hdu 2063 过山车
Go.To.AC的专栏
11-30 424
题目连接：http://acm.hdu.edu.cn/showproblem.php?pid=2063 题目大意： 中文题，不解释。 题目思路： 裸的二分匹配。 匈牙利，网络流。。。随便搞 刚看完匈牙利，测测自己的代码。。。 代码： #include #include #include #include #include #include #in
SYMANTEC+NBU+7.5+安装配置说明(for+linux).docx
12-21
首先，用户需要访问Symantec官方NBU网站()，并使用已注册或新注册的账户下载适用于Linux的NBU 7.5安装介质，推荐下载ESD映像文件，格式为NetBackup_7.5_LinuxR_x86_64.tar.gz。 安装步骤开始于以root权限登录主...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值