HDU 1533

http://acm.hdu.edu.cn/showproblem.php?pid=1533

人和房子数量相同，每个人进房子，费用是人到房子的曼哈顿距离，求最小费用

可用最小费用最大流求解，建立虚拟的源点和汇点即可

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std ;
int ABS(int x)
{
    return x>?x:-x ;
}
struct point{
    int x,y ;
}mm[],hh[] ;
const int INF=0xfffffff ;
struct node{
    int s,t,cap,cost,nxt ;
}e[] ;
int sumflow ;
int n,m,cnt,head[],vis[],dis[],pre[] ;
char M[][] ;
void add(int s,int t,int cap,int cost)
{
    e[cnt].s=s ;e[cnt].t=t ;e[cnt].cap=cap ;e[cnt].cost=cost ;e[cnt].nxt=head[s] ;head[s]=cnt++ ;
    e[cnt].s=t ;e[cnt].t=s ;e[cnt].cap= ;e[cnt].cost=-cost ;e[cnt].nxt=head[t] ;head[t]=cnt++ ;
}
int spfa(int s,int t,int N)
{
    for(int i= ;i<=N ;i++)
        dis[i]=INF ;
    dis[s]= ;
    memset(vis,,sizeof(vis)) ;
    memset(pre,-,sizeof(pre)) ;
    vis[s]= ;
    queue <int> q ;
    q.push(s) ;
    while(!q.empty())
    {
        int u=q.front() ;
        q.pop() ;
        vis[u]= ;
        for(int i=head[u] ;i!=- ;i=e[i].nxt)
        {
            int tt=e[i].t ;
            if(e[i].cap && dis[tt]>dis[u]+e[i].cost)
            {
                dis[tt]=dis[u]+e[i].cost ;
                pre[tt]=i ;
                if(!vis[tt])
                {
                    vis[tt]= ;
                    q.push(tt) ;
                }
            }
        }
    }
    if(dis[t]==INF)return  ;
    return  ;
}
int MCMF(int s,int t,int N)
{
    int flow,minflow,mincost ;
    mincost=flow= ;
    while(spfa(s,t,N))
    {
        minflow=INF ;
        for(int i=pre[t] ;i!=- ;i=pre[e[i].s])
            minflow=min(minflow,e[i].cap) ;
        flow+=minflow ;
        for(int i=pre[t] ;i!=- ;i=pre[e[i].s])
        {
            e[i].cap-=minflow ;
            e[i^].cap+=minflow ;
        }
        mincost+=dis[t]*minflow ;
    }
    sumflow=flow ;//最大流
    return mincost ;
}
int main()
{
    while(scanf("%d%d",&n,&m))
    {
        if(!n && !m)break ;
        cnt= ;
        memset(head,-,sizeof(head)) ;
        for(int i= ;i<n ;i++)
            scanf("%s",M[i]) ;
        int cm,ch ;
        cm=ch= ;
        for(int i= ;i<n ;i++)
        {
            for(int j= ;j<m ;j++)
            {
                if(M[i][j]=='m')
                {
                    mm[cm].x=i ;
                    mm[cm++].y=j ;
                }
                if(M[i][j]=='H')
                {
                    hh[ch].x=i ;
                    hh[ch++].y=j ;
                }
            }
        }
        int SS= ;
        int ST=*ch+ ;
        for(int i= ;i<ch ;i++)
        {
            add(,i+,,) ;
            for(int j= ;j<ch ;j++)
            {
                int temp=ABS(mm[i].x-hh[j].x)+ABS(mm[i].y-hh[j].y) ;
                add(i+,j++ch,,temp) ;
            }
            add(i++ch,ST,,) ;
        }
        printf("%d\n",MCMF(SS,ST,ST+)) ;
    }
    return  ;
}