最朴素的做法 n次最小割 O(mn^2)
“ 关于全局最小割的问题:
首先刚才提到的集合S和T是客观存在的与最优解相关但不知道的量。S+T=全顶点集
那么一定有顶点1属于S或1属于T。
刚才又提到了对任意点s属于S,t属于T,对st求最小割即可得到全局最优解。
当点1属于S时,我们只要枚举另一个点i,强行认为它属于T,求一遍割。那么当i真的属于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;
}