HDU.4700.Flow(构造 最小割树)

题目链接


(Description)

给定(n)以及(n)个点任意两点之间的最大流,求一张无向图满足给定条件。
(nleq100)

(Solution)

有些类似最小割树
我们可以构造一棵树,只要让树上的边成为割边,非树边容量为(0)就可以了。
每次找到当前点集中流量最小的边,设其流量为(c),然后根据(c)将点集分成两个集合,满足两个集合之间的点对的最大流是(c),集合内部的点的最大流(>c)。对于集合内部继续递归做即可。
划分集合的时候也是可以先随便找一个点(x)划分到左集合,将(flow[x][i]>c)的点(i)分到左集合,其余的点分到右集合,再判断一下左右集合是否满足之间的最大流(=c)即可。注意右集合为空时也无解((xdfrac{>c}{}udfrac{ c }{}vdfrac{>c}{}x),这样显然不行)。
复杂度(O(n^3))

注意(A[i])别写成(i)。。


//93MS	1316K
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
const int N=105;

int A[N],tmp[2][N],f[N][N],W[N][N];

inline int read()
{
	int now=0,f=1;register char c=gc();
	for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
	for(;isdigit(c);now=now*10+c-48,c=gc());
	return now*f;
}
bool DFS(int l,int r)
{
	if(l==r) return 1;
	int mn=1<<30,cnt[2]={0,1}; tmp[1][1]=A[l];
	for(int i=l; i<=r; ++i)
		for(int j=i+1; j<=r; ++j) mn=std::min(mn,f[A[i]][A[j]]);
	for(int s=A[l],i=l+1,t; i<=r; ++i) t=f[s][A[i]]>mn, tmp[t][++cnt[t]]=A[i];
	if(!cnt[0]) return 0;
	for(int i=1; i<=cnt[0]; ++i)
		for(int j=1; j<=cnt[1]; ++j)
			if(f[tmp[0][i]][tmp[1][j]]!=mn) return 0;
	for(int i=l,t=1; t<=cnt[0]; ++i) A[i]=tmp[0][t++];
	for(int i=l+cnt[0],t=1; i<=r; ++i) A[i]=tmp[1][t++];
	W[A[l]][A[r]]=mn, W[A[r]][A[l]]=mn;
	return DFS(l,l+cnt[0]-1)&&DFS(l+cnt[0],r);
}

int main()
{
	int n;
	while(~scanf("%d",&n))
	{
		memset(W,0,sizeof W);
		for(int i=1; i<=n; ++i)
			for(int j=1; j<=n; ++j) f[i][j]=read();
		for(int i=1; i<=n; ++i) A[i]=i;
		if(!DFS(1,n)) puts("NO");
		else
		{
			puts("YES");
			for(int i=1; i<=n; ++i)
				for(int j=(W[i][i]=-1,1); j<=n; ++j) printf("%d%c",W[i][j]," 
"[j==n]);
		}
	}

	return 0;
}
原文地址:https://www.cnblogs.com/SovietPower/p/10533310.html