「网络流 24 题」[2] 太空飞行计划

题面

有m个实验和n种器材。

对于每个器材，都有一个购买价格cost，如果购买了一个器材，可以用于多个实验。

对于每个实验，都有一个费用value和需要的器材的集合s，value表示如果完成这个实验可以获得value的报酬，s表示需要完成这个实验的话就需要s中都所有器材都被购买过。

求最多能获得多少钱(报酬-费用)

分析

首先是可以考虑把所有完成实验的钱都拿过来，然后再考虑哪些东西不买会比较赚。

(懂了懂了懂了)

对于一个实验，要么把所有的器材的费用都流满然后取走value的报酬，要么放弃value的报酬。

如果不同的实验用的器材不重合，那么贪心显然是可行的。

那么如果出现可以重合的器材，那么其实也没有关系。

首先考虑器材不重合的时候的网络流，最后的流量=min(value,sum[S]),用完成所有实验的报酬减去之后就是最优的答案。

那么如果器材可能会有重合的时候要怎么就解决呢？

比如有一个器材，当只有一个实验的时候是没有被流满的，但是来了另外一个实验的钱特别多，于是把这个器材的钱也给流满了，那么原先的这个实验的流量回去流它其他没有流满的器材，然后全满了。

这样对应的情况就是，两个实验共用器材，然后对于实验A来说，完成实验变得更好了。

...当然这么说仍然很奇怪，但是只要理解除非流满，否则这个流量可以表示某个实验费用的一部分。

然后最后判断哪些实验被选择的话就是看这个实验后面的器材是不是都被选了。而且还能有流流到这里(有收益)，那么就选择这个器材。

不过可以通过level直接判断s到这个点是否有空余的容量。

然后器材的话也是同理的。

然后问题就解决啦。建图真神奇！

code

#include<bits/stdc++.h>
#define M 205
#define inf 100000000
using namespace std;
struct ed{
	int x,cap,nx;
}e[M*M];
int nx[M],ecnt;
void add(int x,int y,int cap){
	e[ecnt]=(ed){y,cap,nx[x]};
	nx[x]=ecnt++;
	e[ecnt]=(ed){x,0,nx[y]};
	nx[y]=ecnt++;
}
struct Dinic{
	int nnx[M],level[M],Q[M],s,t;
	bool bfs(int x){
		memset(level,0,sizeof(level));
		level[x]=1;
		int l=0,r=0;
		Q[r++]=x;
		for (;l<r;){
			x=Q[l++];
			for (int i=nx[x];~i;i=e[i].nx)if (e[i].cap>0&&!level[e[i].x]){
				level[e[i].x]=level[x]+1;
				if (e[i].x==t)return 1;
				Q[r++]=e[i].x;
			}
		}
		return level[t]>0;
	}
	int dfs(int x,int f){
		if (x==t)return f;
		int sum=0;
		for (int &i=nnx[x];~i;i=e[i].nx)if (e[i].cap>0&&level[x]+1==level[e[i].x]){
			int d=dfs(e[i].x,min(f-sum,e[i].cap));
			e[i].cap-=d;
			e[i^1].cap+=d;
			sum+=d;
			if (f==sum)return f;
		}
		if (sum==0)level[x]=0;
		return sum;
	}
	int solve(int S,int T){
		s=S; t=T;
		int res=0;
		for (;bfs(s);){
			memcpy(nnx,nx,sizeof(nx));
			res+=dfs(s,inf);
		}	
		return res;
	}
}dinic;
int main(){
//	freopen("1.txt","r",stdin);
	int cost,sum=0,i,n,m;
	scanf("%d%d",&m,&n);
	memset(nx,-1,sizeof(nx));
	for (i=1;i<=m;i++){//实验 
		scanf("%d",&cost);
		add(0,i,cost); sum+=cost;
		char tools[10000];
		memset(tools,0,sizeof tools);
		cin.getline(tools,10000);
		int ulen=0,tool;
		while (sscanf(tools+ulen,"%d",&tool)==1){//tool是该实验所需仪器的其中一个      
		    add(i,m+tool,inf);
		    if (tool==0) 
		        ulen++;
		    else {
		        while (tool) {
		            tool/=10;
		            ulen++;
		        }
		    }
		    ulen++;
		}
		//printf("\n");
	}
	for (i=1;i<=n;i++){//器材 
		scanf("%d",&cost);
		add(i+m,n+m+1,cost);
	}
	int tmp=dinic.solve(0,n+m+1);
	for (i=1;i<=m;i++)if (dinic.level[i]>0)printf("%d ",i); printf("\n");
	for (i=1;i<=n;i++)if (dinic.level[i+m]>0)printf("%d ",i); printf("\n");
	printf("%d\n",sum-tmp);
	return 0;
}