网络流——最小割求最大权闭合子图

定义

有一个有向图,每一个点都有一个权值(可以为正或负或0),选择一个权值和最大的子图,使得每个点的后继都在子图里面,这个子图就叫最大权闭合子图。 
如下图: 
 
能选的子图有?,{4},{3,4},{2,4},{1,2,3,4},它们的权值分别为0,-1,5,-6,4. 
所以最大权闭合子图为{3,4},权值为5.

解法

这个问题可以转化为最小割问题,用网络流解决。 
从源点s向每个正权点连一条容量为权值的边,每个负权点向汇点t连一条容量为权值的绝对值的边,有向图原来的边容量全部为无限大。 
 
求它的最小割,割掉后,与源点s连通的点构成最大权闭合子图,权值为(正权值之和-最小割)。

如何理解

割掉一条边的含义

由于原图的边都是无穷大,那么割边一定是与源点s或汇点t相连的。

割掉s与i的边,表示不选择i点作为子图的点; 
割掉i与t的边,表示选择i点为子图的点。

如果s与i有边,表示i存在子图中; 
如果i与t有边,表示i不存在于子图中。

合法性

只有s与t不连通时,才能得到闭合子图。

如果s与t连通,则存在点i,j,使得s到i有边,i到j连通,j到t有边,所以j一定是i的后继,但选择了i,没有选择j,不是闭合子图。

如果s与t不连通,选择了正权点i,一定选择了i后继中的所有负权点。设j是i的后继中的正权点,则割掉s到j的边是没有意义的,最小割不会割掉它,则j一点被选中,所以i的所有后继都被选中,符合闭合图的定义。

最优性

最小割=(不选的正权之和+要选的负权绝对值之和) 
最大权闭合子图=(正权之和-不选的正权之和-要选的负权绝对值之和)=正权值和-最小割 
因为正权值和,是定值,而最小割保证值最小,所以最大权闭合子图一定最优。

例题

POJ2987_Firing

AC代码

#include <iostream>
#include <vector>
#include <queue>
using namespace std;

typedef long long LL;
#define MAX_V 5000 + 16

// 用于表示边的结构体(终点、容量、反向边)
struct edge
{
    int to, rev;
    LL cap;
    edge(int to, LL cap, int rev) :to(to), cap(cap), rev(rev){}
};

vector<edge> G[MAX_V];    // 图的邻接表表示
int level[MAX_V];        // 顶点到源点的距离标号
int iter[MAX_V];        // 当前弧,在其之前的边已经没有用了

// 向图中加入一条从from到to的容量为cap的边
void add_edge(int from, int to, int cap)
{
    G[from].push_back(edge(to, cap, G[to].size() ));
    G[to].push_back(edge(from, 0, G[from].size() - 1));
}

// 通过BFS计算从源点出发的距离标号
void bfs(int s)
{
    memset(level, -1, sizeof(level));
    queue<int> que;
    level[s] = 0;
    que.push(s);
    while (!que.empty())
    {
        int v = que.front(); que.pop();
        for (int i = 0; i < G[v].size(); ++i)
        {
            edge& e = G[v][i];
            if (e.cap > 0 && level[e.to] < 0)
            {
                level[e.to] = level[v] + 1;
                que.push(e.to);
            }
        }
    }
}

// 通过DFS寻找增广路
LL dfs(int v, int t, LL f)
{
    if (v == t)
    {
        return f;
    }
    for (int& i = iter[v]; i < G[v].size(); ++i)
    {
        edge& e = G[v][i];
        if (e.cap > 0 && level[v] < level[e.to])
        {
            LL d = dfs(e.to, t, min(f, e.cap));
            if (d > 0)
            {
                e.cap -= d;
                G[e.to][e.rev].cap += d;
                return d;
            }
        }
    }

    return 0;
}

// 求解从s到t的最大流
LL max_flow(int s, int t)
{
    LL flow = 0;
    for (;;)
    {
        bfs(s);
        if (level[t] < 0)
        {
            return flow;
        }
        memset(iter, 0, sizeof(iter));
        LL f;
        while ((f = dfs(s, t, 0x3f3f3f3f3f3f3f3f)) > 0)
        {
            flow += f;
        }
    }
}

int vertex_count, visited[MAX_V];
// 遍历残余网络
void solve(int v)
{
    ++vertex_count;
    visited[v] = true;
    for (int i = 0; i < int(G[v].size()); ++i)
    {
        const edge &e = G[v][i];
        if (e.cap > 0 && !visited[e.to])
        {
            solve(e.to);
        }
    }
}
///////////////////////////SubMain//////////////////////////////////
int main(int argc, char *argv[])
{
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
#endif
    int n, m, w;
    LL W = 0;
    scanf("%d%d", &n, &m);
    const int s = 0, t = n + 1;
    for (int i = 1; i <= n; i++)
    {
        scanf("%d", &w);
        if (w > 0)
        {
            W += w;
            add_edge(s, i, w);
        }
        if (w < 0)
        {
            add_edge(i, t, -w);
        }
    }

    int u, v;
    for (int i = 0; i < m; ++i)
    {
        scanf("%d%d", &u, &v);
        add_edge(u, v, 0x3f3f3f3f3f3f3f3f);
    }

    LL max_profit = W - max_flow(s, t);
    solve(s);
    printf("%d %I64d\n", --vertex_count, max_profit);
#ifndef ONLINE_JUDGE
    fclose(stdin);
    fclose(stdout);
    system("out.txt");
#endif
    return 0;
}

