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

【问题描述】

W 教授正在为国家航天中心计划一系列的太空飞行。每次太空飞行可进行一系列商业性实验而获取利润。现已确定了一个可供选择的实验集合 E ={ E 1 , E 2 ,…, E m }，和进行这些实验需要使用的全部仪器的集合 I ={  I 1 ,  I 2 ,…, I n  }。实验 Ej  需要用到的仪器是 I 的子集 R j ∈ I 。配置仪器 I k  的费用为 c k  美元。实验 E j  的赞助商已同意为该实验结果支付 p j  美元。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

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

最小割~

从S向实验连边，边权为实验收益；实验向所需器材连边，边权正无穷，保证不会被割；器材向T连边，边权为所需费用。

加边的时候顺便记录ans表示所有实验的总收益。

然后跑最小割，割去S-实验的边表示不做这个实验，减去这个实验的收益；割去器材-T的边表示购买这个器材，多花费器材的价钱，用ans减去就可以了。

……

这道题最厉害的不是思路，而是输入输出……输入判断换行符，换行符还必须是'\r'！出题人肯定是故意的！而且这题因为读入的问题不能在本地测，本地测是RE的！

输出直接判断dis是否为0就可以了，不用再dfs一遍。

dis[u]==dis[k]+1！检查到吐血……

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;
#define inf 999999999

int n,m,x,fi[205],w[25001],ne[25001],v[25001],cnt,dis[205],ans,kk[205];
bool b[205];

void add(int u,int vv,int val)
{
	w[++cnt]=vv;ne[cnt]=fi[u];fi[u]=cnt;v[cnt]=val;
	w[++cnt]=u;ne[cnt]=fi[vv];fi[vv]=cnt;
}

bool bfs()
{
	memset(dis,0,sizeof(dis));
	queue<int> q;q.push(0);dis[0]=1;
	while(!q.empty())
	{
		int k=q.front();q.pop();
		for(int i=fi[k];i;i=ne[i])
		  if(v[i]>0 && !dis[w[i]])
		  {
		  	dis[w[i]]=dis[k]+1;q.push(w[i]);
		  }
	}
	return dis[n+m+1];
}

int findd(int u,int vv)
{
	if(u==n+m+1 || !vv) return vv;
	int kkz,tot=0;
	for(int &i=kk[u];i;i=ne[i])
	  if(dis[w[i]]==dis[u]+1 && v[i]>0 && (kkz=findd(w[i],min(v[i],vv-tot))))
	  {
	  	v[i]-=kkz;v[i^1]+=kkz;tot+=kkz;
	  	if(tot==vv) return tot;
	  }
	return tot;
}

int main()
{
	freopen("shuttle.in","r",stdin);
	freopen("shuttle.out","w",stdout);
	scanf("%d%d",&m,&n);cnt=1;
	for(int i=1;i<=m;i++)
	{
		char ch;
		scanf("%d",&x);add(0,i,x);ans+=x;
		while(scanf("%d%c",&x,&ch))
		{
			add(i,x+m,inf);
			if(ch=='\r') break;
		}
	}
	for(int i=1;i<=n;i++) scanf("%d",&x),add(i+m,n+m+1,x);
	while(bfs())
	{
		memcpy(kk,fi,sizeof(fi));
		ans-=findd(0,inf);
	}
	for(int i=1;i<=m;i++) if(dis[i]) printf("%d ",i);printf("\n");
	for(int i=1;i<=n;i++) if(dis[i+m]) printf("%d ",i);printf("\n");
	printf("%d\n",ans);
	fclose(stdin);fclose(stdout);
	return 0;
}