hdu3488Tour(KM)

题目请戳这里

题目大意：给n个城市，m条有向边，求一些有向环，保证每个点只出现在一个环中并且保证总权值最小。

题目分析：每个点在一个环中，所以对于每个城市都有且仅有一条入边和出边，将一个点拆成2个点，建图跑一下KM即可。因为求最小值，所以所有边权值取反，求最大权完美匹配。因为题目保证有解，所以不需要考虑无解。

trick：有重边。

详情请见代码：

#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N = 205;
const int M = 30005;
const int inf = 0x3f3f3f3f;
int n,m;
int dis[N][N];
int cx[N],cy[N],match[N];
int lx[N],ly[N];
int slack;
bool path(int u)
{
    cx[u] = 1;
    for(int v = 1;v <= n;v ++)
    {
        if(!cy[v])
        {
            int t = lx[u] + ly[v] - dis[u][v];
            if(t)
            {
                if(slack > t)
                    slack = t;
            }
            else
            {
                cy[v] = 1;
                if(match[v] == -1 || path(match[v]))
                {
                    match[v] = u;
                    return true;
                }
            }
        }
    }
    return false;
}

int KM()
{
    int i,j,ret;
    for(i = 1;i <= n;i ++)
    {
        lx[i] = -inf;
        ly[i] = 0;
        for(j = 1;j <= n;j ++)
        {
            if(lx[i] < dis[i][j])
                lx[i] = dis[i][j];
        }
    }
    memset(match,-1,sizeof(match));
    for(int u = 1;u <= n;u ++)
    {
        while(1)
        {
            memset(cx,0,sizeof(cx));
            memset(cy,0,sizeof(cy));
            slack = inf;
            if(path(u))
                break;
            for(i = 1;i <= n;i ++)
            {
                if(cx[i])
                    lx[i] -= slack;
                if(cy[i])
                    ly[i] += slack;
            }
        }
    }
    ret = 0;
    for(i = 1;i <= n;i ++)
        ret += dis[match[i]][i];
    return -ret;
}

int main()
{
    int t;
    int i,j,a,b,c;
    scanf("%d",&t);
    while(t --)
    {
        scanf("%d%d",&n,&m);
        for(i = 1;i <= n;i ++)
            for(j = 1;j <= n;j ++)
                dis[i][j] = -inf;
        for(i = 1;i <= m;i ++)
        {
            scanf("%d%d%d",&a,&b,&c);
            if(dis[a][b] < -c)
                dis[a][b] = -c;
        }
        printf("%d\n",KM());
    }
    return 0;
}
//171MS	400K