tour 【二分图带权最小匹配】

//给一个有向图,
//找出若干环,使得这些环覆盖所有点且每个点只能在一个环中
//问所得的所有环的所有边权值之和的最小值为多少
//对于每一个点只有一个入度和一个出度,那么将每个点拆成
//入度点和出度点,将所有入度点和所有出度点构成一个完备匹配
//由于是完备匹配,所以每个点的出度和入度都有一个不是自己的点
//相连,那么将完备匹配的所有点连接起来一定是若干环
//所得完备匹配的最小匹配即为答案


//因为是最小匹配,所以初始化为负无穷,用负值计算;
//最后得出结果后加负号
#include<stdio.h>
#include<string.h>
#include<string>
#include<algorithm>
#include<iostream>
#include<math.h>
#define inf 0x3f3f3f3f
using namespace std;
const int N=205;

int e[N][N];
int la[N],lb[N];
int visa[N],visb[N];
int slack[N];
int match[N];
int n,m;

int dfs(int a)
{
    int i,j;
    visa[a]=1;
    for(i=1; i<=n; i++)
    {
        if(visb[i])
            continue;
        int temp=la[a]+lb[i]-e[a][i];
        if(temp==0)
        {
            visb[i]=1;
            if(match[i]==-1||dfs(match[i]))
            {
                match[i]=a;
                return true;
            }
        }
        else
            slack[i]=min(slack[i],temp);
    }
    return false;
}

int KM()
{
    int i,j;
    fill(la+1,la+n+1,-inf);
    for(i=1; i<=n; i++)
    {
        for(j=1; j<=n; j++)
            if(e[i][j]>la[i])
                la[i]=e[i][j];
    }
    memset(match,-1,sizeof(match));
    memset(lb,0,sizeof(lb));
    for(i=1; i<=n; i++)
    {
        fill(slack+1,slack+n+1,inf);
        while(1)
        {
            memset(visa,0,sizeof(visa));
            memset(visb,0,sizeof(visb));
            if(dfs(i))
                break;
            int d=inf;
            for(j=1; j<=n; j++)
                if(!visb[j])
                    d=min(slack[j],d);
            for(j=1; j<=n; j++)
            {
                if(visa[j]) la[j]-=d;
                if(visb[j]) lb[j]+=d;
                else slack[j]-=d;
            }
        }
    }
    int ans=0;
    for(i=1; i<=n; i++)
        if(match[i]!=-1)
             ans+=e[match[i]][i];
    return ans;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d %d",&n,&m);
        int i,j,a,b,c;
        for(i=1; i<=n; i++)
            for(j=1; j<=n; j++)
                e[i][j]=-inf;//要求最小值,所以初始化负无穷
        for(i=0; i<m; i++)
        {
            scanf("%d %d %d",&a,&b,&c);
            e[a][b]=max(e[a][b],-c);
        }
        printf("%d\n",-KM());//最后得出的值加-号
    }
    return 0;
}

 

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值