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

时间:2023-03-10 04:46:00
HDU.4700.Flow(构造 最小割树)

题目链接


\(Description\)

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

\(n\leq100\)。

\(Solution\)

有些类似最小割树

我们可以构造一棵树,只要让树上的边成为割边,非树边容量为\(0\)就可以了。

每次找到当前点集中流量最小的边,设其流量为\(c\),然后根据\(c\)将点集分成两个集合,满足两个集合之间的点对的最大流是\(c\),集合内部的点的最大流\(>c\)。对于集合内部继续递归做即可。

划分集合的时候也是可以先随便找一个点\(x\)划分到左集合,将\(flow[x][i]>c\)的点\(i\)分到左集合,其余的点分到右集合,再判断一下左右集合是否满足之间的最大流\(=c\)即可。注意右集合为空时也无解(\(x\dfrac{>c}{}u\dfrac{\ c\ }{}v\dfrac{>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]," \n"[j==n]);
}
} return 0;
}