P2762 太空飞行计划问题（最大权闭合子图）

对每一个实验从源点连一条权值为收益的边，每一个仪器向汇点连一条权值为花费的边，实验到仪器连一条边权为inf的边（确保这条边不会被割掉），最终答案就是所有实验的收益和减最小割。在最后一次跑dinic时没有找到一条从源点到汇点的增广路就是说整个图已经被割成两部分了，这时与源点相连的这个图就是我们需要的实验和仪器了，此时与源点相连的实验的边的容量和就是最大收益。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
const int inf = 0x3f3f3f3f;
const int N=1e3+7,M=2e5+10;

int t,n,m,d,S,T;

int h[N],ne[M],w[M],e[M],idx;
int dist[N],dep[N],st[N];
void add(int a,int b,int c){
	e[idx]=b;ne[idx]=h[a];w[idx]=c;h[a]=idx++;
	e[idx]=a;ne[idx]=h[b];w[idx]=0;h[b]=idx++;
}
bool bfs(int s,int t){//分层
	queue<int> q;
	memset(dep,0,sizeof dep);
	dep[s]=1;
	q.push(s);
	while(q.size()){
		int u=q.front();
		q.pop();
		for(int i=h[u];i!=-1;i=ne[i]){
			int v=e[i];
			if(dep[v]==0&&w[i]>0){
				dep[v]=dep[u]+1;
				q.push(v);
			}
		}
	}
	return dep[t]!=0;
}
int cur[N];
int dfs(int u,int flow,int t){
	if(u==t) return flow;
	int sum=0;
	for(int i=cur[u];i!=-1;i=ne[i]){
		cur[u]=i;
		int v=e[i];
		if(dep[v]==dep[u]+1&&w[i]>0){
			int tmp=dfs(v,min(flow,w[i]),t);
			if(!tmp){
				dep[v]=-1;
				continue;
			}
			flow-=tmp;
			sum+=tmp;
			w[i]-=tmp;
			w[i^1]+=tmp;
			if(flow==0) break;
		}
	}
	return sum;
}
int dinic(int s,int t){
	int sum=0;
	while(bfs(s,t)){
		for(int i=0;i<=t;i++) cur[i]=h[i];
		sum+=dfs(s,0x3f3f3f3f,t);
	}
	return sum;
}
int a[N][N],ans[M],ans_cnt,res[M],res_cnt;
void solve(){
	memset(h,-1,sizeof h);
	scanf("%d%d",&m,&n);
	int sum=0;
	for(int i=1;i<=m;++i){
		int x;
		scanf("%d",&x);
		add(S,i,x);
		sum+=x;
		char op;
		while(scanf("%d%c",&x,&op)){
			add(i,m+x,inf);
			if(op=='\n') break;
		}
	}
	T=n+m+1;
	for(int i=1;i<=n;++i){
		int x;
		scanf("%d",&x); add(i+m,T,x);
	}
	int tmp=dinic(S,T);
	for(int i=1;i<=m;++i) if(dep[i]) printf("%d ",i);
	puts("");
	for(int i=m+1;i<T;++i) if(dep[i]) printf("%d ",i-m);
	puts("");
	printf("%d",sum-tmp);
}

int main(){
    t=1;
	while(t--){
		solve();
	}
	return 0;
}