网络流专题 太空飞行计划（最大权闭合子图）

 

题目描述

W 教授正在为国家航天中心计划一系列的太空飞行。每次太空飞行可进行一系列商业性实验而获取利润。现已确定了一个可供选择的实验集合E={E1，E2，…，Em}，和进行这些实验需要使用的全部仪器的集合I={I1，I2，…In}。实验Ej需要用到的仪器是I的子集RjÍI。配置仪器Ik的费用为ck美元。实验Ej的赞助商已同意为该实验结果支付pj美元。W教授的任务是找出一个有效算法，确定在一次太空飞行中要进行哪些实验并因此而配置哪些仪器才能使太空飞行的净收益最大。这里净收益是指进行实验所获得的全部收入与配置仪器的全部费用的差额。

对于给定的实验和仪器配置情况，编程找出净收益最大的试验计划。

输入输出格式

输入格式：

 

第1行有2 个正整数m和n。m是实验数，n是仪器数。接下来的m 行，每行是一个实验的有关数据。第一个数赞助商同意支付该实验的费用；接着是该实验需要用到的若干仪器的编号。最后一行的n个数是配置每个仪器的费用。

 

输出格式：

 

第1 行是实验编号；第2行是仪器编号；最后一行是净收益。

 

输入输出样例

输入样例#1： 复制

2 3
10 1 2
25 2 3
5 6 7

输出样例#1： 复制

1 2
1 2 3
17

说明

感谢@FlierKing 提供spj

n,m<=50

这道题数据是在windows生成的，输入数据中所有的换行都是'\r\n'而不是'\n'
读入某实验需要用到的仪器编号的时候，可以这么读入。(感谢@zhouyonglong的提供)

题目大意：给出一些实验，做这些实验有一个收益，同时需要买一些仪器每个仪器有一个花费。求最大收益并且输出方案。

解题思路：每个实验如果选的话就要买对应的仪器，一起可以供另一些实验使用。仪器与实验存在依赖关系。

转化为网络流模型：如果把每个实验与仪器之间连一条边的话，如果选则该实验则该实验所连的边对应的点也要选，这明显符合闭合子图的定义。

闭合子图点这里：提示：最大权闭合子图

关于建图：建立一个源点和汇点，每个实验与原点连一条边，权值为实验的收益，实验与仪器之间连一条容量为inf的边，仪器与

汇点连一条容量为花费的边。然后这个题就是让求一个最大权闭合子图。

一个定理：最大权闭合子图的权值等于所有正权点之和减去最小割。证明在上边的链接里。

根据我们建的图可知，最小割把原来的二分图划分成两个集合S,T，割边的时候一定不会割掉实验与仪器之间的边。

总的来说就是：最小割也对应了一个闭合子图，这个子图就是S集合去点源点，并且这个子图的权值最大。

这个题读入有点坑,我借鉴了一下别人是怎么读入的。

这个我用的是EK求最大流并且重新bfs了一下源点能到达的点（这些就是S），如果用dinic去做的话，最后一次分层有标号的点就是S集合。

#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define sca(x) scanf("%d",&x)
#define per(i,j,k) for(int i=j;i>=k;i--)
#define inf 0x3f3f3f3f
#define LL long long
#define N 40005
#define inf 0x3f3f3f3f
const LL mod = 998244353;
int flag;
int n,m;

struct edg
{
    int to,w,nt;
}g[N];

struct node
{
    int id,v;
}pre[N];

int vis[N],head[N];

bool bfs(int s,int t)
{
    memset(vis,false,sizeof(vis));
    queue<int>q;
    vis[s]=0;
    q.push(s);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=head[u];i!=-1;i=g[i].nt)
        {
            int to=g[i].to;
            if(!vis[to]&&g[i].w>0)
            {
                vis[to]=true;
                pre[to].v=u;
                pre[to].id=i;
                if(to==t)return true;
                q.push(to);
            }
        }
    }
    return false;
}

int EK(int s,int t)
{
    int ans=0;
    while(bfs(s,t))
    {
        int mini=inf;
        for(int i=t;i!=s;i=pre[i].v)
        {
            mini=min(mini,g[pre[i].id].w);
        }
        for(int i=t;i!=s;i=pre[i].v)
        {
            g[pre[i].id].w-=mini;
            g[pre[i].id^1].w+=mini;
        }
        ans+=mini;
    }
    return ans;
}

int tot;
void addedg(int u,int v,int w)
{
    g[tot].to=v;
    g[tot].w=w;
    g[tot].nt=head[u];
    head[u]=tot++;
}

void get_set(int s,int t)
{
    int Set[205];
    memset(vis,false,sizeof(vis));
    memset(Set,0,sizeof(Set));
    queue<int>q;
    vis[s]=true;
    q.push(s);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        Set[u]=1;
        for(int i=head[u];i!=-1;i=g[i].nt)
        {
            int to=g[i].to;
            if(!vis[to]&&g[i].w>0)
            {
                vis[to]=1;
                q.push(to);
            }
        }
    }
    for(int i=1;i<=n;i++)
    {
        if(Set[i])printf("%d ",i);
    }
    printf("\n");
    for(int i=n+1;i<=110;i++)
    {
        if(Set[i]&&i!=t)printf("%d ",i-n);
    }
    printf("\n");
}

int main()
{
    scanf("%d%d",&n,&m);
    int sum=0;
    memset(head,-1,sizeof(head));
    rep(i,1,n)
    {
       int tmp;
       scanf("%d",&tmp),sum+=tmp;
       addedg(0,i,tmp);
       addedg(i,0,0);
       char c;
       while(1)
       {
           scanf("%d%c",&tmp,&c);
           addedg(i,n+tmp,inf);
           addedg(n+tmp,i,0);
           if(c=='\n'||c=='\r')break;
       }
    }
    rep(i,1,m)
    {
       int tmp;
       sca(tmp);
       addedg(n+i,m+n+1,tmp);
       addedg(n+m+1,n+i,0);
    }
    int ans=EK(0,n+m+1);
    get_set(0,n+m+1);
    printf("%d\n",sum-ans);
}