poj 2195 km算法初步

题意很明显，就是指在一个地图上面有n个人，有n个房间，然后要使得这些人能够到达一个房间，并且移动的步数最短。

看这道题目的人相信应该都是正在学习km算法的人。。

我推荐文库当中的这个资料，说实话，比较精炼，看起来不至于繁琐。。

转载：http://wenku.baidu.com/view/072015d133d4b14e852468a9.html

说说我对km算法的看法吧。。  之前一直想要搞清楚这个是什么样的原理，但是现在发现其实观察一下细节你就能发现这个算法的思想，特别是这个定标的设置，说实话我现在也不知道这个的原理是什么，但是基本上我的体会就是讲一个二分图当中的两堆顶点设为X、Y集合，因为是为了求其中的边权值最大的情况，所以只要选出X集合出发的最大的边权，那么记录为lx[i]，表示第i个顶点最大的边权，那么最优解绝对要小于等于lx[i]的所有之和。然后先找大的边权，然后再找次大，一直推下去，到最后一定能找到最优解。我对这个算法感觉就是这样。我明白我差不多只学会小用一下算法，里面的东西肯定还有很多思想要完善，所以自己只能再接再厉了。。。

代码如下：

#include<iostream>
#include<string.h>
#include<math.h>
#define Inf 0x7fffffff 
#define N 110
using namespace std;
int map[N][N],n,m,man_Count,home_Count;
int lx[N],ly[N],visx[N],visy[N],match[N],lack;
char temp[N];
struct node
{
    int x,y;       
}home[N],man[N];
bool find(int u)
{
     visx[u]=1;
     for(int i=0;i<home_Count;i++)
     {
         if(!visy[i])
         {
             int t=map[u][i]-lx[u]-ly[i];
             if(t==0)
             {
                 visy[i]=1;
                 if(match[i]==-1||find(match[i]))
                 {
                    match[i]=u;
                    return true;
                 }
             }
             else
             if(lack>t)lack=t;//这个是为找小一条边权做准备
         }
     }
     return false;
}
void km()
{
       memset(match,-1,sizeof(match));
       memset(ly,0,sizeof(ly));
       for(int i=0;i<man_Count;i++)
       {
          while(1)//只有找到最好的匹配的时候才能跳出循环
          {
             memset(visx,0,sizeof(visx));
             memset(visy,0,sizeof(visy));
             lack=Inf;
             if(find(i))break;
             for(int i=0;i<man_Count;i++)if(visx[i])lx[i]+=lack;
             for(int i=0;i<home_Count;i++)if(visy[i])ly[i]-=lack;
          }
       }
       return ;  
}
int main()
{
   while(scanf("%d%d",&n,&m))
   {
      if(n==0&&m==0)break;
      man_Count=home_Count=0;
      for(int i=0;i<n;i++)
      {
         scanf("%s",temp);
         for(int j=0;j<m;j++)
         {
            if(temp[j]=='H'){home[home_Count].x=i+1;home[home_Count++].y=j+1;}
            if(temp[j]=='m'){man[man_Count].x=i+1;man[man_Count++].y=j+1;}
         }
      }
      for(int i=0;i<man_Count;i++)
      {
         lx[i]=220;
         for(int j=0;j<home_Count;j++)
         {
            map[i][j]=abs(man[i].x-home[j].x)+abs(man[i].y-home[j].y);
            if(lx[i]>map[i][j])lx[i]=map[i][j];
         }
         for(int i=0;i<home_Count;i++)ly[i]=0;
      }
      int answer=0;
      km();
      for(int i=0;i<home_Count;i++)answer+=map[match[i]][i];
      printf("%d\n",answer);
   }
   return 0;    
}