最大权闭合子图_洛谷3410

题目：
题目描述
小B有n个下属，现小B要带着一些下属让别人拍照。

有m个人，每个人都愿意付给小B一定钱让n个人中的一些人进行合影。如果这一些人没带齐那么就不能拍照，小B也不会得到钱。

注意：带下属不是白带的！！！对于每个下属，如果他带了那么小B需要给他一些钱，保证当他拍照时配合。

请问，小B的净收益最多是多少。

输入格式
第1行有2个正整数m和n（0<m,n<=100）。接下来的m行，每行是一个要求拍照的人的有关数据。第一个数是他同意支付该合影的费用；接着是该合影需要的若干下属的编号，以一个0作为行的结束标记。最后一行的n个数是带每个下属的费用。

输出格式
一个数，表示最大收益。小B可以一个人也不带。

输入输出样例
输入 #1 复制
2 3
10 1 2 0
25 2 3 0
5 6 7
输出 #1 复制
17
说明/提示
对于10%的数据每个人都要求让全部n个人合影

对于30%的数据n<=15 m<=15

另有10%的数据答案为0

对于50%的数据n<=40 m<=40

另有10%的数据每个人只愿意拍一个人

对于100%的数据m,n<=100

最大权闭合子图：
有一个有向图，每一个点都有一个权值（可以为正或负或0），选择一个权值和最大的子图，使得每个点的后继都在子图里面，这个子图就叫最大权闭合子图。

解题思路1：
有一些点，有些点的权值是正的有一些点的权值是负数，然后所有的正的点都指向一些权值是负的点，这道题让我们去得到收益的最大值，也就是让我们去找到一个最大权闭合子图，求出这个子图的收益之和。所以我们可以让s 指向m个点，边权分别为对应的收益，然后这m个然分别指向对应的n个人中的一些人，边权为inf， 最后让这个n个人指向t，边权为这n个人的花费，最后 答案 ans = 有效的收益 - 最小割。

解题思路（相当于解题思路一中结论的证明）：
首先构图和上面的想法是一样的，然后我们不妨这样想，我们需要在A地和B地里种菜，这些菜要么只能在A地里获得收益，要么只能在B地里获得收益，然后我们去求这个图的最小割，即用最小的代价让s 和 t 不连通。那么我们最后 割下的边要么是题目中正收益的边，要么是负收益的边；被割下的正收益的边表示我选这条边反而会是我的收益下降所以割下了，割下的负收益的边表示这些边是我选取正收益的边本来就应该支付的代价，那么用我们的全部的正收益之和 - 最小的代价和 就是我们想要得到的答案。

注意点：
1.记得给flow 赋值，即while(flow = dinic(s, inf)) 。
AC代码：

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

using namespace std;

const int N = 300, M = 30000;
const int inf = 0x3f3f3f3f;

int n, m, ans, maxflow, s, t;
int h[N], w[M], e[M], ne[M], now[M], idx;
int d[N];

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;
}

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