【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=500;
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]=1;
    for (int i=1;i<=numm;i++)
    {
        int wt=fx[u]+fy[i]-g[u][i];
        if (!visy[i] && wt==0)
        {
            visy[i]=1;
            if (lk[i]==-1 || Hungary_dfs(lk[i]))
            {
                lk[i]=u;
                return 1;
            }
        }
        else if (wt<slack[i]) slack[i]=wt; 
    }
    return 0;
}

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

void init()
{
    numm=numh=0;
    for (int i=0;i<n;i++)
    {
        char c;
        getchar();
        for (int j=0;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=1;i<=numm;i++)
        for (int j=1;j<=numh;j++)
        {
            g[i][j]=-1*(abs(mx[i]-hx[j])+abs(my[i]-hy[j]));
        }
}

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