[最大权闭合子图]太空飞行计划网络流24题（2/24）

最新推荐文章于 2021-03-26 23:05:47 发布

Mr_Alice

最新推荐文章于 2021-03-26 23:05:47 发布

阅读量176

点赞数

分类专栏：网络流

本文链接：https://blog.csdn.net/Mr_Alice/article/details/102593227

版权

网络流专栏收录该内容

6 篇文章

订阅专栏

本文深入探讨了最大权闭合子图问题的解决方法，通过详细解释算法流程和代码实现，帮助读者理解如何在给定图中找到与源点连通且权重最大的子图，特别适用于网络流和图论的学习。

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

无序依赖。

注意最后输出方案时，即要输出最后的最大权闭合子图。

最大权闭合子图即最后还与S连通的点。

最后与S连通的点d一定不等于0（or-1，看写法了

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define INF (1e9)
const ll maxn=115,maxm=30005;
ll m,n,cnt,S,T;
ll w[maxn],cost[maxn],head[maxn];
struct node{
    ll to,ne,w;
}e[maxm];
void add(ll u,ll v,ll w){
    e[cnt]=(node){v,head[u],w};head[u]=cnt++;
    e[cnt]=(node){u,head[v],0};head[v]=cnt++;
}
queue<ll>q;
ll d[maxn];
ll bfs(){
    memset(d,0,sizeof(d));
    while(!q.empty()) q.pop();
    q.push(S);
    d[S]=1;
    while(!q.empty()){
        ll u=q.front();q.pop();
        if(u==T) break;//优化
        for(ll i=head[u];~i;i=e[i].ne){
            ll v=e[i].to;
            if(!e[i].w) continue;
            if(!d[v]) d[v]=d[u]+1,q.push(v);
        }
    }
    return d[T];
}

ll dfs(ll u,ll remain){
    if(u==T) return remain;
    ll use=0;
    for(ll &i=head[u];~i;i=e[i].ne){
        ll v=e[i].to;
        if(e[i].w&&remain&&d[v]==d[u]+1){
            ll flow=dfs(v,min(remain,e[i].w));
            use+=flow;remain-=flow;e[i].w-=flow;e[i^1].w+=flow;
        }
    }
    if(!use) d[u]=0;//优化
    return use;
}
ll rechead[maxn];
ll dinic(){
    ll ans=0;
    memcpy(rechead,head,sizeof(head));
    while(bfs()){
        while(1){
            ll tmp=dfs(S,INF);
            if(!tmp) break;
            ans+=tmp;
        }
        memcpy(head,rechead,sizeof(head));
    }
    return ans;
}
int main(){
    scanf("%lld%lld",&m,&n);
    memset(head,-1,sizeof(head));
    S=m+n+1;T=S+1;
    ll x,sum=0;
    for(ll i=1;i<=m;i++){
        scanf("%lld",&w[i]);
        sum+=w[i];
        add(S,i,w[i]);
        while(getchar()==' '){
            scanf("%lld",&x);
            add(i,x+m,INF);
        }
    }
    for(ll i=1;i<=n;i++){
        scanf("%lld",&x);
        add(i+m,T,x);
    }
    sum-=dinic();
    memcpy(head,rechead,sizeof(head));
    //求方案。选择的仪器必在最大权闭合子图内，故一定与S连通，所以d一定不是0
     for(ll i=1;i<=m+n;i++){
        if(d[i]){
            printf("%lld ",i>m?i-m:i);
        }
        if(i==m||i==n+m) printf("\n");
    }
    printf("%lld\n",sum);
    return 0;
}