原文地址:https://www.cnblogs.com/shuaihui520/p/9157921.html

时间: 2024-10-25 19:31:27

网络流——最小割求最大权闭合子图的相关文章

BZOJ1565 NOI 2009 植物大战僵尸 topo+最小割(最大权闭合子图)

题目链接:https://www.luogu.org/problemnew/show/P2805(bzoj那个实在是有点小小的辣眼睛...我就把洛谷的丢出来吧...) 题意概述:给出一张有向图,这张有向图上的每个点都有一个点权,想要访问某个点必须要先访问这个点所能够访问(遍历)到的所有点,在访问到一个点之后将会得到这个点的权值(可正可负).问访问这张图可以得到的最大点权和. 原题说过来说过去实际上是描述了一个植物之间的保护关系,也就是说明了植物之间的先后访问顺序之间的关系.可以描述为要"要访问点

【HDU 5855】Less Time, More profit(网络流、最小割、最大权闭合子图)

Less Time, More profit Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Problem Description The city planners plan to build N plants in the city which has M shops. Each shop needs products from some plants to make p

HDU 4971 A simple brute force problem.(最小割,最大权闭合图)

http://acm.hdu.edu.cn/showproblem.php?pid=4971 A simple brute force problem. Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 182    Accepted Submission(s): 115 Problem Description There's a com

[HihoCoder1398]网络流五&#183;最大权闭合子图

题目大意:有$N$项活动$M$个人,每个活动$act_i$有一个正的权值$a_i$,每个人$stu_i$有一个负的权值$b_i$.每项活动能够被完成当且仅当该项活动所需的所有人到场.如何选择活动使最终权值总和最大?即对于给定的有向无环图,求出最大权闭合子图的权值. 结论:最大权闭合子图的权值等于所有正权点之和减去最小割. 思路:引理:1.最小割一定是简单割:2.简单割一定和一个闭合子图对应.即最小割一定对应一个闭合子图,且就是最大权闭合子图.证明(摘自HihoCoder):首先有割的容量C(S,

[BZOJ 1497][NOI 2006]最大获利(最大权闭合子图)

题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1497 分析: 这是在有向图中的问题,且边依赖于点,有向图中存在点.边之间的依赖关系可以考虑最大权闭合子图 假设a与b之间有权值为c的边(根据题意是双向边) 那么我们可以建一个新节点,点的权值为c,并指向a点和b点(单向),同时断掉原本a,b之间的双向边,a,b的点的权值是它们的花费(负的) 那么对于原问题就转化成了求最大权闭合子图的问题了 ——————————————————————

刷题总结——太空飞行计划(最大权闭合子图用最大流解决)

题目: 题目描述 W 教授正在为国家航天中心计划一系列的太空飞行.每次太空飞行可进行一系列商业性实验而获取利润.现已确定了一个可供选择的实验集合 E={E1,E2,-,Em},和进行这些实验需要使用的全部仪器的集合 I={I1, I2,-In}. 实验 Ej 需要用到的仪器是 I 的子集 Rj∈I.配置仪器 Ik 的费用为 Ck 美元.实验 Ej 的赞助商已同意为该实验结果支付 Pj 美元.W 教授的任务是找出一个有效算法, 确定在一次太空飞行中要进行哪些实验并因此而配置哪些仪器才能使太空飞行的

FZU - 2295 Human life (最大权闭合子图)

题目链接 FZU - 2295 Human life 题目分析 题意:你在玩一个游戏,在其中你可以通过学习一些技能,但是学习某些技能之前,可能还要学习一些其他的技能,并且学习任何技能都有一定的花费: 而我们可以通过掌握某些工作以获取报酬,为了掌握这一工作,我们必须学会特定的技能. 不过有些工作彼此之间是冲突的,简单来说:如果你掌握了工作A,那么将无法掌握工作B 思路: 由于技能之间也存在依赖关系,但实际上如果要获取某一工作的报酬,那么必须选择这个工作的前置技能以及前置技能的前置技能, 那么显然,

hiho一下 第119周 #1398 : 网络流五&#183;最大权闭合子图 【最小割-最大流--Ford-Fulkerson 与 Dinic 算法】

#1398 : 网络流五·最大权闭合子图 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 周末,小Hi和小Ho所在的班级决定举行一些班级建设活动. 根据周内的调查结果,小Hi和小Ho一共列出了N项不同的活动(编号1..N),第i项活动能够产生a[i]的活跃值. 班级一共有M名学生(编号1..M),邀请编号为i的同学来参加班级建设活动需要消耗b[i]的活跃值. 每项活动都需要某些学生在场才能够进行,若其中有任意一个学生没有被邀请,这项活动就没有办法进行. 班级建设的活

bzoj 1497 [NOI2006]最大获利【最大权闭合子图+最小割】

不要被5s时限和50000点数吓倒!大胆网络流!我一个5w级别的dinic只跑了1s+! 看起来没有最大权闭合子图的特征--限制,实际上还是有的. 我们需要把中转站看成负权点,把p看成点权,把客户看成正权点,把c看成点权,然后把中转站点a.b作为客户点的依赖点 s点向所有正权点连边,流量为点权:所有负权点向t连边,流量为负点权(即正数!) 对于所有有依赖关系的点,由客户点向中转站点连边,流量为inf,也就是最大权闭合子图中的向其依赖点连边 连边的意义详见:http://www.cnblogs.c