OpenJudge 3765(最大权闭合图，最小割

最新推荐文章于 2019-01-07 19:11:37 发布

aduqw1129

最新推荐文章于 2019-01-07 19:11:37 发布

阅读量73

点赞数

原文链接：http://www.cnblogs.com/Cw-trip/p/4733795.html

版权

探讨了解决最大权闭合图问题的一种算法实现，该问题涉及在公司裁员场景中，如何通过解雇特定员工及其下属以达到最大收益。文章详细介绍了使用最大流最小割算法来解决这一问题的方法，并提供了完整的C++代码实现。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目大意：需要解雇若干人，解雇每个人有一个损失或收益，解雇一个人就必须结果他的所有下属，问最大收益。

思路：裸最大权闭合图问题，对于权值为q的点，负权向汇点连-p的边，正权源点向其连p的边，上下级关系上级向下级连oo的边，然后求最大流最小割，最大收益即为总正点权减最小割，需要裁员的人即为残余图中从源点可达的点。

参考：http://blog.youkuaiyun.com/scorpiocj/article/details/6085637

#include <bits/stdc++.h>
#define pb push_back
#define se second
#define fs first
#define sq(x) (x)*(x)
#define eps 0.000000001
#define LB lower_bound
#define IINF (1<<29)
#define LINF (1ll<<59)
using namespace std;
typedef long long ll;
typedef pair<ll,ll> P;
const int maxv=5010;
struct EDGE{
    int to,rev;
    ll cap;
    EDGE(int t,ll c,int r):to(t),cap(c),rev(r){}
};
vector<EDGE> G[maxv];
void addedge(int from,int to,ll cap){
    G[from].pb(EDGE(to,cap,G[to].size()));
    G[to].pb(EDGE(from,0,G[from].size()-1));
}
int level[maxv];
queue<int> Q;
void bfs(int s){
    memset(level,-1,sizeof level);
    level[s]=0;
    Q.push(s);
    while(!Q.empty()){
        int v=Q.front();Q.pop();
        for(int i=0;i<G[v].size();i++){
            EDGE &e=G[v][i];
            if(level[e.to]==-1&&e.cap>0){
                level[e.to]=level[v]+1;
                Q.push(e.to);
            }
        }
    }
}
int iter[maxv];
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];
        ll d;
        if(e.cap>0&&level[e.to]>level[v]){
            d=dfs(e.to,t,min(e.cap,f));
            if(d>0){
                e.cap-=d;
                G[e.to][e.rev].cap+=d;
                return d;
            }
        }
    }
    return 0;
}
ll dinic(int s,int t){
    ll flow=0;
    for(;;){
        bfs(s);
        if(level[t]==-1){
            return flow;
        }
        ll f;
        memset(iter,0,sizeof iter);
        while((f=dfs(s,t,LINF))>0) flow+=f;
    }
}
int cc=0;
bool vis[maxv];
void cont(int v){
    cc++;
    vis[v]=1;
    for(int i=0;i<G[v].size();i++){
        EDGE &e=G[v][i];
        if(!vis[e.to]&&e.cap>0){
            cont(e.to);
        }
    }
}
const int s=5008,t=5009;
int n,m;
int pro[maxv];
int main(){
    freopen("/home/files/CppFiles/in","r",stdin);
    cin>>n>>m;
    ll sum=0;
    for(int i=1;i<=n;i++){
        ll p;
        scanf("%lld",&p);
        pro[i]=p;
        if(p>0){
            addedge(s,i,p);
            sum+=p;
        }
        else 
            addedge(i,t,-p);
    }
    for(int i=0;i<m;i++){
        ll x,y;
        scanf("%lld%lld",&x,&y);
        addedge(x,y,LINF);
    }
    ll ans=sum-dinic(s,t);
    cont(s);
    cout<<cc-1<<" "<<ans<<endl;
}