P2762 太空飞行计划问题

题目描述

有m个实验，每个实验只可以进行一次，但会获得相应的奖金，有n个仪器，每个实验都需要一定的仪器，每个仪器可以运用于多个实验，但需要一定的价值，问奖金与代价的差的最大值是多少？

思路

最大权闭合图模型题，最大流最小割。从起点往每个权值为正的点建立一条边，容量为点权值。每个权值为负的点往终点建立一条边，容量为权值的绝对值。如果选A就必须选B 则就从A建立一条往B的边，容量为正无穷。然后正权值加起来减去最大流就是能选出来最大权闭合图所有点加起来的值，最大权闭合图的点就是从起点开始广搜，权值为0的点不走，能走到的点就是被选中的点。dinic最后一次bfs的d数组正好可以用来判断这个条件。

洛谷上有个图不错：

跑最大流，根据最大流等于最小割的定理来求出最小割，最小割就是买仪器和亏本的实验的金钱的和，用总的赞助减去最小割就是净利润，跑BFS最后一遍出来的dep数组中的值不为0的话，就是选择了这个实验（这个实验净利润>0）和实验需要的仪器（实验与仪器连的边为INF，正好选上）

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define ll long long
using namespace std;
const int N = 1e5 + 10;
int n, m, s, t, h[N], cnt, sum, dep[N], cur[N];
char c;
struct node {
    int v, w, nt;
} no[N];
void add(int u, int v, int w) {
    no[cnt] = node{v, w, h[u]};
    h[u] = cnt++;
}
int bfs() {
    queue<int> q;
    memset(dep, 0, sizeof dep);
    dep[s] = 1;
    q.push(s);
    while(!q.empty()) {
        int u = q.front();
        q.pop();
        for(int i = h[u]; ~i; i = no[i].nt) {
            int v = no[i].v;
            if(!dep[v] && no[i].w) {
                dep[v] = dep[u] + 1;
                q.push(v);
            }
        }
    }
    return dep[t] > 0;
}
int dfs(int u, int flow) {
    if(u == t)
        return flow;
    for(int i = cur[u]; ~i; i = no[i].nt) {
        int v = no[i].v;
        if(dep[v] == dep[u] + 1 && no[i].w) {
            int res = dfs(v, min(flow, no[i].w));
            if(res > 0) {
                no[i].w -= res;
                no[i ^ 1].w += res;
                return res;
            }
        }
    }
    return 0;
}
int dinic() {
    int res = 0;
    while(bfs()) {
        for(int i = 0; i <= cnt; i++)
            cur[i] = h[i];
        while(int d = dfs(s, INF))
            res += d;
    }
    return res;
}
int main() {
    memset(h, -1, sizeof h);
    scanf("%d%d", &n, &m);
    t = n + m + 1;
    for(int x, i = 1; i <= n; i++) {
        scanf("%d", &x);
        sum += x;
        add(s, m + i, x);
        add(m + i, s, 0);
        while(1) {
            scanf("%d%c", &x, &c);
            add(m + i, x, INF);
            add(x, m + i, 0);
            if(c == '\n' || c == '\r')
                break;
        }
    }
    for(int w, i = 1; i <= m; i++) {
        scanf("%d", &w);
        add(i, t, w);
        add(t, i, 0);
    }
    sum -= dinic();
    for(int i = 1; i <= n; i++)
        if(dep[i + m])
            printf("%d ", i);
    puts("");
    for(int i = 1; i <= m; i++)
        if(dep[i])
            printf("%d ", i);
    puts("");
    printf("%d\n", sum);
    return 0;
}