HDU3691 Nubulsa Expo【全局最小割】

最新推荐文章于 2020-06-17 20:09:20 发布

Master.Yi

最新推荐文章于 2020-06-17 20:09:20 发布

阅读量172

点赞数 1

分类专栏：网络流

本文链接：https://blog.csdn.net/C20181220_xiang_m_y/article/details/103324375

版权

网络流专栏收录该内容

35 篇文章 0 订阅

订阅专栏

题目描述：

n个点,m条边的无向图，给定源点，请选择汇点使得最小割最小，输出最小割的值。
n<=300,m<=50000,边权<=5000000.

题目分析：

源点一定在全局最小割的某一侧，只要汇点选在另一侧就可以了，所以问题就是求全局最小割。

全局最小割StoerWagner算法详解

论文：《浅谈无向图最小割问题的一些算法及应用》：
在这里插入图片描述

做n次最大邻接搜索，复杂度为 $O (n E l o g n)$ ，若不用堆优化为 $O(n^3)$

Code：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define maxn 305
using namespace std;
int n,m;
bool del[maxn],vis[maxn];
LL c[maxn][maxn],sum[maxn];
void merge(int s,int t){
	for(int i=1;i<=n;i++) c[s][i]+=c[t][i],c[i][s]+=c[i][t];
	del[t]=1;
}
LL StoerWagner(int pt){
	memset(sum,0,(n+1)<<3),memset(vis,0,(n+1));
	int s=0,t=0;
	for(int i=1;i<=n;i++) if(!del[i]) {t=i;break;}
	while(--pt){
		vis[s=t]=1,t=0;
		for(int i=1;i<=n;i++) if(!del[i]&&!vis[i]) sum[i]+=c[s][i];
		for(int i=1;i<=n;i++) if(!del[i]&&!vis[i]&&sum[i]>sum[t]) t=i; 
	}
	merge(s,t);
	return sum[t];
}
int main()
{
	while(scanf("%d%d%*d",&n,&m),n){
		memset(del,0,(n+1));
		for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) c[i][j]=0;
		for(int i=1,x,y,z;i<=m;i++) scanf("%d%d%d",&x,&y,&z),c[x][y]+=z,c[y][x]+=z;
		LL ans=1ll<<60;
		for(int i=n;i>1;i--) ans=min(ans,StoerWagner(i));
		printf("%lld\n",ans);
	}
}