티스토리 뷰

알고리즘 관련/BOJ

BOJ)3640 제독

JASON 자손9319 2017. 3. 9. 15:41

문제: icpc.me/3640


경로가 겹치지 않게 두가지 방법으로 목적지 까지 갈 때 사용되는 최소 비용을 구하는 문제이다.


경로가 겹치지 않는 두가지 방법을 구하는 모델링은 네트워크 플로우 이론에서 정점 분할을 한 뒤 간선과 정점 사이의 모든 capacity를 1으로 준 뒤 플로우를 흘려주는 것으로 가능하다. 


이 때 두가지 방법으로 이동할 때 드는 최소 비용을 구하는 문제이므로 MCMF를 이용하여 minimum cost를 구하되 sink까지 플로우를 두번만 흘려주면 된다.



1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
#include <cstdio>
#include <algorithm>
#include <vector>
#include <queue>
#define INF 1987654321
using namespace std;
int n, m, x, y, z;
struct Edge {
    int v, cost, cap, rev;
    Edge(int v,int cost,int cap,int rev)
        :v(v),cost(cost),cap(cap),rev(rev){}
};
vector<vector<Edge>> vt;
vector<int> pv, pe;
void addEdge(int here, int there, int cost, int cap) {
    vt[here].emplace_back(there, cost, cap, vt[there].size());
    vt[there].emplace_back(here, -cost, 0, vt[here].size() - 1);
}
bool spfa(int src, int sink) {
    vector<int> v(* n + 30);
    vector<int> dist(* n + 3, INF);
    queue<int> qu;
    qu.push(src);
    dist[src] = 0;
    v[src] = 1;
    while (qu.size()) {
        int here = qu.front();
        v[here] = 0;
        qu.pop();
        for (int i = 0; i < vt[here].size(); i++) {
            int there = vt[here][i].v;
            int cap = vt[here][i].cap;
            if (cap&&dist[there] > dist[here] + vt[here][i].cost) {
                dist[there] = dist[here] + vt[here][i].cost;
                pv[there] = here;
                pe[there] = i;
                if (!v[there]) {
                    v[there] = 1;
                    qu.push(there);
                }
            }
        }
    }
    return dist[sink] != INF;
}
int solve(int src, int sink) {
    int ret = 0;
    while (spfa(src, sink)) {
        int minFlow = INF;
        for (int i = sink; i != src; i = pv[i]) {
            int prev = pv[i];
            int idx = pe[i];
            minFlow = min(minFlow, vt[prev][idx].cap);
        }
        for (int i = sink; i != src; i = pv[i]) {
            int prev = pv[i];
            int idx = pe[i];
            vt[prev][idx].cap -= minFlow;
            vt[i][vt[prev][idx].rev].cap += minFlow;
            ret += vt[prev][idx].cost*minFlow;
        }
    }
    return ret;
}
int main() {
    while (scanf("%d%d"&n, &m) != EOF) {
        int src = * n + 1, sink = * n + 2;
        vt.clear();
        vt.resize(* n + 3);
        pv.assign(* n + 3-1);
        pe.assign(* n + 3-1);
        for (int i = 1; i <= n; i++) {
            if (i == || i == n)
                addEdge(i, i + n, 02);
            else
                addEdge(i, i + n, 01);
        }
        for (int i = 0; i < m; i++) {
            scanf("%d%d%d"&x, &y, &z);
            addEdge(x + n, y, z, 1);
        }
        addEdge(src, 102);
        addEdge(n, sink, 02);
        printf("%d\n", solve(src, sink));
    }
    return 0;
}
cs


'알고리즘 관련 > BOJ' 카테고리의 다른 글

BOJ)11407 책 구매하기3  (0) 2017.03.11
BOJ)14428 수열과 쿼리 16  (0) 2017.03.09
BOJ)3640 제독  (0) 2017.03.09
BOJ)11409 열혈강호 6  (0) 2017.03.09
BOJ)11408 열혈강호5  (0) 2017.03.09
BOJ)9370 미확인 도착지  (0) 2017.03.08
TAG
댓글
댓글쓰기 폼