题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1533
大意: 在一个N*M的方格中, 有些人和房子, 让每个人都进入房子, 每个人每走一步的花费为1, 每个房子只能进一个人,但可以经过多次。问所有的人都进入房子的最小花费。
问题分析:这是网络流的最大流时的最小花费。这个建图模型就是经典的方格上的建图。我们可以建一个超级源点和超级汇点, 把源点和所有人都连一条边,容量为1, 费用为0, 把所有的房子和超级汇点都建一条边, 容量为1, 费用为0, 然后再把方格中的每一个小方格和其相邻的四个方向的小方格都建一条边,容量为无穷大, 费用为1,建图就完成了。
#include<stdio.h>
#include<string.h>
#include<vector>
#include<queue>
#include<algorithm>
using namespace std;
const int maxn = 10000 + 10;
const int inf = 0x3f3f3f3f;
int dir[2][4] = {{0, -1, 0, 1}, {-1, 0, 1, 0}};
struct Edge
{
int from, to, cap, flow, cost;
Edge(int u, int v, int c, int f, int s)
: from(u), to(v), cap(c), flow(f), cost(s) {}
};
struct MCMF
{
vector<int> G[maxn];
vector<Edge> edges;
bool inq[maxn];
int d[maxn], p[maxn], a[maxn];
int n, m, s, t;
void inin(int n)
{
this->n = n;
for(int i = 0; i <= n; i++)
G[i].clear();
edges.clear();
}
void AddEdge(int from, int to, int cap, int cost)
{
edges.push_back(Edge(from, to, cap, 0, cost));
edges.push_back(Edge(to, from, 0, 0, -cost));
m = edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}
bool spfa(int& flow, int& cost)
{
memset(inq, false, sizeof(inq));
memset(d, 0x3f, sizeof(d));
d[s] = 0, a[s] = inf, p[s] = s;
queue<int> Q;
Q.push(s);
inq[s] = true;
while(!Q.empty())
{
int x = Q.front();
Q.pop();
inq[x] = false;
for(int i = 0; i < G[x].size(); i++)
{
Edge& e = edges[G[x][i]];
if(e.cap>e.flow && d[e.to]>d[x]+e.cost)
{
d[e.to] = d[x] + e.cost;
p[e.to] = G[x][i];
a[e.to] = min(a[x], e.cap-e.flow);
if(!inq[e.to])
{
Q.push(e.to);
inq[e.to] = true;
}
}
}
}
if(d[t] == inf)
return false;
flow += a[t];
cost += d[t] * a[t];
int u = t;
while(u != s)
{
edges[p[u]].flow += a[t];
edges[p[u]^1].flow -= a[t];
u = edges[p[u]].from;
}
return true;
}
int min_cost(int s, int t)
{
this->s = s;
this->t = t;
int flow = 0, cost = 0;
while(spfa(flow, cost));
return cost;
}
};
bool Pan(int x, int y, int n, int m)
{
if(x>=0 && x<n && y>=0 &&y<m)
return true;
return false;
}
int main()
{
MCMF solve;
int n, m;
char Map[101][101];
while(scanf("%d%d", &n, &m), n||m)
{
int s = n*m, t = n*m+1;
solve.inin(t);
for(int i = 0; i < n; i++)
{
getchar();
for(int j = 0; j < m; j++)
scanf("%c", &Map[i][j]);
}
for(int i = 0; i < n; i++)
for(int j = 0; j < m; j++)
{
if(Map[i][j] == 'm')
solve.AddEdge(s, i*m+j, 1, 0);
else if(Map[i][j] == 'H')
solve.AddEdge(i*m+j, t, 1, 0);
for(int k = 0; k < 4; k++)
{
int x = i + dir[0][k];
int y = j + dir[1][k];
if(Pan(x, y, n, m))
solve.AddEdge(i*m+j, x*m+y, inf, 1);
}
}
printf("%d\n", solve.min_cost(s, t));
}
return 0;
}