hdu2255 (二分图最优匹配)

Kuhn Munkras算法

纯模板题。具体实现还不太懂，先贴模板。有空再研究。

很好的参考资料：http://www.cnblogs.com/skyming/archive/2012/02/18/2356919.html

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <map>
#include <math.h>
#include <algorithm>
using namespace std;
#define N 310
int mp[N][N];
bool visx[N],visy[N];
int lx[N],ly[N];
int match[N];
int n;

bool Hungary(int u){//匈牙利算法
    visx[u]=1;
    for(int i=0;i<n;i++){
        if(!visy[i] && lx[u]+ly[i]==mp[u][i]){
            visy[i]=1;
            if(match[i] == -1||Hungary(match[i])){
                match[i]=u;
                return true;
            }//if
        }//if
    }//for
    return false;
}

void KM_perfect_match(){
    int temp;
    memset(lx,0,sizeof(lx));
    memset(ly,0,sizeof(ly));
    for(int i=0;i<n;i++)//lx[i]为权值最大的边
        for(int j=0;j<n;j++)
            lx[i]=max(mp[i][j],lx[i]);
    for(int i=0;i<n;i++){//对n个点匹配
        while(1){
            memset(visx,0,sizeof(visx));
            memset(visy,0,sizeof(visy));
            if(Hungary(i)) break;//匹配成功
            else{//匹配失败，找最小值
                temp = INT_MAX;
                for(int j=0;j<n;j++) if(visx[j])//x在交错树中
                    for(int k=0;k<n;k++)
                    if(!visy[k]&&temp>lx[j]+ly[k]-mp[j][k])
                        temp=lx[j]+ly[k]-mp[j][k];
                for(int j=0;j<n;j++){//更新定标
                    if(visx[j]) lx[j]-=temp;
                    if(visy[j]) ly[j]+=temp;
                }
            }//else
        }//while
    }
}

int main(){
    int ans;
    while(scanf("%d",&n)!=EOF){
        ans=0;
        memset(match,-1,sizeof(match));
        for(int i=0;i<n;i++) for(int j=0;j<n;j++) scanf("%d",&mp[i][j]);

        KM_perfect_match();
        for(int i=0;i<n;i++) ans+=mp[match[i]][i];
        printf("%d\n",ans);
    }
    return 0;
}