最小费用最大流

土豪红红突然兴致大发准备免费送大家一套房子，现在他有一张N行M列的地图，地图中m表示一个人，H表示一幢房子，每个房子只能被一个人选。红红决定为了方便起见让所有人到自己房子的距离和最小(人到房子的距离为人到房子所走的步数)，请问这个距离和是多少？
输入
输入包含多组测试样例。n和m都为零时表示程序结束。

每组测试样例输入第一行包含两个整数n和m，表示这个地图的行数和列数。

接下来n行m列表示红红的地图。

输出
输出题目要求的最小距离和。
样例输入
2 2
.m
H.
5 5
HH…m
…
…
…
mm…H
7 8
…H…
…H…
…H…
mmmHmmmm
…H…
…H…
…H…
0 0
样例输出
2
10
28

#include<bits/stdc++.h>
using namespace std;
const int maxn=1000;
const int INF=0x3f3f3f3f;

struct edge{
   int to,cap,w,next;
}edge[maxn*maxn];
int head[maxn];
int cnt;

char Map[maxn][maxn];
int n,m;

int vis[maxn],dis[maxn],pre[maxn];

struct info{
    int num,x,y;
};

 void init(){
    memset(head,-1,sizeof(head));
    cnt=0;
 }
void addedge(int u,int v,int c,int w){
    //正向边
     edge[cnt].to=v;
     edge[cnt].cap=c;
     edge[cnt].w=w;
     edge[cnt].next=head[u];
     head[u]=cnt++;
     //反向边
     edge[cnt].to=u;
     edge[cnt].cap=0;
     edge[cnt].w=-w;
     edge[cnt].next=head[v];
     head[v]=cnt++;
}

bool spfa(int s,int e){

   memset(vis,0,sizeof(vis));
   memset(dis,INF,sizeof(dis));
   memset(pre,-1,sizeof(pre));
   queue<int> q;
   q.push(s);
   vis[s]=1;
   dis[s]=0;
   while(!q.empty()){
      int u=q.front();
      q.pop();
      vis[u]=0;
      for(int i=head[u];~i;i=edge[i].next){
          int v=edge[i].to;
          int cap=edge[i].cap;
          if(dis[v]>dis[u]+edge[i].w&&cap>0){
              dis[v]=dis[u]+edge[i].w;
              pre[v]=i;
              if(!vis[v]){
                vis[v]=1;
                q.push(v);
              }
          }
      }
   }
   if(dis[e]==INF)
      return false;
    else
      return true;
}

int MincostMaxflow(int s,int e,int &cost){
    int maxflow=0;
    while(spfa(s,e)){
      int minflow=INF;
        for(int i=pre[e];i!=-1;i=pre[edge[i^1].to])  //要求正向边的起点，就求反向边的终点  ^1 偶数+1，奇数-1
          minflow=min(minflow,edge[i].cap);
        for(int i=pre[e];i!=-1;i=pre[edge[i^1].to]){
          edge[i].cap-=minflow;
          edge[i^1].cap+=minflow;
          cost+=minflow*edge[i].w;
        }
        maxflow+=minflow;
    }
    return maxflow;
}

int main(){
    while((cin>>n>>m)&&(n||m)){
        init();
        int start=0,end=n*m+1;
        for(int i=0;i<n;i++){
           scanf("%s",Map[i]);
        }
        vector<info> people,house;
        for(int i=0;i<n;i++){
           for(int j=0;j<m;j++){
              if(Map[i][j]=='.')  continue;
              if(Map[i][j]=='H'){
                 int num=i*m+j+1;
                 people.push_back((info){num,i,j});
              }
              if(Map[i][j]=='m'){
                 int num=i*m+j+1;
                 house.push_back((info){num,i,j});
              }
           }
        }
         for(int i=0;i<people.size();i++){
             info ipeople=people[i];
             addedge(start,ipeople.num,1,0);
         }
         for(int i=0;i<house.size();i++){
             info ihouse=house[i];
             addedge(ihouse.num,end,1,0);
         }
        for(int i=0;i<people.size();i++){
            for(int j=0;j<house.size();j++){
              info ipeople=people[i];
              info ihouse=house[j];
              int dis=fabs(ipeople.x-ihouse.x)+fabs(ipeople.y-ihouse.y);
              addedge(ipeople.num,ihouse.num,1,dis);
            }
        }
        int cost=0;
        MincostMaxflow(start,end,cost);
        cout<<cost<<endl;
    }
    return 0;
}