洛谷 P2762 太空飞行计划问题 【最小割--最大权闭合子图】

最大权闭合子图

定义：

如果对于一个点集合，其中任何一个点都不能到达此集合以外的点，这就叫做闭合子图。每个点都有一个权值，那么最大权闭合子图就是权值最大的那个闭合子图。
(或者说对于一个点集，这个点集中所有点的出边所指向的点都在此点集中)

求解

超级源点向每个权值为正的点连边，容量为该点权值
每个点权为负的点向超级汇点连边，容量为该点权值相反数
原图中的边，容量为inf
然后跑最小割
最后用正点权的总和-最大流即为最大权闭合子图的权值

证明

蒟蒻自己的理解：
首先按这样的方式建图最小割一定是简单割(只隔与s,t关联的边)
割与s关联的边，表示不选这个点，要付出这个点正权值的代价
割与T关联的边，表示选这个点，要付出这块点的负权值的代价
所以最大的闭合图权值就是总正权值-最小割

形式化的证明，参见——胡伯涛《最小割模型在信息学竞赛中的应用》

---

时空限制 1000ms / 128MB

题目描述

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行是仪器编号；最后一行是净收益。

---

题目分析

超源向每个实验连边，容量为其能得的支付费用
每个实验向其对应的器材连边，容量为inf
每个器材向超源连边，容量为其费用

在最小割求出后在建立一次层次图
有层次的就是选择的实验与仪器
(卡死在输入系列+_+)

---

#include<iostream>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int inf=1e9;
const int maxn=500;
int n,m,s=0,t;
struct node{int v,f,nxt;}E[maxn*200];
int head[maxn],tot=1;
int lev[maxn],sum,maxf;

void add(int u,int v,int f)
{
    E[++tot].nxt=head[u];
    E[tot].v=v; E[tot].f=f;
    head[u]=tot;
}

int bfs()
{
    memset(lev,-1,sizeof(lev)); lev[s]=0;
    queue<int> q; q.push(s);
    while(!q.empty())
    {
        int u=q.front(); q.pop();
        for(int i=head[u];i;i=E[i].nxt)
        {
            int v=E[i].v;
            if(lev[v]==-1&&E[i].f)
            {
                lev[v]=lev[u]+1;
                if(v==t) return 1;
                q.push(v);
            }
        }
    }
    return 0;
}

int dfs(int u,int cap)
{
    if(u==t) return cap;
    int flow=cap;
    for(int i=head[u];i;i=E[i].nxt)
    {
        int v=E[i].v;
        if(lev[v]==lev[u]+1&&E[i].f&&flow)
        {
            int f=dfs(v,min(E[i].f,flow));
            flow-=f;
            E[i].f-=f; E[i^1].f+=f; 
        }
    }
    return cap-flow;
}

int main()
{
    n=read();m=read();t=n+m+1;
    for(int i=1;i<=n;++i)
    {
    	int k=read(); sum+=k;
    	add(s,i,k); add(i,s,0);
    	
    	char ss=getchar();
        while(ss!='\n'&&ss!='\r')//这个输入自行体会
        {
            int x=0;
            while(ss>='0'&&ss<='9') x=x*10+ss-'0',ss=getchar();
            if(x) add(i,x+n,inf),add(x+n,i,0);
            if(ss!='\n'&&ss!='\r') ss=getchar();
        }
    } 
    for(int i=1;i<=m;++i)
    add(i+n,t,read()),add(t,i+n,0);
    
    while(bfs())
    maxf+=dfs(s,inf);
    
    bfs();
    for(int i=1;i<=n;++i)
    if(lev[i]!=-1) printf("%d ",i);
    printf("\n");
    
    for(int i=n+1;i<=n+m;++i)
    if(lev[i]!=-1) printf("%d ",i-n);
    printf("\n%d",sum-maxf);
    return 0;
}

---

领三倍经验
洛谷P4174 [NOI2006]最大获利
洛谷P3410 拍照