UVA 10746 Crime Wave - The Sequel

最小费用最大流 源点->警察->bank->汇点

剩下的模板就可以

#include <map>
#include <set>
#include <list>
#include <cmath>
#include <ctime>
#include <deque>
#include <stack>
#include <queue>
#include <cctype>
#include <cstdio>
#include <string>
#include <vector>
#include <climits>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define LL long long
#define PI 3.1415926535897932626
using namespace std;
int gcd(int a, int b) {return a % b ==  ? b : gcd(b, a % b);}
#define MAXN 50
#define INF 100000000
const double eps = 1e- ;
int flow[MAXN][MAXN],cap[MAXN][MAXN];
double cost[MAXN][MAXN];
int N,M,src,tag;
double c;
int ans;
void read()
{
    memset(cap,,sizeof(cap));
    memset(flow,,sizeof(flow));
    for (int i = ; i < MAXN; i++) for (int j = ; j < MAXN; j++) cost[i][j] = INF;
    src =  ; tag = N + M + ;
    for (int i = ; i <= M; i++) {cap[][i] = ; cost[][i] = ;}
    for (int i = ; i <= N; i++) {cap[M + i][tag] = ; cost[M + i][tag] = ;}
    for (int i = ; i <= N; i++)
        for (int j = ; j <= M; j++)
    {
        double tmp;
        scanf("%lf",&tmp);
        cost[j][i + M] = tmp;
        cost[i + M][j] = -tmp;
        cap[j][i + M] = ;
    }
}
void slove()
{
    queue<int>q;
    bool inq[MAXN];
    int p[MAXN];
    double d[MAXN];
    ans = ;
    c = 0.0;
    while (true)
    {
        memset(inq,false,sizeof(inq));
        for (int i = ; i < MAXN; i ++) d[i] = INF;
        d[src] = ;
        q.push(src);
        while (!q.empty())
        {
            int u = q.front(); q.pop();
            inq[u] = false;
            for (int v =  ; v <= tag; v++)
              if (cap[u][v] > flow[u][v] && d[v] > d[u] + cost[u][v] + eps)
              {
                  d[v] = d[u] + cost[u][v];
                  p[v] = u;
                  if (!inq[v])
                  {
                      inq[v] = true;
                      q.push(v);
                  }
              }
        }
        if (d[tag] == INF)  break;
        int a = INF;
        for (int u = tag; u != src; u = p[u])  a = min(a,cap[p[u]][u] - flow[p[u]][u]);
        for (int u = tag; u != src; u = p[u])
        {
            flow[p[u]][u] += a;
            flow[u][p[u]] -= a;
        }
        c += d[tag] * a;
        ans += a;
    }
}
int main()
{
    //freopen("sample.txt","r",stdin);
    //freopen("output.txt","w",stdin);
    while (scanf("%d%d",&N,&M) != EOF)
    {
        if (N ==  && M == ) break;
        read();
        slove();
        printf("%.2lf\n",c / N + eps);
    }
    return ;
}