北京集训队 2016 Day4 alarm

原题目

问题描述

  小H很喜欢玩电子游戏。在通关了“超级跳”以后,他又开始玩别的游戏。
  这次,小H对一个经营策略类游戏产生了兴趣。
  游戏中,小H控制一个平面营地,在营地上进行建设和发展,还要防备来自敌方的打击。尤其是要防备敌方的空袭。营地可以抽象成一个无向图。图中的有些节点是小H建造的防空避难所,其他的则是住房、农田、矿场等种种建筑,节点之间通过一些道路连接。每当空袭来临时,小H要将所有的农民(包括工人等)从各自所在的工作地点沿着道路转移到防空避难所中,以减少人员伤亡。
  经过许多次尝试,小H已经估算出从空袭警报到空袭降临的平均时间间隔,并据此重新部署了防空避难所。现在,已知每个避难所的人员容量,经过每条道路所需的时间,以及在每个节点工作的农民数目。假设任意数量的人可以同时在一条道路上移动,并假设所有的农民在空袭警报发出时同时开始移动。小H想知道将所有农民都从工作地点转移至防空所所需的最少时间。

输入格式

  第一行,两个整数n,m,表示无向图的点数和边数。
  接下来n行,每行两个整数wi和ci,依次描述节点1~n。wi表示每个节点上的农民数;如果ci =0,表示这不是一个防空避难所;否则ci表示防空避难所的人员容量。
  接下来m行,每行三个整数ui,vi和ti,表示存在一条节点ui和节点vi之间的道路,任何人经过这条道路的时间耗费为ti。

输出格式

  输出一行,一个整数,表示最少要花费的时间。如果不可能将所有农民都分配到防空避难所中,输出-1。

样例输入

3 3
5 0
0 4
2 6
1 2 4
1 3 12
2 3 7

样例输出

11

样例说明

  空袭警报来临时,节点3上的两个农民可以立刻躲进防空避难所中,这样3号点上的防空所的容量还剩4个人。节点1上的5个农民先花费4单位时间移动到点2。其中的4个人可以躲进防空所中,剩余的1个人继续移动到点3并躲进防空所。
  最后一个人需要11单位的时间才能躲进防空所。而这也是最优的方案。

数据规模和约定

  20% 的数据中,1 ≤ n≤5,wi之和不超过10。
  60% 的数据中,1 ≤ n≤50.
  100%的数据中,1 ≤n≤200,1≤m, ti≤10^4,0≤wi, ci≤10^4。无向图没有自环,但可能有重边。

题目大意

给定一个无向图,在图上有若干难民和若干避难所,每个人通过一条边需要 ti 的时间,一条边同时可以通过无数人。问在最短多少时间内能让所有难民进入避难所。

解答

首先看数据规模,n极小,且答案也很小,又由于题目的“从哪到哪”很像网络流,所以想到二分答案+网络流。
首先将两点间距离用Floyd预处理出来,时间复杂度为 O(n3)
之后二分答案,根据答案在每一对能在答案时间内到达的点对(难民位置,避难所)建边,再利用供给需求问题的思路建源到难民和避难所到汇的边,之后跑一遍Dinic,如果Dinic的结果恰好等于难民总数,则当前答案可以到达,把右边界降下来,如果不能到达,把左边界降下来。
【注意】:
(1)可能有重边,在建邻接矩阵时注意不要直接赋值而是比较大小后赋值。
(2)需要拆点

参考代码

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>

using namespace std;

const int oo = 2147483600;

class Edge
 {
    public:
        int to;
        int len;
        int next;
};

Edge es[100000];
int head[500], h[500];
int dis[250][250];
int et = -1;
int n, m;
int c[400], w[400];
int sumw;
int fl[500];
int st, th;

void addEdge(int u, int v, int l) {
    et++;
    es[et].to = v;
    es[et].len = l;
    es[et].next = h[u];
    h[u] = et;
}

void Floyd() {
    for (int k = 1; k <= n; k++)
        for (int i = 1; i <= n; i++)
            for (int j = 1; j <= n; j++)
                dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
}

void readin() {
    scanf("%d %d", &n, &m);
    for (int i = 1; i <= n; i++) {
        scanf("%d %d", &w[i], &c[i]);
        sumw += w[i];
    }
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= n; j++)
            if (i == j)
                dis[i][j] = 0;
            else
                dis[i][j] = 2000000;
    int x, y, z;
    for (int i = 0; i < m; i++) {
        scanf("%d %d %d", &x, &y, &z);
        dis[x][y] = dis[y][x] = min(z, dis[y][x]);
    }
}

bool bfs() {
    queue<int> ser;
    memset(fl, -1, sizeof(fl));
    fl[st] = 0;
    ser.push(st);
    int rt;
    bool flag = false;
    while (!ser.empty()) {
        rt = ser.front();
        for (int i = h[rt]; i != -1; i = es[i].next) {
            if (fl[es[i].to] == -1 && es[i].len) {
                fl[es[i].to] = fl[rt] + 1;
                ser.push(es[i].to);
                if (es[i].to == th) {
                    flag = true;
                }
            }
        }
        ser.pop();
    }
    return flag;
}

int dfs(int x, int flow) {
    if (x == th) {
        return flow;
    }
    int cflow = flow;
    int ret;
    for (int &i = head[x]; i != -1; i = es[i].next)
        if (fl[es[i].to] == fl[x] + 1 && es[i].len) {
            ret = dfs(es[i].to, min(es[i].len, flow));
            flow -= ret;
            es[i].len -= ret;
            es[i^1].len += ret;
            if (flow == 0)
                return cflow;
        }
    return cflow - flow;
}

int dinic()
 {
    int ans = 0;
    while (bfs()) {
        for (int i = 0; i <= th; i++)
            head[i] = h[i];
        ans += dfs(st, oo);
    }
    return ans;
}

bool check(int mt) {
    memset(h, -1, sizeof(h));
    et = -1;
    for (int i = 1; i <= n; i++) {
        if (w[i]) {
            addEdge(0, i*2, w[i]);
            addEdge(i*2, 0, 0);
        }
        if (c[i]) {
            addEdge(i*2+1, th, c[i]);
            addEdge(th, i*2+1, 0);
        }
        if (w[i])
            for (int j = 1; j <= n; j++)
                if (c[j] && dis[i][j] <= mt) {
                    addEdge(i*2, j*2+1, oo);
                    addEdge(j*2+1, i*2, 0);
                }
    }
    if (dinic() == sumw)
        return true;
    else
        return false;
}

void work() {
    int al = 0, ar = 1000000;
    int mid;
    st = 0, th = 2*n+3;
    while (al < ar) {
        mid = (al + ar) >> 1;
        if (check(mid)) {
            ar = mid;
        } else {
            al = mid+1;
        }
    }
    if (ar == 10000 && al == 10001)
        printf("-1");
    else
        printf("%d", ar);
}

int main() {
    freopen("alarm.in", "r", stdin);
    freopen("alarm.out", "w", stdout);
    readin();
    Floyd();
    work();
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值