2017CCPC湘潭邀请赛B.Roads（全局最小割+最小生成树）

2017CCPC湘潭邀请赛B.Roads（全局最小割+最小生成树）

题目大意：

给你一个n点m边的图，现在要在图中删除一些边，使得最小生成树的值严格地增大（不连通则最小生成树的值视为INF），问你最少需要删除多少条边？

思路：

这个题让我了解了什么是全局最小割：

无向图的全局最小割即为将图G(V)中的点划分为两个部分S,T使得任意一条从S到T的边的权值之和最小。

也可以理解为将这个无向图分为两个不相链接的部分所需要的最小割的花费。

 

此题中，假设新图的点集为S，边集为E，开始时S = {},E = {}把权值第一小的所有候选边和点加入S和E，求出所有连通点集的最小割取最小值，再接着依次把第i小的候选边和点加入S和E，求最小割取最小，也就是模拟一个最小生成树的过程，不过每次加入一些权值相同的边集合，我们就得用全局最小割SW算法求一下当前的割的代价是多少，然后每次取个最小值就好了。

全局最小割sw算法模板

//stoer-wagner算法模板开始
int edge[MAXN][MAXN],dist[MAXN];
bool vis[MAXN],bin[MAXN];
void sw_init()
{
    memset(edge,0,sizeof(edge));
    memset(bin,0,sizeof(bin));
}
int merge(int &s,int &t)
{
    memset(dist,0,sizeof(dist));
    memset(vis,0,sizeof(vis));
    int k,mincut,maxc;
    for(int i=1;i<=n;i++)
    {
        k=-1, maxc=-1;
        for(int j=1;j<=n;j++) if(!bin[j] && !vis[j] && dist[j] > maxc) {k=j;maxc=dist[j];}
        if(k == -1) return mincut;
        vis[k]=true;
        s=t; t=k;
        mincut=maxc;
        for(int j=1;j<=n;j++) if(!bin[j] && !vis[j]) dist[j]+=edge[k][j];
    }
    return mincut;
}
int stoer_wagner()
{
    int mincut=INF,s,t,ans;
    for(int i=1;i<=n-1;i++)
    {
        ans=merge(s,t);
        bin[t]=1;
        if(ans!=0 && ans<mincut) mincut=ans;
        for(int j=1;j<=n;j++) if(!bin[j]) edge[s][j]=(edge[j][s]+=edge[j][t]);
    }
    return mincut;
}
//stoer-wagner算法模板结束

此题代码：

#include<cstdio>
#include<cstring>
#include<map>
#include<vector>
#include<algorithm>
#define MAXN 53
#define MAXM MAXN*MAXN
#define INF 0x3f3f3f3f
using namespace std;
​
int n,m;
struct Node{
    int u,v,w;
};
//并查集模板开始
int par[MAXN],ran[MAXN];
void ufs_init(){for(int i=1;i<=n;i++) par[i]=i,ran[i]=0;}
int find(int x)
{
    if(par[x] == x) return x;
    else return( par[x] = find(par[x]) );
}
void unite(int x,int y)
{
    x=find(x),y=find(y);
    if(x == y) return;
    if(ran[x] < ran[y]) par[x]=y;
    else
    {
        par[y]=x;
        if(ran[x] == ran[y]) ran[x]++;
    }
}
bool isSame(int x,int y){return( find(x) == find(y) );}
//并查集模板结束
//stoer-wagner算法模板开始
int edge[MAXN][MAXN],dist[MAXN];
bool vis[MAXN],bin[MAXN];
void sw_init()
{
    memset(edge,0,sizeof(edge));
    memset(bin,0,sizeof(bin));
}
int merge(int &s,int &t)
{
    memset(dist,0,sizeof(dist));
    memset(vis,0,sizeof(vis));
    int k,mincut,maxc;
    for(int i=1;i<=n;i++)
    {
        k=-1, maxc=-1;
        for(int j=1;j<=n;j++) if(!bin[j] && !vis[j] && dist[j] > maxc) {k=j;maxc=dist[j];}
        if(k == -1) return mincut;
        vis[k]=true;
        s=t; t=k;
        mincut=maxc;
        for(int j=1;j<=n;j++) if(!bin[j] && !vis[j]) dist[j]+=edge[k][j];
    }
    return mincut;
}
int stoer_wagner()
{
    int mincut=INF,s,t,ans;
    for(int i=1;i<=n-1;i++)
    {
        ans=merge(s,t);
        bin[t]=1;
        if(ans!=0 && ans<mincut) mincut=ans;
        for(int j=1;j<=n;j++) if(!bin[j]) edge[s][j]=(edge[j][s]+=edge[j][t]);
    }
    return mincut;
}
//stoer-wagner算法模板结束
bool cmp(Node a,Node b){return a.w<b.w;}
int main()
{
    Node E[MAXM];
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        for(int i=1;i<=m;i++)
        {
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            E[i].u=a, E[i].v=b, E[i].w=c;
        }
        sort(E+1,E+m+1,cmp);
        int ans=INF;
        int edge_cnt=0;
        ufs_init();
        //模拟最小生成树的建立过程，每次加入最小的相同权值的边集合
        for(int i=1;i<=m;)
        {
            int now_i=i,now_w=E[i].w;
            sw_init();
            for(;now_w==E[i].w && i<=m;i++)//建图，加入最小的相同权值的边集合
            {
                if(!isSame(E[i].u,E[i].v))
                {
                    edge[find(E[i].u)][find(E[i].v)]++;
                    edge[find(E[i].v)][find(E[i].u)]++;
                }
            }
            //加入完最小边集合后求全局最小割，答案取最小值
            ans=min(ans,stoer_wagner());
​
            for(int i=now_i;now_w==E[i].w && i<=m;i++)//逐步建立最小生成树
            {
                if(!isSame(E[i].u,E[i].v))
                {
                    unite(E[i].u,E[i].v);
                    edge_cnt++;//最多选n-1条
                }
            }
            if(edge_cnt>=n-1) break;
        }
        printf("%d\n",ans);
    }
}
​