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
*/