hdu3488Tour KM算法

//给一个有向图,
//找出若干环,使得这些环覆盖所有点且每个点只能在一个环中
//问所得的所有环的所有边权值之和的最小值为多少
//对于每一个点只有一个入度和一个出度,那么将每个点拆成
//入度点和出度点,将所有入度点和所有出度点构成一个完备匹配
//由于是完备匹配,所以每个点的出度和入度都有一个不是自己的点
//相连,那么将完备匹配的所有点连接起来一定是若干环
//所得完备匹配的最小匹配即为答案
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std ;
const int maxn = 210 ;
const int inf = 0x3f3f3f3f ;
int lx[maxn] , ly[maxn] , slack[maxn] ;
int visx[maxn] , visy[maxn] , match[maxn] ;
int map[maxn][maxn] ;
int n  , m ;
int find(int x)
{
    visx[x] = 1 ;
    for(int i = 1;i <= n;i++)
    {
        if(visy[i])continue ;
        int tmp = lx[x] + ly[i] - map[x][i] ;
        if(tmp == 0)
        {
            visy[i] = 1 ;
            if(match[i] == -1 || find(match[i]))
            {
                match[i] = x ;
                return true ;
            }
        }
        else slack[i] = min(slack[i] , tmp) ;
    }
    return false ;
}
int KM()
{
    memset(ly , 0 , sizeof(ly)) ;
    memset(match , -1 , sizeof(match)) ;
    for(int i = 1;i <= n;i++)
    {
        lx[i] = -inf ;
        for(int j = 1;j <= n;j++)
        lx[i] = max(lx[i] , map[i][j]) ;
    }
    for(int i = 1;i <= n;i++)
    {
        for(int j = 1;j <= n;j++)
        slack[j] = inf ;
        while(1)
        {
            memset(visx , 0 ,sizeof(visx)) ;
            memset(visy , 0 ,sizeof(visy)) ;
            if(find(i))break ;
            int d = inf ;
            for(int j = 1;j <= n;j++)
            if(!visy[j])
            d = min(d , slack[j]) ;
            for(int j = 1;j <= n;j++)
            if(visx[j])
            lx[j] -= d ;
            for(int j = 1;j <= n;j++)
            if(visy[j])
            ly[j] += d ;
            else slack[j] -=d ;
        }
    }
    int ans = 0 ;
    for(int i = 1;i <= n;i++)
    ans += map[match[i]][i] ;
    return ans ;
}
int main()
{
    int t ;
    scanf("%d" , &t) ;
    while(t--)
    {
        scanf("%d%d" , &n , &m) ;
        for(int i = 1;i <= n;i++)
          for(int j = 1;j <= n;j++)
           map[i][j] = -inf ;
        while(m--)
        {
            int u , v , w ;
            scanf("%d%d%d" , &u , &v , &w) ;
            if(u != v && -w > map[u][v])
            map[u][v]  = -w ;
        }
        int ans = KM() ;
        printf("%d\n" , -ans) ;
    }
    return  0 ;
}

你可能感兴趣的:(hdu3488Tour KM算法)