最大获利（最大闭合图做法）：及其概念辨析和证明

最大获利

新的技术正冲击着手机通讯市场，对于各大运营商来说，这既是机遇，更是挑战。

THU 集团旗下的 CS&T 通讯公司在新一代通讯技术血战的前夜，需要做太多的准备工作，仅就站址选择一项，就需要完成前期市场研究、站址勘测、最优化等项目。

在前期市场调查和站址勘测之后，公司得到了一共 $N$ 个可以作为通讯信号中转站的地址，而由于这些地址的地理位置差异，在不同的地方建造通讯中转站需 要投入的成本也是不一样的，所幸在前期调查之后这些都是已知数据：

建立第 $i$ 个通讯中转站需要的成本为 $P_i(1\le i\le N)$。

另外公司调查得出了所有期望中的用户群，一共 $M$ 个。

关于第 $i$ 个用户群的信息概括为 $A_i,B_i$ 和 $C_i$：这些用户会使用中转站 $A_i$ 和中转站 $B_i$ 进行通讯，公司可以获益 $C_i。（1\le i\le M,1\le A_i,B_i\le N）$

THU 集团的 CS&T 公司可以有选择的建立一些中转站（投入成本），为一些用户提供服务并获得收益（获益之和）。

那么如何选择最终建立的中转站才能让公司的净获利最大呢？（净获利 = 获益之和 – 投入成本之和）

输入格式

第一行有两个正整数 $N$ 和 $M$。

第二行中有 $N$ 个整数描述每一个通讯中转站的建立成本，依次为 $P_1,P_2,\dots ,P_N$ 。

以下 $M$ 行，第 $(i+2)$ 行的三个数 $A_i,B_i$ 和 $C_i$ 描述第 $i$ 个用户群的信息。

所有变量的含义可以参见题目描述。

输出格式

输出一个整数，表示公司可以得到的最大净获利。

数据范围

$1 \le N \le 5000$, $1 \le M \le 50000$, $0 \le C_i,P_i \le 100$

输入样例：

5 5
1 2 3 4 5
1 2 3
2 3 4
1 3 3
1 4 2
4 5 3

输出样例：

4

样例解释

选择建立 $1、2、3$ 号中转站，则需要投入成本 $6$，获利为 $10$，因此得到最大收益 $4$。

概念解析及本道题思路

• 概念就是点的集合，集合里面的点与点的边不能指向集合外面，这样就叫作闭合图（本质上是点的集合）。

像这个就是闭合图（集合外面的指向里面的可以）。

那么此时我们可以建立如上图所示的流网络，我们是把正的权值跟源点连边，负的权值跟汇点连边，中间的边设为无穷。为什么呢？因为原图中间的限制我们不变，就设为无穷。

简单割概念：简单割就是割在如上图橙色部分，不在中间部分的成为简单割。

我们就得知道最小割是不是简单割，发现是，为什么呢？因为最小割=最大流，我们知道可行流是有具体数值的，因此最小割就等于简单割。

一定要注意简单割我们是按定义来的，不要形象理解，要不然会头脑风暴。

又因为这个关系，所以我们可以求所有割的最小割，此时就是简单割了。

• 现在就证明简单割和闭合图的关系：

我要证明原问题能构造出一种割是简单割，我们可以按上图所示建立源点集合和汇点集合，其中 $V^{\prime }$ 指的是原问题的点集，我们发现原问题是闭合子图，因此 $V^{\prime }$ 无论怎么流都只能回到 $V^{\prime }$ ，因此它不可能从 $S$ 走到 $T$，那么同理 $T$ 也不能走到 $S$。

那么此时我们就能证明闭合子图等价于简单割。

那么我们该怎么求简单割呢？此时我们设简单割为 $[S,T]$，闭合子图为 $V1,V2$。此时满足： V 1 = S − { s } V_1=S-\{s\} V1​=S−{s}，此时我们就能画出：

因为 $V_1,V_2$ 为内部点（这边就直接把前者取成权值为正的点，后者为负数），不能相互连接， $s,t$ 之间也不能相互连接，因此：

此时我们就是求最小割即可。

注意1：简单割是针对最大权闭合图而言的，别的没说。

拓展：对于证明原图跟流网络的关系，如果在理论上能相互推出，那么数量上一般能相互对应。（就像上述证明，最终就是求出割与简单割的对应关系）

本道题的题意：

代码

//也就是说，每个用户必须跟两个基站绑定
//把用户标正，中转站标负
//而且所有用户只能在给定的基站内相互通信，因此这就是一个闭合子图

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>

using namespace std;

const int N = 55010,M = (50000*3+N)*2+10,INT = 1e8;

int e[M],ne[M],f[M],h[N],idx;
int cur[N],d[N];
int n,m,S,T;
int tot;

void add(int a,int b,int c){
    e[idx]=b,f[idx]=c,ne[idx]=h[a],h[a]=idx++;
    e[idx]=a,f[idx]=0,ne[idx]=h[b],h[b]=idx++;
}

bool bfs(){
    queue<int>q;
    q.push(S);

    memset(d,-1,sizeof d);

    d[S]=0,cur[S]=h[S];

    while(q.size()){
        int t=q.front();
        q.pop();

        for(int i=h[t];~i;i=ne[i]){
            int ver=e[i];

            if(d[ver]==-1&&f[i]){
                d[ver]=d[t]+1;
                cur[ver]=h[ver];
                if(ver==T)return true;
                q.push(ver);
            }
        }
    }
    return false;
}

int find(int u,int lim){
    if(u==T)return lim;

    int flow=0;

    for(int i=cur[u];~i&&flow<lim;i=ne[i]){
        int ver=e[i];

        cur[u]=i;

        if(d[ver]==d[u]+1&&f[i]){
            int t=find(ver,min(f[i],lim-flow));
            if(!t)d[ver]=-1;
            f[i]-=t,f[i^1]+=t,flow+=t;
        }
    }
    return flow;
}

int dinic(){
    int r=0,flow;
    while(bfs())while(flow=find(S,INT))r+=flow;

    return r;
}

int main(){
    cin>>n>>m;
    memset(h,-1,sizeof h);

    S=0,T=n+m+1;

    for(int i=1;i<=n;i++){//这个是消耗品
        int x;
        cin>>x;
        add(m+i,T,x);
    }

    for(int i=1;i<=m;i++){
        int a,b,c;
        cin>>a>>b>>c;
        add(S,i,c);
        add(i,a+m,INT);
        add(i,m+b,INT);
        tot+=c;
    }

    cout<<tot-dinic();

    return 0;
}