hdu 2255 km入门题

#include<stdio.h>  
#include<string.h>  
#include<stdlib.h>  
#define N 305  
#define inf 0x3fffffff  
int n,map[N][N],lx[N],ly[N],match[N];  
bool s[N],t[N];//s[]、t[]为左/右第i个点是否在匈牙利树中
int min(int a,int b)  
{  
    if(a>b)return b;  
    return a; 
}  
int max(int a,int b)  
{  
    if(a>b)return a;  
    return b;  
}  
bool find(int u)//匈牙利,匹配(找增广路)  
{  
    int i;  
     s[u]=true;  
    for(i=1;i<=n;i++)  
    {  
        if(!t[i]&&lx[u]+ly[i]==map[u][i])  
        {  
            t[i]=true;  
            if(match[i]==-1||find(match[i]))  
            {  
                match[i]=u;//u跟i匹配 
                return true;  
            }  
        }  
    }  
    return false;  
}  
int KM()  
{  
    int ans=0;  
    memset(match,-1,sizeof(match));  
    memset(lx,0,sizeof(lx));  
    memset(ly,0,sizeof(ly));  
    for(int i=1;i<=n;i++)  
        for(int j=1;j<=n;j++)  
           lx[i]=max(lx[i],map[i][j]);//初始化S顶标为最大权   
        for(i=1;i<=n;i++)//匹配每一个点    
        {  
            while(1)  
            {  
                memset(s,false,sizeof(s));  
                memset(t,false,sizeof(t));  
                if(find(i))break;//匹配成功  
                else//匹配失败,找最小值  
                {  
                    int a=inf;  
                    for(int j=1;j<=n;j++)  
                        if(s[j])//j在匈牙利树中 
                        {  
                            for(int k=1;k<=n;k++)  
                                if(!t[k])//k在匈牙利树外  
                                    a=min(a,lx[j]+ly[k]-map[j][k]);//顶标修改量,使得权值减小最小  
                        }  
                        for(j=1;j<=n;j++)//修改顶标  
                        {  
                            if(s[j])lx[j]-=a;//保证至少有一条边可以加入 
                            if(t[j])ly[j]+=a;//保证原来匹配的边修改后依然可以匹配  
                        }  
                }  
            }  
        }  
        for(i=1;i<=n;i++)//权值相加   
            ans+=map[match[i]][i];  
        return ans;  
}                 
int main()  
{  
    int i,j;  
    while(scanf("%d",&n)!=-1)  
    {  
        for(i=1;i<=n;i++)  
            for(j=1;j<=n;j++)  
                scanf("%d",&map[i][j]);  
            printf("%d\n",KM());  
    }  
    return 0;  
}  

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值