hdu1533 Going Home KM算法求最小费用

二分图是这样一种图：所有顶点可以分成两个集：X和Y，其中X和Y中的任意两个在同一个集中的点都不相连，而来自X集的顶点与来自Y集的顶点有连线。当这些连线被赋于一定的权重时，这样的二分图便是带权二分图。

KM算法则就是对带权二分图进行操作  使两个集合在最大匹配的情况下 权值的和最大

推荐讲解KM算法的博客：http://philoscience.iteye.com/blog/1754498  http://www.cnblogs.com/wenruo/p/5264235.html

对于hdu1533这道题 人是一个集合 房子是一个集合 每个人都要走一定的距离（也可以看作是花费）走到一间房子中  我们要求最小距离 只需要把距离取反 这样就完全是一个km算法的模型 最后得出的结果再取反即可 

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<string>
#include<stack>
#include<queue>
#include<cmath>
#include<stack>
#include<list>
#include<map>
#include<set>
typedef long long ll;
using namespace std;
#define M 310
#define inf 0x3f3f3f3f

int nx,ny;
int link[M],lx[M],ly[M],slack[M];    //lx,ly为顶标，  nx,ny分别为x点集y点集的个数
int visx[M],visy[M],w[M][M];
char ch[105][105];
int cnt=0;
int DFS(int x)
{
    visx[x] = 1;
    for (int y = 1;y <= ny;y ++)
    {
        if (visy[y])
            continue;
        int t = lx[x] + ly[y] - w[x][y];
        if (t == 0)       
        {
            visy[y] = 1;
            if (link[y] == -1||DFS(link[y]))
            {
                link[y] = x;
                return 1;
            }
        }
        else if (slack[y] > t)  //不在相等子图中slack 取最小的
            slack[y] = t;
    }
    return 0;
}
int KM()
{
    int i,j;
    memset (link,-1,sizeof(link));
    memset (ly,0,sizeof(ly));
    for (i = 1;i <= nx;i ++)            //lx初始化为与它关联边中最大的
        for (j = 1,lx[i] = -inf;j <= ny;j ++)
            if (w[i][j] > lx[i])
                lx[i] = w[i][j];
    
    for (int x = 1;x <= nx;x ++)
    {
        for (i = 1;i <= ny;i ++)
            slack[i] = inf;
        while (1)
        {
            memset (visx,0,sizeof(visx));
            memset (visy,0,sizeof(visy));
            if (DFS(x))     //若成功（找到了增广轨），则该点增广完成，进入下一个点的增广
                break;  //若失败（没有找到增广轨），则需要改变一些点的标号，使得图中可行边的数量增加。
            //方法为：将所有在增广轨中（就是在增广过程中遍历到）的X方点的标号全部减去一个常数d，
            //所有在增广轨中的Y方点的标号全部加上一个常数d
            int d = inf;
            for (i = 1;i <= ny;i ++)
                if (!visy[i]&&d > slack[i])
                    d = slack[i];
            for (i = 1;i <= nx;i ++)
                if (visx[i])
                    lx[i] -= d;
            for (i = 1;i <= ny;i ++)  //修改顶标后，要把所有不在交错树中的Y顶点的slack值都减去d
                if (visy[i])
                    ly[i] += d;
                else
                    slack[i] -= d;
        }
    }
    int res = 0;
    for (i = 1;i <= ny;i ++)
        if (link[i] > -1)
            res += w[link[i]][i];
    return res;
}
struct node
{
    int x,y;
}man[100*105],hou[100*105];
int step(node a,node b)
{
    return abs(a.x-b.x)+abs(a.y-b.y);
}
int main ()
{
    int i,j,n,m;
    while (scanf ("%d%d",&n,&m)!=EOF)
    {
        cnt=0;
        getchar();
        if(m==0&&n==0)break;
        int cnt1=0,cnt2=0;
        
        for(i=0;i<n;i++)
        {
            for(j=0;j<m;j++)
            {
                scanf("%c",&ch[i][j]);
                if(ch[i][j]=='m')
                {
                    cnt1++;
                    man[cnt1].x=i;
                    man[cnt1].y=j;
                    
                }
                else if(ch[i][j]=='H')
                {
                    cnt2++;
                    hou[cnt2].x=i;
                    hou[cnt2].y=j;
                  
                }
            }
            getchar();
        }
        cnt=cnt1;
        nx=ny=cnt;
        for(i=1;i<=cnt;i++)
        {
            for(j=1;j<=cnt;j++)
            {
                w[i][j]=-step(man[i],hou[j]);//取反
            }
        }
        printf("%d\n",-KM());
        
    }
    return 0;
}