poj 2195 最小费用最大流

///题意：n个人和m个房子，每个房子中要安排一个人，每个人移动一步费用为1，
    ///所有人安排房子后需要的金钱的最小值。

///曼哈顿距离：例如在平面上，坐标（x1, y1）的i点与坐标（x2, y2）的j点的曼哈顿距离为：
    ///  d(i,j)=|X1-X2|+|Y1-Y2|.

///建图：网络流，建一个超源 0，超汇  E
   /// 1> 超源  0  到 所有人建一条边  ，容量为  1   费用为  0
  ///  2> 所有人到所有房子  建一条边   容量为1，正向费用为用： abs(x人-x房子)+abs(y人-y房子)。
///反向为:   -正向
  ///  3> 所有房子  到  超汇(E) 建一条边 ， E = 所有人个数 + 所有边个数 + 1 。容量为 1 ，费用为 0

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <algorithm>
#include <vector>
#include <queue>
#define Max 210
using namespace std;
struct loc
{
    int x,y;
}coordinate[Max];
struct Edge
{
    int to,cap,cost,rev;  ///终点，容量，费用，反向边
};
vector<Edge> g[Max];  ///图的邻接表
queue<int> Q;
char map[Max][Max];
int dist[Max];   ///记录最短路径
int prevv[Max];  ///最短路的前驱点
int preve[Max];  ///前驱点的后继点
bool vis[Max];
int n,m;
int sh,sm;
int E;

///向图中增加一条从from到to容量为cap费用为cost的边
void add_edge(int from,int to,int cap,int cost)
{
    g[from].push_back((Edge){to,cap,cost,g[to].size()});
    g[to].push_back((Edge){from,0,-cost,g[from].size()-1});
}
int min_cost_flow(int f)///spfa求最短路，并记录路径
{
	int res=0;
    while(f>0)
    {
        memset(dist,0x3f,sizeof(dist));
        memset(vis,0,sizeof(vis));
        dist[0]=0;
        vis[0]=true;
        Q.push(0);
        while(!Q.empty())
        {
            int u=Q.front();
            Q.pop();
            vis[u]=false;
            for(int v=0;v<g[u].size();v++)
            {
                Edge &e=g[u][v];
                if(e.cap>0&&dist[e.to]>dist[u]+e.cost)
                {
                    dist[e.to]=dist[u]+e.cost;
                    prevv[e.to]=u;
                    preve[e.to]=v;
                    if(!vis[e.to])
                    {
                        Q.push(e.to);
                        vis[e.to]=1;
                    }
                }
            }
        }
		f-=1;
        res+=dist[E];
        for(int v=E;v!=0;v=prevv[v])    ///更改路径流量
        {
            Edge &e = g[prevv[v]][preve[v]];
            e.cap -= 1;
			g[v][e.rev].cap += 1;
            ///g[e.to][e.rev].cap += 1;
        }
    }
    return res;
}
int main()
{
	while(scanf("%d%d",&n,&m) && (n+m))
	{
		int i,j,k;
		getchar();
		sh = 0,sm = 0;
		for(i=0;i<Max;i++)
		{
			g[i].clear();
		}
		for(i=0;i<n;i++)
		{
			gets(map[i]);
			for(j=0;j<m;j++)
			{
				if(map[i][j] == 'm')/// 找到房子 计数
				{
					sm ++ ;
					coordinate[sm].x = i;
					coordinate[sm].y = j;
					add_edge(0,sm,1,0);///超级源点 连接  每个人
				}
			}
		}
		for(i=0;i<n;i++)
		{
			for(j=0;j<m;j++)
			{
				if(map[i][j] == 'H')
				{
					sh ++ ;
					for(k=1;k<=sh;k++)
					{
						add_edge(k,sh+sm,1,abs(coordinate[k].x-i)+abs(coordinate[k].y-j));///每一个人连接 每一个房子
					}
				}
			}
		}
		E = sh+sm+1;
		for(i=sh+1;i<E;i++)
		{
			add_edge(i,E,1,0);///每一个房子连接 超级汇点
		}
		int f = min(sh,sm);
		printf("%d\n",min_cost_flow(f));
	}
	return 0;
}

/*
2 2
.m
H.
5 5
HH..m
.....
.....
.....
mm..H
7 8
...H....
...H....
...H...
mmmHmmmm
...H....
...H....
...H....
0 0
*/