专栏首页数据结构与算法洛谷P2762 太空飞行计划问题(最大权闭合图)

洛谷P2762 太空飞行计划问题(最大权闭合图)

题意

有$m$个实验,$n$中器材,每个实验需要使用一些器材

每个实验有收入,每个器材有花费

最大化收入 - 花费

Sol

最大权闭合图的经典应用

从$S$向每个实验连流量为该实验收入的边

从每个器材箱$T$连流量为花费的边

每个实验向其需要其器材连边权为$INF$的边

答案为:总收入 - 最小割

考虑如何统计方案

在最小割中,割去实验表示不选该实验。

那么我们从源点出发,不经过边权为$0$的边,走到的就是需要选的。

这正好是Dinic最后一次增光的deep数组

#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
const int MAXN = 1e5 + 10, INF = 1e9 + 10;
char c;
inline int read() {
    if(c == '\r') return 0;
    c = getchar(); int x = 0, f = 1; 
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int M, N, S, T;
struct Edge {
    int u, v, f, nxt;
}E[MAXN];
int head[MAXN], cur[MAXN], num;
inline void add_edge(int x, int y, int f) {
    E[num] = (Edge){x, y, f, head[x]};
    head[x] = num++;
}
inline void AddEdge(int x, int y, int z) {
    add_edge(x, y, z);
    add_edge(y, x, 0);
}
int sum = 0, deep[MAXN];
bool BFS() {
    queue<int> q; q.push(S);
    memset(deep, 0, sizeof(deep)); deep[S] = 1;
    while(!q.empty()) {
        int p = q.front(); q.pop();
        for(int i = head[p]; i != -1; i = E[i].nxt) {
            int to = E[i].v;
            if(!deep[to] && E[i].f) {
                deep[to] = deep[p] + 1;
                q.push(to);
            }
        }
    }
    return deep[T] > 0;
}
int DFS(int x, int flow) {
    if(x == T) return flow;
    int ansflow = 0;
    for(int &i = cur[x]; i != -1; i = E[i].nxt) {
        int to = E[i].v;
        if(deep[to] == deep[x] + 1 && E[i].f) {
            int nowflow = DFS(to, min(flow, E[i].f));
            E[i].f -= nowflow; E[i ^ 1].f += nowflow;
            ansflow += nowflow; flow -= nowflow;
            if(flow <= 0) break;
        }
    }
    return ansflow;
}
int Dinic() {
    int ans = 0;
    while(BFS()) {
        memcpy(cur, head, sizeof(head));
        ans += DFS(S, INF);
    }
    return ans;
}
int main() {
    memset(head, -1, sizeof(head));
    scanf("%d %d\n", &M, &N); S = 0; T = N + M + 1;
    int ans = 0;
    for(int i = 1; i <= M; i++) {
        c = '+';
        int val = read(), x; ans += val;
        AddEdge(S, i, val);
        while(x = read()) 
            AddEdge(i, x + M, INF);
    }
    for(int i = 1; i <= N; i++) {
        int x; scanf("%d", &x);
        AddEdge(i + M, T, x);
    }
    int cut = Dinic();
    for(int i = head[S]; i != -1; i = E[i].nxt)
        if(E[i].f)
            printf("%d ", E[i].v); puts("");
    for(int x = M + 1; x <= N + M; x++)
        for(int i = head[x]; i != -1; i = E[i].nxt)
            if(E[i].f)
                {printf("%d ", x - M); break;}
    puts("");
    printf("%d", ans - cut);
    return 0;
}

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • UOJ#310. 【UNR #2】黎明前的巧克力(FWT)

    attack
  • 洛谷P3763 [TJOI2017]DNA(后缀数组 RMQ)

    具体来说,每次比较当前后缀和\(S_0\)的lcp,如果长度\(< N\)的话就从不合法的位置继续匹配

    attack
  • P1018 乘积最大

    题目描述 今年是国际数学联盟确定的“2000――世界数学年”,又恰逢我国著名数学家华罗庚先生诞辰90周年。在华罗庚先生的家乡江苏金坛,组织了一场别开生面的数学智...

    attack
  • 【UVALive - 6534 】Join two kingdoms (树的直径的期望)

    给两棵树,分别有 n,m 个节点(1 ≤ N, Q ≤ 4 × 10^4),等概率连接属于不同树的两个节点,求新树的直径(最远两点的距离)的期望。

    饶文津
  • UOJ#310. 【UNR #2】黎明前的巧克力(FWT)

    attack
  • BUPT2017 wintertraining(15) #2 题解

    ​ 因为2520%pre_lcm0,所以x%pre_lcm(x%2520)%pre_lcm

    饶文津
  • 【GPLT】L2-005 集合相似度

    版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。 ...

    喜欢ctrl的cxk
  • P2015 二叉苹果树 树状dp

    用户2965768
  • 2013年第四届蓝桥杯C/C++B组省赛题目解析

    Zoctopus
  • poj-1218 THE DRUNK JAILER 喝醉的狱卒

    就是一个狱卒喝醉了,他第一趟吧所有的监狱都带开,第二趟把能把二整除的监狱关闭,第三趟操作能把三整除的监狱;

    瑾诺学长

扫码关注云+社区

领取腾讯云代金券