全局最小割 学习总结

全局最小割的意思是在一个无向图中任取S和T,求最小割的最小值

还有一种描述是删掉无向图中的边使得其不连通的最小代价

当然，这种题目可以用分治+最小割来求解

但是时间复杂度大约在O(n^4)左右

有一种更好的求解方法可以在O(n^3)的时间复杂度内求解

做法是这样的：

首先对于图中任意两点S->T

要么S和T不在一个集合里时是答案，答案显然是S和T的最小割

否则S和T在一个集合里，我们可以将S和T缩成一个点，不难证明这样是等效的

我们模拟这个过程，每次任取S和T跑最小割，时间复杂度大概跟分治+最小割差不多OwO

但是我们注意到这一步的S和T都是任取的，也就是我们只需要求出一组解就可以了

这样就有了一种构造算法：

1、一开始A集合为空，我们任取一点加入A集合中

2、定义W(A,x)表示A集合所有点到x的边权和

3、每次寻找W(A,x)最大的点(相同的话任取),加入A集合中并更新其他的W值

4、最后加入的两个点则分别是S和T，其最小割为最后加入时W(A,T)

这个过程直接模拟显然是O(n^2)的，加上缩点最多进行O(n)次

所以总时间复杂度O(n^3)，可以用堆优化，但是稠密图跑得很慢很慢

 

BZOJ 3345

 
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm> 
using namespace std;
 
const int maxn=510;
const int oo=0x7fffffff;
int n,m,u,v,d;
int S,T,cut;
int f[maxn][maxn];
int w[maxn];
bool com[maxn],vis[maxn];
 
void Get_ans(){
    cut=oo;S=T=-1;
    memset(vis,false,sizeof(vis));
    memset(w,0,sizeof(w));
    for(int i=1;i<=n;++i){
        int mx=-oo,tmp;
        for(int j=1;j<=n;++j){
            if(!com[j]&&!vis[j]&&w[j]>mx){
                mx=w[j];tmp=j;
            }
        }
        if(T==tmp)return;
        S=T;T=tmp;
        cut=mx;vis[tmp]=true;
        for(int j=1;j<=n;++j)if(!com[j]&&!vis[j])w[j]+=f[tmp][j];
    }return;
}
int Get_SW(){
    memset(com,false,sizeof(com));
    int ans=oo;
    for(int i=1;i<n;++i){
        Get_ans();
        ans=min(ans,cut);
        if(ans==0)return 0;
        com[T]=true;
        for(int j=1;j<=n;++j){
            if(!com[j]){f[S][j]+=f[T][j];f[j][S]+=f[j][T];}
        }
    }return ans;
}
 
int main(){
    scanf("%d%d",&n,&m);
    memset(f,0,sizeof(f));
    for(int i=1;i<=m;++i){
        scanf("%d%d%d",&u,&v,&d);
        f[u][v]+=d;f[v][u]+=d;
    }
    printf("%d\n",Get_SW());
    return 0;
}

hdu 3691

注意到这里的S是给定的

但是不难发现这并没有什么卵用，答案还是全局最小割

因为设最后全局最小割中一定有一个T跟S不属于一个集合OwO

又因为是全局最小割，所以这组S->T的最小割就是全局最小割

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;

typedef long long LL;
const int maxn=310;
const LL oo=1LL<<60;
int n,m,k;
int u,v,w;
int S,T;
LL ans,cut;
LL f[maxn][maxn];
LL W[maxn];
bool com[maxn],vis[maxn];

void read(int &num){
    num=0;char ch=getchar();
    while(ch<'!')ch=getchar();
    while(ch>='0'&&ch<='9')num=num*10+ch-'0',ch=getchar();
}
void Get_ans(){
    cut=oo;S=T=-1;
    memset(vis,false,sizeof(vis));
    memset(W,0,sizeof(W));
    for(int i=1;i<=n;++i){
        LL mx=-oo;int tmp;
        for(int j=1;j<=n;++j){
            if(!com[j]&&!vis[j]&&W[j]>mx){mx=W[j];tmp=j;}
        }
        if(T==tmp)return;
        S=T;T=tmp;
        cut=mx;vis[tmp]=true;
        for(int j=1;j<=n;++j){
            if(!com[j]&&!vis[j]){W[j]+=f[tmp][j];}
        }
    }return;
}
LL Get_SW(){
    ans=oo;memset(com,0,sizeof(com));
    for(int i=1;i<n;++i){
        Get_ans();
        ans=min(ans,cut);
        if(ans==0)return ans;
        com[T]=true;
        for(int j=1;j<=n;++j){
            if(!com[j]){f[S][j]+=f[T][j];f[j][S]+=f[j][T];}
        }
    }return ans;
}

int main(){
    while(scanf("%d%d%d",&n,&m,&k)==3){
        if(!n&&!m&&!k)break;
        memset(f,0,sizeof(f));
        for(int i=1;i<=m;++i){
            read(u);read(v);read(w);
            f[u][v]+=w;f[v][u]+=w;
        }
        printf("%d\n",Get_SW());
    }return 0;
}