kuangbin专题十 HDU3488 (二分图最优匹配+有向环覆盖)

题意：
有n个节点，m条有权单向路，要求用一个或者多个环覆盖住所有的节点。每个节点只能出现在一个环中，每个环中至少有两个节点。问最小边权花费为多少？
题解：
听说这道题可以用最优匹配去做，看了一下要求最小权值，那么肯定是边取负数进图，得出的值再取反就行了，这里还有个坑点就是，重边。。。MD。至于为什么可以用最优匹配来做有向环覆盖呢？你可以想一下，一个有向环图每个点是不是都有一个入度和出度呢？那么我们就可以把它拆成两个点i和i’，分别表示入度和出度，这样就构成了两份图了。不相信的话我们可以看一下例子，1->2,2->3,3->1,最终得出的是什么图呢？就是1与2’匹配，2与3’匹配，3与1’匹配，那么合起来之后不就是1->2->3->1了吗?是不是很神奇呢？我也觉得，然后我也理解了题目为啥会说一个或者多少环覆盖了。。原来就是想让你可以用这个算法去做这道题。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define INF 0x3f3f3f3f
const int MAXN=200+7;
int map[MAXN][MAXN];
int lx[MAXN],ly[MAXN];
int slack[MAXN];
bool visx[MAXN],visy[MAXN];
int cx[MAXN],cy[MAXN];
int nx,ny,n,m;
bool dfs(int u)
{
    visx[u]=true;
    for(int v=1;v<=ny;v++)
    {
        if(visy[v])
        continue;
        int tmp=lx[u]+ly[v]-map[u][v];
        if(tmp==0)
        {
            visy[v]=true;
            if(cy[v]==-1||dfs(cy[v]))
            {
                cx[u]=v;
                cy[v]=u;
                return true;
            }
        }
        else
        slack[v]=min(slack[v],tmp);
    }
    return false;
}
int KM()
{
    memset(cx,-1,sizeof(cx));
    memset(cy,-1,sizeof(cy));
    memset(ly,0,sizeof(ly));
    for(int i=1;i<=nx;i++)
    {
        lx[i]=-INF;
        for(int j=1;j<=ny;j++)
        if(map[i][j]>lx[i])
        lx[i]=map[i][j];
    }
    for(int i=1;i<=nx;i++)
    {
        for(int j=1;j<=ny;j++)
        slack[j]=INF;
        while(true)
        {
            memset(visx,false,sizeof(visx));
            memset(visy,false,sizeof(visy));
            if(dfs(i))
            break;
            int d=INF;
            for(int j=1;j<=ny;j++)
            if(!visy[j])
            d=min(d,slack[j]);
            for(int j=1;j<=nx;j++)
            if(visx[j])
            lx[j]-=d;
            for(int j=1;j<=ny;j++)
            {
                if(visy[j]) ly[j]+=d;
                else slack[j]-=d;
            }
        }
    }
    int res=0;
    for(int i=1;i<=ny;i++)
    if(cy[i]!=-1)
    res+=map[cy[i]][i];
    return -res;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        fill(map[0],map[0]+MAXN*MAXN,-INF);//二维数组，fill的赋值方法。 
    //  memset(map,-INF,sizeof(map));
        scanf("%d%d",&n,&m);
        for(int i=1;i<=m;i++)
        {
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            map[u][v]=max(map[u][v],-w);
        }
        ny=nx=n;
        int ans=KM();
        printf("%d\n",ans);
    }
}