COGS727 [网络流24题] 太空飞行计划

【问题描述】

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

【输出文件示例】shuttle.out

1 2
1 2 3
17

图论网络流特别的读入技巧

题目本身就是个裸的最大权闭合子图，没什么好说的

但是这个输入数据很好玩啊，每行没有终止标记。

但是读入也没什么难的，把读入优化魔改了一下就可以直接用了。

交上去发现还是T了，下回来数据发现——数据文件前两行有两个换行符，然后才是m和n，于是17行那样的写法就直接炸了

多加了一个enter标记，换成22行和105行那样的写法就妥了。

中途因为忘了改文件操作就交，又炸了两次

本来以为能1A来着……获得技能[迷之自信]　@布吕歇尔自信姬

  1 /*by SilverN*/
  2 #include<algorithm>
  3 #include<iostream>
  4 #include<cstring>
  5 #include<cstdio>
  6 #include<cmath>
  7 #include<vector>
  8 #include<queue>
  9 using namespace std;
 10 const int INF=0x3f3f3f3f;
 11 const int mxn=100010;
 12 bool enter=0;
 13 int read(){
 14     int x=0,f=1;char ch=getchar();
 15     while(ch<'0' || ch>'9'){
 16         if(ch=='-')f=-1;
 17 //        if(ch=='
' || ch=='
')return -1;//并没有卵用 
 18         ch=getchar();
 19     }
 20     while(ch>='0' && ch<='9'){
 21         x=x*10+ch-'0';ch=getchar();
 22         if(ch=='
' || ch=='
')enter=1;
 23     }
 24     return x*f;
 25 }
 26 inline int min(int a,int b){return a<b?a:b;}
 27 struct edge{
 28     int v,nxt,f;
 29 }e[mxn<<1];
 30 int hd[mxn],mct=1;
 31 void add_edge(int u,int v,int f){
 32     e[++mct].v=v;e[mct].nxt=hd[u];e[mct].f=f;hd[u]=mct;return;
 33 }
 34 void insert(int u,int v,int c){
 35     add_edge(u,v,c);add_edge(v,u,0);return;
 36 }
 37 int S,T;
 38 int d[300];
 39 queue<int>q;
 40 bool BFS(){
 41     memset(d,0,sizeof d);
 42     q.push(S);d[S]=1;
 43     while(!q.empty()){
 44         int u=q.front();q.pop();
 45         for(int i=hd[u],v;i;i=e[i].nxt){
 46             v=e[i].v;
 47             if(e[i].f && !d[v]){
 48                 d[v]=d[u]+1;
 49                 
 50                 q.push(v);
 51             }
 52         }
 53     }
 54     return d[T];
 55 }
 56 int DFS(int u,int lim){
 57     if(u==T)return lim;
 58     int f=0,tmp;
 59     for(int i=hd[u],v;i;i=e[i].nxt){
 60         v=e[i].v;
 61         if(e[i].f && d[v]==d[u]+1 && (tmp=DFS(v,min(e[i].f,lim)))){
 62             e[i].f-=tmp;
 63             e[i^1].f+=tmp;
 64             lim-=tmp;
 65             f+=tmp;
 66             if(!lim)return f;
 67         }
 68     }
 69     d[u]=0;
 70     return f;
 71 }
 72 int Dinic(){
 73     int res=0;
 74     while(BFS())res+=DFS(S,INF);
 75     return res;
 76 }
 77 int n,m;
 78 int ans=0;
 79 void solve(){
 80     ans-=Dinic();
 81     for(int i=1;i<=m;i++)
 82         if(d[i])printf("%d ",i);
 83     puts("");
 84     for(int i=1;i<=n;i++){
 85         if(d[i+m])printf("%d ",i);
 86     }
 87     puts("");
 88     printf("%d
",ans);
 89     return;
 90 }
 91 int main(){
 92 //    freopen("in.txt","r",stdin);
 93     freopen("shuttle.in","r",stdin);
 94     freopen("shuttle.out","w",stdout);
 95     int i,j,x;
 96     m=read();n=read();
 97     S=0;T=m+n+1;
 98     for(i=1;i<=m;i++){
 99         x=read();
100         ans+=x;
101         insert(S,i,x);
102         enter=0;//
103         while(!enter){
104             x=read();
105             if(x==-1)break;
106             insert(i,x+m,INF);
107         }
108     }
109     for(i=1;i<=n;i++){
110         x=read();
111         insert(i+m,T,x);
112     }
113     solve();
114     return 0;
115 }