BOJ)10319 좀비 아포칼립스

문제: icpc.me/10319

감염이 되기전에 병원으로 이동가능한 최대한의 사람을 구하는 문제이다.

문제에서 요구하는 조건이 곧 시작지점부터 병원까지의 최대유량을 구하는 문제인데 

문제는 간선에 capacity말고도 단위시간이라는 제약 조건이 붙는다.

이걸 해결해주기 위해 단위시간만큼 정점을 분할 시켜주어 각 시간별로 주어진 정점들로 이루어진 그래프로 새로 재구성 시켜준 뒤 maximum flow를 구해주면 된다.

#include <cstdio>
#include <algorithm>
#include <vector>
#include <queue>
#define INF 987654321
using namespace std;
struct dinic {
    int n;
    struct Edge {
        int v, cap, rev;
        Edge(int v, int cap, int rev) :v(v), cap(cap), rev(rev) {}
    };
    vector<vector<Edge>> vt;
    vector<int> level, work;
    void addEdge(int s, int e, int c) {
        vt[s].emplace_back(e, c, vt[e].size());
        vt[e].emplace_back(s, 0, vt[s].size() - 1);
    }
    dinic(int n) :n(n) {
        vt.resize(n + 1);
    }
    bool bfs(int src, int sink) {
        level.assign(n + 1, -1);
        level[src] = 0;
        queue<int> qu;
        qu.push(src);
        while (qu.size()) {
            int here = qu.front();
            qu.pop();
            for (auto next : vt[here]) {
                if (next.cap&&level[next.v] == -1) {
                    level[next.v] = level[here] + 1;
                    qu.push(next.v);
                }
            }
        }
        return level[sink] != -1;
    }
    int dfs(int here, int crtcap, int target) {
        if (here == target)return crtcap;
        for (int &i = work[here]; i < vt[here].size(); i++) {
            int next = vt[here][i].v;
            int cap = vt[here][i].cap;
            if (cap&&level[here] + 1 == level[next]) {
                int flow = dfs(next, min(crtcap, cap), target);
                if (flow) {
                    vt[here][i].cap -= flow;
                    vt[next][vt[here][i].rev].cap += flow;
                    return flow;
                }
            }
        }
        return 0;
    }
    int stream(int src, int sink) {
        int totalflow = 0;
        while (bfs(src, sink)) {
            work.assign(n + 1, 0);
            while (1) {
                int flow = dfs(src, INF, sink);
                if (!flow)break;
                totalflow += flow;
            }
        }
        return totalflow;
    }
};
int t, n, src, g, s, m, r, x, y, z, w;
int main() {
    scanf("%d", &t);
    while (t--) {
        scanf("%d%d%d%d", &n, &src, &g, &s);
        dinic f((n + 1)*(s + 1) + 3);
        f.addEdge((n + 1)*(s + 1) + 1, src*(s + 1), g);
        for (int i = 1; i <= n; i++) {
            for (int j = 0; j < s; j++) 
                f.addEdge(i*(s + 1) + j, i*(s + 1) + j + 1, INF);
        }
        scanf("%d", &m);
        for (int i = 0; i < m; i++) {
            scanf("%d", &x);
            f.addEdge(x*(s + 1) + s, (n + 1)*(s + 1) + 2, INF);
        }
        scanf("%d", &r);
        for (int i = 0; i < r; i++) {
            scanf("%d%d%d%d", &x, &y, &z, &w);
            for (int i = 0; i + w <= s; i++)
                f.addEdge(x*(s + 1) + i, y*(s + 1) + i + w, z);
        }
        printf("%d\n", f.stream((s + 1)*(n + 1) + 1, (s + 1)*(n + 1) + 2));
    }
    return 0;
}