太空飞行计划

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

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

##解题思路：
最大权闭合子图：源点连正点权（即实验收益），容量为正；负点权（器材花费）连汇点，容量为负的负（花费较总收益来说原本就是负的，连得时候连正的就行），有关系的实验和器材之间连边容量为正无穷，答案为全部实验收益 - 上面图跑出来的最小割
解释：正点之间与负点连边为INF最小割跑出来，要么割正边，意义就是舍弃这个实验，要么割负的负边，意义为买器材花钱，无论怎样，要是使得放弃实验钱与花的器材钱最少（也就是最小割），最终答案就是最优的

##源程序：

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
const int inf=1e9;
int n,m,p[51],c[51],x,p1[51][51],S=0,T,sum=0,ans,mi=inf,yu=0;
int to[6001],w[6001],w1[6001],next[6001],first[201],e1=0;
int d[201],f[6001],cur[201],e[201],now,check[51];
int read()
{
    char ch=getchar();yu=0;
    while((ch<'0')||(ch>'9')) {if(ch=='\n') return 0;ch=getchar();}
    while((ch>='0')&&(ch<='9')) yu=yu*10+ch-'0',ch=getchar();
    if(ch=='\n') return 0;
    return yu;
}
void insert(int o,int p,int q)
{
    to[e1]=p;w[e1]=w1[e1]=q;next[e1]=first[o];first[o]=e1++;
    to[e1]=o;w[e1]=w1[e1]=0;next[e1]=first[p];first[p]=e1++;
}
int sap()
{
    int yus=0;
    for(int i=0;i<=T;i++) cur[i]=first[i];
    memset(e,-1,sizeof(e));
    memset(d,0,sizeof(d));
    memset(f,0,sizeof(f));
    e[S]=now=S;f[0]=T;
    while(d[S]<T)
    {
        bool bo=false;
        for(int i=cur[now];i!=-1;i=next[i])
            if((d[now]==(d[to[i]]+1))&&(w1[i]>0))
            {
                cur[now]=i;
                e[to[i]]=now;
                now=to[i];
                bo=true;
                mi=min(mi,w1[i]);
                if(now==T)
                {
                    yus+=mi;
                    while(now!=S)
                    {
                        w1[cur[e[now]]]-=mi;
                        w1[cur[e[now]]^1]+=mi;
                        now=e[now];
                    }
                    mi=inf;
                }
                break;
            }
        if(bo) continue;int mmi=T;
        for(int i=first[now];i!=-1;i=next[i])
            if((w1[i]>0)&&(mmi>d[to[i]])) mmi=d[to[i]],cur[now]=i;
        if(!(--f[d[now]])) break;
        f[d[now]=(mmi+1)]++;
        now=e[now];
    }
    return yus;
}
int main()
{
    memset(first,-1,sizeof(first));
    scanf("%d%d",&m,&n);
    for(int i=1;i<=m;i++)
    {
        scanf("%d",&p[i]);sum+=p[i];
        while(true) 
        {
            x=read();
            p1[i][++p1[i][0]]=yu;
            if(!x)break;
        }
        S+=p1[i][0];
    }
    for(int i=1;i<=n;i++) scanf("%d",&c[i]);
    S+=n+m+1;T=S+1;
    for(int i=1;i<=m;i++) 
    {
        insert(S,i,p[i]);
        for(int j=1;j<=p1[i][0];j++) insert(i,m+p1[i][j],inf);
    }
    for(int i=1;i<=n;i++) insert(i+m,T,c[i]);
    ans=sap();
    memset(check,0,sizeof(check));
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<=e1;j++) w1[j]=w[j];
        w1[first[m+i]]=0;
        if(ans-sap()==w[first[m+i]]) check[i]=true;
    }
    for(int i=1;i<=m;i++)
    {
        bool b1=true;
        for(int j=1;j<=p1[i][0];j++)
            if(!check[p1[i][j]]) {b1=false;break;}
        if(b1) cout<<i<<' ';
    }
    cout<<endl;
    for(int i=1;i<=n;i++)
        if(check[i]) cout<<i<<' ';
    cout<<endl;
    cout<<sum-ans<<endl;
    return 0;
}