【网络流24题】太空飞行计划问题

（网络流24题大多需要spj，所以需要一个有spj的oj，本系列代码均在www.oj.swust.edu.cn测试通过）
最大权闭合子图裸题，先将所有收益加起来，源点向每个方案连接一条收益的流，每个方案向对应需要的仪器连一条INF，每个仪器向汇点连一条花费的流。
输出方案只需要在最后的阻塞流中看还剩哪几个点就好了

#include<cstdio>
#include<cstdlib>
#include<iomanip>
#include<iostream>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
#define INF 100000000
struct bian
{
    int l,r,f,lei;
}a[1000000];
int fir[1000000];
int nex[1000000];
int d[1000000];
int tot=1;
void add_edge(int l,int r,int f,int lei=0)
{
    a[++tot].l=l;
    a[tot].r=r;
    a[tot].f=f;
    a[tot].lei=lei;
    nex[tot]=fir[l];
    fir[l]=tot;
}
int s=0,t=999999;
bool bfs()
{
    static int dui[1000000];
    int top=1,my_final=2;
    memset(d,-1,sizeof(d));
    d[s]=0;
    dui[top]=s;
    while(top<my_final)
    {
        int u=dui[top];
        for(int o=fir[u];o!=0;o=nex[o])
        {
            if(d[a[o].r]!=-1 || !a[o].f) continue;
            dui[my_final++]=a[o].r;
            d[a[o].r]=d[u]+1;
            if(a[o].r==t) return true;
        }
        top++;
    }
    return false;
}
int dinic(int u,int flow)
{
    if(u==t) return flow;
    int left=flow;
    for(int o=fir[u];o!=0&&left;o=nex[o])
    {
        if(d[a[o].r]==d[u]+1 && a[o].f)
        {
            int temp=dinic(a[o].r,min(left,a[o].f));
            if(temp==0) d[a[o].r]=-1;
            left-=temp;
            a[o].f-=temp;
            a[o^1].f+=temp;
        }
    }
    return flow-left;
}
int main()
{
    int m,n;
    scanf("%d%d",&m,&n);
    int ans=0;
    for(int i=1;i<=m;i++)
    {
        char c;
        while(c<'0' || c>'9') scanf("%c",&c);
        int x=0;
        while(c<='9' && c>='0') x=x*10+c-'0',scanf("%c",&c);
        add_edge(s,i,x,1);
        add_edge(i,s,0);
        ans+=x;
        while(c!='\n' && c!='\r')
        {
            int y=0;
            while(c<'0' || c>'9')
            {
                if(c=='\n' || c=='\r') break;
                scanf("%c",&c);
            }
            if(c=='\n' || c=='\r') break;
            while(c<='9' && c>='0') y=y*10+c-'0',scanf("%c",&c);
            add_edge(i,m+y,INF);
            add_edge(m+y,i,0);
        }
    }
    for(int i=1;i<=n;i++)
    {
        int x;
        scanf("%d",&x);
        add_edge(m+i,t,x,2);
        add_edge(t,m+i,0);
    }
    while(bfs())
    {
        int t=dinic(s,INF);
        ans-=t;
    }
    for (int i=1;i<=m;i++)
        if (d[i]!=-1)
            printf("%d ",i);
    putchar('\n');
    for (int i=m+1;i<=m+n;i++)
        if (d[i]!=-1)
            printf("%d ",i-m);
    cout<<endl;
    cout<<ans<<endl;
    return 0;
}