hdu 6081 度度熊的王国战略

#度度熊的王国战略

Problem Description

度度熊国王率领着喵哈哈族的勇士，准备进攻哗啦啦族。

哗啦啦族是一个强悍的民族，里面有充满智慧的谋士，拥有无穷力量的战士。

所以这一场战争，将会十分艰难。

为了更好的进攻哗啦啦族，度度熊决定首先应该从内部瓦解哗啦啦族。

第一步就是应该使得哗啦啦族内部不能同心齐力，需要内部有间隙。

哗啦啦族一共有n个将领，他们一共有m个强关系，摧毁每一个强关系都需要一定的代价。

现在度度熊命令你需要摧毁一些强关系，使得内部的将领，不能通过这些强关系，连成一个完整的连通块，以保证战争的顺利进行。

请问最少应该付出多少的代价。

Input

本题包含若干组测试数据。

第一行两个整数n，m，表示有n个将领，m个关系。

接下来m行，每行三个整数u,v,w。表示u将领和v将领之间存在一个强关系，摧毁这个强关系需要代价w

数据范围：

2<=n<=3000

1<=m<=100000

1<=u,v<=n

1<=w<=1000

Output

对于每组测试数据，输出最小需要的代价。

Sample Input

2 1
1 2 1
3 3
1 2 5
1 2 4
2 3 3

Sample Output

1
3


思路：
全局最小割

参考博客：无向图全局最小割

代码：

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
#include <queue>
using namespace std;
typedef long long ll;
const int MAXV = 3005;
const int MAXE = 200005;
const int INF = 0x3f3f3f3f;

struct edge
{
    int v,w,next;
} e[MAXE];
int head[MAXV],edgeNum;

void insert_edge(int u,int v,int w)
{
    e[edgeNum].v=v;
    e[edgeNum].w=w;
    e[edgeNum].next=head[u];
    head[u]=edgeNum++;
}


int fa[MAXV],link[MAXV];

int find(int x)
{
    return x==fa[x]?x:fa[x]=find(fa[x]);
}

void merge(int u,int v)
{
    int p=u;
    while(~link[p])
        p=link[p];
    link[p]=v;
    fa[v]=u;
}


int val[MAXV];
bool vis[MAXV];
int N,M;

int MCP(int cnt,int &s,int &t)
{
    memset(val,0,sizeof(val));
    memset(vis,0,sizeof(vis));
    priority_queue<pair<int,int> > que;
    t=1;
    while(--cnt)
    {
        s=t;
        vis[s]=true;
        for(int u=s; ~u; u=link[u])
        {
            for(int p=head[u]; ~p; p=e[p].next)
            {
                int v=find(e[p].v);
                if(!vis[v])
                    que.push(make_pair(val[v]+=e[p].w,v));
            }
        }
        t=0;
        while(!t)
        {
            if(que.empty())
                return 0;
            pair<int,int> tp=que.top();
            que.pop();
            if(val[tp.second]==tp.first)
                t=tp.second;
        }
    }
    return val[t];
}

int StoerWagner()
{
    int res=INF;
    for(int i=N,s,t; i>1; --i)
    {
        res=min(res,MCP(i,s,t));
        if(res==0)
            break;
        merge(s,t);
    }
    return res;
}

void init()
{
    memset(head,-1,sizeof(head));
    memset(link,-1,sizeof(link));
    for(int i=1; i<=N; ++i)
    {
        fa[i]=i;
    }
    edgeNum=0;
}

int main()
{
    while(~scanf("%d%d",&N,&M))
    {
        init();
        int u,v,w;
        for(int i=0; i<M; ++i)
        {
            scanf("%d%d%d",&u,&v,&w);
            insert_edge(u,v,w);
            insert_edge(v,u,w);
        }
        printf("%d\n",StoerWagner());
    }
    return 0;
}