【KM】POJ2195/HDU1533-Going home

//最近没什么时间quq据说长得帅的人都在切八中，然而长得丑的人只能水水裸题

【题目大意】

给出一张地图及人和房屋的位置，求出每个人回到不同房屋所具有的最小代价和。

【思路】

最小权匹配，先O(n^2)求出每个人到任意房屋的距离后取反，就变成了裸的KM。数据范围很小，O（n^4）可以0ms过。

【错误点】

1.slack在每次循环中药重新设置为无穷大。2.内循环和外循环的循环变量重名了。3.变量名取太乱了用错了一次。

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cmath>
 #include<algorithm>
 using namespace std;
 const int INF=0x7fffffff;
 const int MAXN=;
 int n,m;
 int mx[MAXN],my[MAXN],hx[MAXN],hy[MAXN],numm,numh;
 int g[MAXN][MAXN];
 int visx[MAXN],visy[MAXN],fx[MAXN],fy[MAXN],slack[MAXN],lk[MAXN];

 int Hungary_dfs(int u)
 {
     visx[u]=;
     for (int i=;i<=numm;i++)
     {
         int wt=fx[u]+fy[i]-g[u][i];
         if (!visy[i] && wt==)
         {
             visy[i]=;
             if (lk[i]==- || Hungary_dfs(lk[i]))
             {
                 lk[i]=u;
                 return ;
             }
         }
         else if (wt<slack[i]) slack[i]=wt;
     }
     return ;
 }

 int KM()
 {
     memset(lk,-,sizeof(lk));
     for (int i=;i<=numm;i++)
     {
         fx[i]=-INF;
         fy[i]=;
         for (int j=;j<=numm;j++) fx[i]=max(fx[i],g[i][j]);
     }

     for (int i=;i<=numm;i++)
     {
         memset(visx,,sizeof(visx));
         memset(visy,,sizeof(visy));
         memset(slack,,sizeof(slack));//slack每次都要初始化为无穷大
         while (!Hungary_dfs(i))
         {
             int delta=INF;
             for (int j=;j<=numm;j++) if (!visy[j]) delta=min(delta,slack[j]);
             for (int j=;j<=numm;j++)//这里忘掉了外面还套了一重循环，重复用了i！
             {
                 if (visx[j])
                 {
                     visx[j]=;
                     fx[j]-=delta;
                 }
                 if (visy[j])
                 {
                     visy[j]=;
                     fy[j]+=delta;
                 }
             }
         }
     }

     int ret=;
     for (int i=;i<=numm;i++) ret+=fx[i]+fy[i];//这里numm写成了n
     return ret;
 }

 void init()
 {
     numm=numh=;
     for (int i=;i<n;i++)
     {
         char c;
         getchar();
         for (int j=;j<m;j++)
         {
             scanf("%c",&c);
             if (c=='m')
             {
                 mx[++numm]=i;
                 my[numm]=j;
             }
             else if (c=='H')
             {
                 hx[++numh]=i;
                 hy[numh]=j;
             }
         }
     }
     for (int i=;i<=numm;i++)
         for (int j=;j<=numh;j++)
         {
             g[i][j]=-*(abs(mx[i]-hx[j])+abs(my[i]-hy[j]));
         }
 }

 int main()
 {
     while (scanf("%d%d",&n,&m))
     {
         if (n==m && m==) break;
         init();
         cout<<KM()*(-)<<endl;
     }
     return ;
 }