[全局最小割] BZOJ 3345

最朴素的做法 n次最小割 O(mn^2)

“ 关于全局最小割的问题：

首先刚才提到的集合S和T是客观存在的与最优解相关但不知道的量。S+T=全顶点集
那么一定有顶点1属于S或1属于T。
刚才又提到了对任意点s属于S，t属于T，对st求最小割即可得到全局最优解。
当点1属于S时，我们只要枚举另一个点i，强行认为它属于T，求一遍割。那么当i真的属于T的时候我们就得到最优解了。否则得到了一个没什么用的解。 ”

——hzqtrue 

Orrz

Stoer-Wagner

http://blog.sina.com.cn/s/blog_700906660100v7vb.html

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define cl(x) memset(x,0,sizeof(x))
#define oo 1<<30
using namespace std;

int n,m,ret=oo;
int M[505][505];
int vst[505],dis[505],v[505];

inline void Stoer_Wagner(int n)
{
	for (int i=1;i<=n;i++)
		v[i]=i;
	while (n>1)
	{
		int k,pre=1;
		cl(dis); cl(vst); vst[1]=1;
		for (int i=1;i<n;i++)
		{
			k=0;
			for (int j=1;j<=n;j++)
				if (!vst[v[j]])
				{
					dis[v[j]]+=M[v[pre]][v[j]];
					if (dis[v[k]]<dis[v[j]])
						k=j;
				}
			vst[v[k]]=1;
			if (i==n-1)
			{
				ret=min(ret,dis[v[k]]);
				for (int j=1;j<=n;j++)
				{
					M[v[j]][v[pre]]+=M[v[j]][v[k]];
					M[v[pre]][v[j]]+=M[v[j]][v[k]];
				}
				v[k]=v[n--];
			}
			pre=k;
		}
	}
}

int main()
{
	int _u,_v,_w;
	freopen("t.in","r",stdin);
	freopen("t.out","w",stdout);
	while (~scanf("%d%d",&n,&m))
	{
		for (int i=1;i<=m;i++)
		{
			scanf("%d%d%d",&_u,&_v,&_w);
			M[_u][_v]+=_w;
			M[_v][_u]+=_w;
		}
		ret=oo;
		Stoer_Wagner(n);
		printf("%d\n",ret);
	}
	return 0;
}