cogs 727. [网络流24题] 太空飞行计划

1. [网络流24题] 太空飞行计划
   ★★☆ 输入文件：shuttle.in 输出文件：shuttle.out 简单对比
   时间限制：1 s 内存限制：128 MB
   【问题描述】
   W 教授正在为国家航天中心计划一系列的太空飞行。每次太空飞行可进行一系列商业性实验而获取利润。现已确定了一个可供选择的实验集合E={E1,E2,…,Em}，和进行这些实验需要使用的全部仪器的集合I={ I1, I2,…,In }。实验Ej 需要用到的仪器是I的子集Rj∈I。配置仪器Ik 的费用为ck 美元。实验Ej 的赞助商已同意为该实验结果支付pj 美元。W教授的任务是找出一个有效算法，确定在一次太空飞行中要进行哪些实验并因此而配置哪些仪器才能使太空飞行的净收益最大。这里净收益是指进行实验所获得的全部收入与配置仪器的全部费用的差额。
   【编程任务】
   对于给定的实验和仪器配置情况，编程找出净收益最大的试验计划。
   【数据输入】
   第1行有2个正整数m和n(m,n <= 100)。m是实验数，n是仪器数。接下来的m行，每行是一个实验的有关数据。第一个数赞助商同意支付该实验的费用；接着是该实验需要用到的若干仪器的编号。最后一行的n个数是配置每个仪器的费用。
   【结果输出】
   第1行是实验编号；第2行是仪器编号；最后一行是净收益。
   【输入文件示例】shuttle.in
   2 3
   10 1 2
   25 2 3
   5 6 7

【输出文件示例】shuttle.out
1 2
1 2 3
17

---

【分析】
裸最大权闭合子图，有套路，有证明。
http://blog.csdn.net/qq_33229466/article/details/51018916

---

【代码】

//cogs [网络流24题] 太空飞行计划
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#define inf 1e9+7
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
queue <int> q;
const int mxn=100005;
int n,m,cnt,ans,t=10000;
int head[mxn],dis[mxn],get[mxn];
struct node {int to,next,flow;} f[mxn<<2]; 
inline void add(int u,int v,int flow)
{
    f[++cnt].to=v,f[cnt].next=head[u],f[cnt].flow=flow,head[u]=cnt;
    f[++cnt].to=u,f[cnt].next=head[v],f[cnt].flow=0,head[v]=cnt;
}
inline bool bfs()
{
    int i,j,u,v,flow;
    memset(dis,-1,sizeof dis);
    q.push(0);
    dis[0]=0;
    while(!q.empty())
    {
        u=q.front();
        q.pop();
        for(i=head[u];i;i=f[i].next)
        {
            v=f[i].to,flow=f[i].flow;
            if(dis[v]==-1 && flow>0)
              dis[v]=dis[u]+1,q.push(v);
        }
    }
    if(dis[t]>0) return 1;
    return 0;
}
inline int find(int u,int low)
{
    int v,i,j,a=0,sum=0,flow;
    if(u==t) return low;
    for(i=head[u];i;i=f[i].next)
    {
        v=f[i].to,flow=f[i].flow;
        if(flow>0&& dis[v]==dis[u]+1 && (a=find(v,min(low-sum,flow))))
        {
            f[i].flow-=a;
            if(i&1) f[i+1].flow+=a;
            else f[i-1].flow+=a;
            sum+=a;
        }
    }
    if(!sum) dis[u]=-1;
    return sum;
}
int main()
{
    freopen("shuttle.in","r",stdin);
    freopen("shuttle.out","w",stdout);
    char ch;
    int i,j,u,v,x,tans,flow,now,sum=0;
    scanf("%d%d",&m,&n);  //m个实验,n个仪器 
    fo(i,1,m)   //辣鸡读入 
    {
        scanf("%d",&x);
        sum+=x;get[i]=x;
        add(0,i,x);
        ch=getchar();
        while(ch!='\r' && ch!='\n')
        {
            now=0;
            while(ch>='0'&&ch<='9')
            {
                x=ch-'0';
                now=now*10+x;
                ch=getchar();
            }
            add(i,now+300,inf);
            if(ch=='\r' || ch=='\n') break;
            ch=getchar();
        }
    }
    fo(i,1,n) 
    {
        scanf("%d",&x);
        add(i+300,t,x);
    }
    while(bfs())
      while(tans=find(0,0x7f7f))
        ans+=tans;
    bfs();
    fo(u,1,m) if(dis[u]>0) printf("%d ",u);printf("\n");
    fo(u,301,n+300) if(dis[u]>0) printf("%d ",u-300);printf("\n");
    printf("%d\n",sum-ans);
    return 0;
}
//2 3
//1 1 2
//25 2 3
//5 6 7