洛谷P2762_最大权闭合子图

题目：
题目描述
W 教授正在为国家航天中心计划一系列的太空飞行。每次太空飞行可进行一系列商业性实验而获取利润。现已确定了一个可供选择的实验集合E={E1，E2，…，Em}，和进行这些实验需要使用的全部仪器的集合I={I1，I2，…In}。实验Ej需要用到的仪器是I的子集RjÍI。配置仪器Ik的费用为ck美元。实验Ej的赞助商已同意为该实验结果支付pj美元。W教授的任务是找出一个有效算法，确定在一次太空飞行中要进行哪些实验并因此而配置哪些仪器才能使太空飞行的净收益最大。这里净收益是指进行实验所获得的全部收入与配置仪器的全部费用的差额。

对于给定的实验和仪器配置情况，编程找出净收益最大的试验计划。

输入格式
第1行有2 个正整数m和n。m是实验数，n是仪器数。接下来的m 行，每行是一个实验的有关数据。第一个数赞助商同意支付该实验的费用；接着是该实验需要用到的若干仪器的编号。最后一行的n个数是配置每个仪器的费用。

输出格式
第1 行是实验编号；第2行是仪器编号；最后一行是净收益。

输入输出样例
输入 #1 复制
2 3
10 1 2
25 2 3
5 6 7
输出 #1 复制
1 2
1 2 3
17
说明/提示
感谢@FlierKing 提供spj

n,m<=50

这道题数据是在windows生成的，输入数据中所有的换行都是’\r\n’而不是’\n’
读入某实验需要用到的仪器编号的时候，可以这么读入。(感谢@zhouyonglong的提供)

解题思路：
构图：让源点指向实验项目，边权为每个实验项目的收益；然后让没个实验项目指向各自的器材，边权为inf， 最后让每个器材指向汇点，边权为每个器材的消耗。最后跑一边最大流。然后让正的总收益 - 最小割。

注意点：
1.输入可以自己写一个read
2.最后一次bfs中d数组有值的节点就是我们最后选到的实验项目和器材

AC代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>

using namespace std;

const int N = 200, M = 1e4;
const int inf = 0x3f3f3f3f;

int n, m, s, t, ans, maxflow;
int h[N], e[M], w[M], ne[M], now[M], idx;
int d[N];
bool flag;
 
inline void add(int a, int b, int c) {
	e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx ++;
	e[idx] = a, w[idx] = 0, ne[idx] = h[b], h[b] = idx ++;
}

inline bool bfs(void) {
	queue<int> q;
	while(q.size()) q.pop();
	memset(d, 0, sizeof d);
	
	d[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(w[i] && !d[v]) {
				d[v] = d[u] + 1;
				q.push(v);
			}
		}
	}
	
	if(d[t] == 0) return false;
	return true;
}

inline int dinic(int u, int flow) {
	int i;
	
	if(u == t) return flow;
	
	for(i = now[u]; i != -1; i = ne[i]) {
		int v = e[i];
		if(d[v] == d[u] + 1 && w[i]) {
			int k = dinic(v, min(flow, w[i]));
			if(k) {
				w[i] -= k;
				w[i ^ 1] += k;
				return k;
			} else d[v] = 0;
		}
		
		now[u] = i;
	}
	
	return 0;
}

inline void read(int &x) {
	char w = getchar(); x = 0;
	while(w < '0' && w > '9') {
		w = getchar();
	}
	while(w >= '0' && w <= '9') {
		x = x * 10 + w - '0';
		w = getchar();
	} 
	
	if(w == ' ') flag = true;
	else flag = false;
}

int main(void) {
//	freopen("in.txt", "r", stdin);
	
	scanf("%d%d", &m, &n);
	memset(h, -1, sizeof h);
	s = n + m + 1, t = s + 1;
	
	for(int i = 1; i <= m; i ++) {
		int cost, b;
		scanf("%d", &cost);
		getchar();
		add(s, i, cost);
		ans += cost;
		
		flag = true;
		while(flag) {
			read(b);
			add(i, b + m, inf);
		} 
		
	}
	
	for(int i = 1; i <= n; i ++) {
		int cost; scanf("%d", &cost);
		add(m + i, t, cost);
	}
	
	int flow;
	while(bfs()) {
		memcpy(now, h, sizeof h);
		while(flow = dinic(s, inf)) maxflow += flow;
	}

	int cnt = 0, tmp = 0;
	for(int i = 1; i <= m; i ++)
		if(d[i]) cnt++;
	for(int i = 1; i <= m; i ++)
		if(d[i] && tmp < cnt - 1) printf("%d ", i), tmp ++;
		else if(d[i]) printf("%d\n", i);
		
	cnt = tmp = 0;
	for(int i = 1; i <= n; i ++)
		if(d[i + m]) cnt++;
	for(int i = 1; i <= n; i ++)
		if(d[i + m] && tmp <  cnt - 1) printf("%d ", i), tmp++;
		else if(d[i + m]) printf("%d\n", i);

	printf("%d\n",  ans - maxflow);
	
//	fclose(stdin);
	return 0;	
}