标签：
Going Home28
题意：m个人和m个房子，每个房子中要安排一个人，每个人移动一步费用为1，
所有人安排房子后需要的金钱的最小值。
建图：网络流，建一个超源 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 <queue>
using namespace std;
#define Max 202
#define INF 0x3f3f3f3f
bool vis[Max];
int dist[Max];
int pre[Max];
int cost[Max][Max];
int cap[Max][Max];
int E, sh, sm, res;
struct Node
{
int x, y;
}hos[Max], man[Max];
bool spfa()//求最短路
{
int i, u;
memset(dist,0x3f,sizeof(dist));
memset(vis,false,sizeof(vis));
queue <int> q;
dist[0] = 0;
vis[0] = true;
while (!q.empty())
q.pop();
q.push(0);
while (!q.empty())
{
u = q.front();
q.pop();
vis[u] = false;
for (i = 1; i <= E; i++)
{
if (cap[u][i] && dist[i] > dist[u] + cost[u][i])
{
dist[i] = dist[u] + cost[u][i];
pre[i] = u;
if (!vis[i])
{
vis[i] = true;
q.push(i);
}
}
}
}
if(dist[E] != INF)
{
return true;
}
return false;
}
void change()//更新 流量
{
int i, Min;
Min = INF;
for (i = E; i != 0; i = pre[i])
{
Min = min(Min, cap[pre[i]][i]);
}
for (i = E; i != 0; i = pre[i])
{
cap[pre[i]][i] -= Min;
cap[i][pre[i]] += Min;
res += cost[pre[i]][i] * Min;
}
}
int main()
{
int i, j, n, m;
while(scanf("%d%d",&n,&m) && (n+m))
{
char c;
sh = 0; sm = 0;
memset(cap, 0, sizeof(cap));
for (i = 0; i < n; i++)
{
for (j = 0; j < m; j++)
{
scanf("\n%c", &c);
if (c == 'H')
{
sh++;
hos[sh].x = i;
hos[sh].y = j;
}
else if (c == 'm')
{
sm++;
man[sm].x = i;
man[sm].y = j;
}
}
}
E = sh + sm + 1;//超汇
for (i = 1; i <= sm; i++) //超源为0 0到 所有人
{
cap[0][i] = 1; //容量为 1
cost[0][i] = 0; //费用为 0
}
for (i = 1; i <= sm; i++) //每个人到 每个房子
{
for (j = 1; j <= sh; j++)
{
cap[i][sm+j] = 1;
cost[i][sm+j] = abs(man[i].x-hos[j].x) + abs(man[i].y-hos[j].y); //正向边费用
cost[sm+j][i] = -cost[i][sm+j]; // 反向边费用
}
}
for (i = 1; i <= sh; i++) //所有房子 到 E 超汇为E
{
cap[sm+i][E] = 1;
cost[sm+i][E] = 0;
}
res = 0;
while(spfa())
{
change();
}
printf("%d\n", res);
}
return 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;
int u,v;
while(f>0)
{
memset(dist,0x3f,sizeof(dist));
memset(vis,0,sizeof(vis));
dist[0]=0;
vis[0]=1;
Q.push(0);
while(!Q.empty())
{
u=Q.front();Q.pop();vis[u]=0;
for(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(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] == 'H')
{
sh ++ ;
coordinate[sh].x = i;
coordinate[sh].y = j;
add_edge(0,sh,1,0);//超级源点 连接 每个人
}
}
}
for(i=0;i<n;i++)
{
for(j=0;j<m;j++)
{
if(map[i][j] == 'm')
{
sm ++ ;
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;
}
标签：
原文地址：http://blog.csdn.net/huangshuai147/article/details/51338695
