CodeChef Course Selection（最小割）

题目链接

题意：
$\ \ \ \ $有$n$项课程，每项课程可能在$m$个学期中的其中至少一个学期开设，对应得到的分数$x$，每项课程必须也只能选一次，还有一些完成课程的前后关系，为完成这些课程最大的得分。

思路：
$\ \ \ \ $参考$16$年姜志豪的论文《网络流的一些建模方法》，最大化得分， 就是最小化失去的分， 用最小割模型：
$1$ 对于每个课程$i$，学期$j$，建立点$(i， j)$
$2$ 源点$S$连向$(i, 1)$一条边， 容量为$200 - x_{i,j}$，如果$x_{i,j}$为$-1$则设为$inf$
$3$ 对于$2 \leqslant j \leqslant m$,所有$(i,j- 1)$连向$(i,j)$，容量为$200 - x_{i,j}$
$4$ $(i,m)$连向汇点$T$，容量为$inf$
$5$ 对于有前后关系的课程$a$和$b$，源点$S$连向$(b, 1)$一条边， 容量为$inf$，对于$2 \leqslant j \leqslant m$,所有$(a,j- 1)$连向$(b,j)$，容量为$inf$。
$\ \ \ \ $那么求最小割的意义就是，割边代表课程选择时候损失的得分，这样损失越小得分越大，求下最小割$tot$，然后$200 * n - tot$就是最大得分了，平均值直接除以 n n <script type="math/tex" id="MathJax-Element-5880">n</script>即可。

#include<bits/stdc++.h>
typedef long long ll;
const int INF = 1e6 + 10;
const int maxn = 1e4 + 10;
using namespace std;

struct st {
    int to, cap, re;
    st(int t = 0, int c = 0, int r = 0) : to(t), cap(c), re(r) {}
};
int n, m, s, t, k, T, kase = 1;
vector<st> G[maxn];
int it[maxn], lv[maxn];

void add(int f, int t, int c) {
    //if(c == INF) printf("%d -> %d : INF\n", f, t);
   // else printf("%d -> %d : %d\n", f, t, c);
    G[f].push_back(st(t, c, G[t].size()));
    G[t].push_back(st(f, 0, G[f].size() - 1));
}

void bfs() {
    memset(lv, -1, sizeof(lv));
    queue<int> q;
    lv[s] = 0;
    q.push(s);
    while(!q.empty()) {
        int u = q.front(); q.pop();
        for(int i = 0; i < G[u].size(); i++) {
            st &e = G[u][i];
            if(e.cap > 0 && lv[e.to] < 0) {
                lv[e.to] = lv[u] + 1;
                q.push(e.to);
            }
        }
     }
}

int dfs(int v, int t, int f) {
    if(v == t) return f;
    for(int &i = it[v]; i < G[v].size(); i++) {
        st &e = G[v][i];
        if(e.cap > 0 && lv[v] < lv[e.to]) {
            int d = dfs(e.to, t, min(f, e.cap));
            if(d > 0) {
                e.cap -= d;
                G[e.to][e.re].cap += d;
                return d;
            }
        }
    }
    return 0;
}

int maxflow() {
    int f = 0;
    while(1) {
        bfs();
        if(lv[t] < 0) return f;
        memset(it, 0, sizeof(it));
        int fl;
        while((fl = dfs(s, t, INF)) > 0) f += fl;
    }
}

int x[110][110], mp[110][110];;
int a[maxn], b[maxn];

int solve() {
    s = 0; t = n * m + 1;
    for(int i = s; i <= t; i++) G[i].clear();
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            int to = mp[i][j], from = ((j == 1) ? s : mp[i][j - 1]);
            if(x[i][j] == -1) add(from, to, INF);
            else add(from, to, 200 - x[i][j]);
        }
        add(mp[i][m], t, INF);
    }
    for(int i = 1; i <= k; i++) {
        int tx = a[i], ty = b[i];
        add(s, mp[ty][1], INF);
        for(int j = 1; j < m; j++) add(mp[tx][j], mp[ty][j + 1], INF);
    }
    int tot = maxflow();
    //cout << "tot = " << tot << endl;
    return 200 * n - tot;
}


int main() {
    while(scanf("%d %d %d", &n, &m, &k) != EOF) {
        int cnt = 1;
        for(int i = 1; i <= n; i++)
            for(int j = 1; j <= m; j++) {
                scanf("%d", &x[i][j]);
                mp[i][j] = cnt++;
            }
        for(int i = 1; i <= k; i++) scanf("%d %d", &a[i], &b[i]);
        int ans = solve();
        printf("%.2f\n", (double)ans / n);
    }
    return 0;
}