KM算法板子2

//HDU2426

//

#include <iostream>
#include<stdio.h>
#include<string.h>
#include<string>
using namespace std;
const int maxn=550;
int visx[maxn],visy[maxn];
int link[maxn];
int lx[maxn],ly[maxn];
int w[maxn][maxn];
int slack[maxn];
int link1[maxn];
int n,m;
const int inf=999999999;
int find(int x)
{
    visx[x]=1;
    for(int i=0;i<m;i++)
    {
        if(w[x][i]==inf||visy[i])
            continue;
        int t=lx[x]+ly[i]-w[x][i];
        if(t==0)
        {
            visy[i]=1;
            if(link[i]==-1||find(link[i]))
            {
                link[i]=x;
                link1[x]=i;
                return 1;
            }
        }
        else
        {
            if(slack[i]>t)
                slack[i]=t;
        }
    }
    return 0;
}
int sk()
{
    memset(link,-1,sizeof(link));
    memset(link1,-1,sizeof(link1));
    memset(ly,0,sizeof(ly));
    for(int i=0;i<n;i++)
    {
        lx[i]=-inf;
        for(int j=0;j<m;j++)
        {
            if(w[i][j]!=inf)
                lx[i]=max(lx[i],w[i][j]);
        }
    }
    for(int i=0;i<n;i++)
    {
        //cout<<"m:"<<i<<endl;
        for(int j=0;j<m;j++)
            slack[j]=inf;
        while(true)
        {
            memset(visx,0,sizeof(visx));
            memset(visy,0,sizeof(visy));
            if(find(i))
                break;
           // cout<<"sb"<<endl;
            int d=inf;
           
            for(int j=0;j<m;j++)
            {
                if(!visy[j]&&d>slack[j])
                    d=slack[j];
            }
             if(d==inf)
               return -1;
            for(int j=0;j<n;j++)
            {
                if(visx[j])
                    lx[j]-=d;
            }
            for(int j=0;j<m;j++)
            {
                if(visy[j])
                    ly[j]+=d;
                else
                {
                   slack[j]-=d;
                }
            }
        }
    }
    int ans=0;
    int flag=0;
    for(int i=0;i<n;i++)
    {
        if(link1[i]==-1)
        {
            flag=1;
            break;
        }
        else
        {
            ans+=w[i][link1[i]];
           // cout<<"ans"<<ans<<endl;
        }
    }
    if(flag)
     ans=-1;
    return ans;
}
int main()
{
    int e;
    int tot=1;
    while(scanf("%d%d%d",&n,&m,&e)!=EOF)
    {
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<m;j++)
                w[i][j]=inf;
        }
        for(int i=1;i<=e;i++)
        {
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            if(c>=0)
            {
                w[a][b]=c;
              //  w[a][b]=c;
            }
        }
        printf("Case %d: %d\n",tot++,sk());
    }
    return 0;
}


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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值