uva 10746 Crime Wave - The Sequel（最小费用流）

题目：

链接：点击打开链接

题意：

一天n个银行被抢了，m个值班的警车在不同的位置，n个这样的警车应该被派往每个银行，以便平均到达n个银行的时间最少。

思路：

尽管题目很裸，还是不会做，看来还是差很多啊，。

单源，源点s=0,汇点t=n+m+1.警车编号1到m。银行编号m+1到n+m+1。

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
#define INF 100000000

const int N = 50;
const double eps = 1e-8;

//cap表示边(u,v)的容量。flow表示实际运送的流量。cost为时间花费
int n,m;
int cap

,flow

;
double cost

;
int s,t;

void init()
{
double c;
s = 0;
t = n + m + 1;
memset(cap,0,sizeof(cap));
memset(flow,0,sizeof(flow));
for(int i=0; i<=t; i++)
{
for(int j=0; j<=t; j++)
cost[i][j] = INF;
}
for(int i=1; i<=m; i++)//源点和警车相连
{
cap[s][i] = 1;
cost[s][i] = 0;
}
for(int i=m+1; i<t; i++)//银行和汇点相连
{
cap[i][t] = 1;
cost[i][t] = 0;
}
for(int i=m+1; i<t; i++)
{
for(int j=1; j<=m; j++)
{
scanf("%lf",&c);
cap[j][i] = 1;
cost[i][j] = -c;
cost[j][i] = c;
}
}
}

void EK()
{
queue<int> q;
int p
;
int a;
int vis
;
double d
,c = 0;
for(;;)
{
memset(vis,0,sizeof(vis));
for(int i=0; i<=t; i++)
d[i] = INF;
d[s] = 0;
q.push(s);
while(!q.empty())//BFS找增广路
{
int u = q.front();
q.pop();
vis[u] = 0;
for(int v=s; v<=t; v++)
{
if(cap[u][v] > flow[u][v] && d[v] > d[u] + cost[u][v] + eps)//找到新结点v
{
d[v] = d[u] + cost[u][v];
p[v] = u;//记录v的父亲，
if(!vis[v])
{
q.push(v);//并加入FIFO队列
vis[v] = 1;
}
}
}
}
if(d[t] == INF)
break;
a = INF;
for(int u=t; u != s; u=p[u])//从汇点往回走
{
a = min(a,cap[p[u]][u] - flow[p[u]][u]);//a为s-v路径上的最小残量
}
for ( int u = t; u != s; u = p[u] )
{
flow[p[u]][u] += a;//更新正向流量
flow[u][p[u]] -= a;//更新反向流量
}
c += d[t] * a;//从S流出的总流量
}
c /= n;//平均
printf("%.2lf\n",c + eps);
}

int main()
{
//freopen("input.txt","r",stdin);
while(scanf("%d%d",&n,&m) != EOF && (n+m))
{
init();
EK();
}
return 0;
}

-----------------------------------------------------------------------

收获：

最小费用最大流的建图：

cap边的容量，flow为实际的流量，还有花费cost。。。。。。。

应运EK增广路算法是要结合题目：

。。。。。。。。。。。。。。。。。。。。。。。

---------------------------------------------------------------------------

战斗，永不退缩；奋斗，永不停歇、、、、、、、、、、