poj 2987 Firing (最大流、最大權閉合圖)

傳送門：http://poj.org/problem?id=2987

題目大意：

  公司裁員，每裁一個員工有相應的收益(正負都有可能)，當一個員工被裁，他的下屬也全會被裁。求最大的收益以及此時被裁的最少人數。

思路：

  首先增加一個源點S和彙點T，S指向所有收益>0的點，所有收益<0的點指向T，再将每個員工/上司關系以上司指向員工建圖，這樣來建構一個二分圖，求它的最大權閉合圖。

  由于最大權閉合圖=最大權之和(所有正的權相加)-最小割，而根據最大流最小割定理可以知道最大流就等于最小割，是以可以轉變為求最大流。

  統計人數隻需要從S搜尋沒有滿流的那一塊的點數即可。

  另外需要注意開long long。。

AC代碼：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
#include<utility>
#include<algorithm>
#include<utility>
#include<queue>
#include<vector>
#include<set>
#include<stack>
#include<cmath>
#include<map>
#include<ctime>
#define P pair<int,int>
#define ll long long
#define ull unsigned long long
#define lson id*2,l,mid
#define rson id*2+1,mid+1,r
#define ls id*2
#define rs id*2+1
#define Mod(a,b) a<b?a:a%b+b 
using namespace std;

const ll M = 1e9 + 7;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const int N = 5010;
const double e = 10e-6;

struct edge 
{ 
	int to;
	ll cap;
	int rev;
};
vector<edge>G[N];
int level[N], iter[N];

void addEdge(int from, int to, ll cap)
{
	G[from].push_back(edge{ to,cap,(int)G[to].size() });
	G[to].push_back(edge{ from,0,(int)G[from].size() - 1 });
}

void bfs(int s)
{
	memset(level, -1, sizeof(level));
	queue<int> que;
	level[s] = 0;
	que.push(s);
	while (!que.empty()) {
		int v = que.front(); que.pop();
		int len = G[v].size();
		for (int i = 0; i < len; i++) {
			edge &e = G[v][i];
			if (e.cap > 0 && level[e.to] < 0) {
				level[e.to] = level[v] + 1;
				que.push(e.to);
			}
		}
	}
}

ll dfs(int v, int t, ll f)
{
	if (v == t)
		return f;
	int len = G[v].size();
	for (int &i = iter[v]; i < len; i++) {
		edge &e = G[v][i];
		if (e.cap > 0 && level[v] < level[e.to]) {
			ll d = dfs(e.to, t, min(f, e.cap));
			if (d > 0) {
				e.cap -= d;
				G[e.to][e.rev].cap += d;
				return d;
			}
		}
	} 
	return 0;
}

ll maxFlow(int s, int t)
{
	ll flow = 0;
	while (1) {
		bfs(s);
		if (level[t] < 0)
			return flow;
		memset(iter, 0, sizeof(iter));
		ll f;
		while ((f = dfs(s, t, INF)) > 0)
			flow += f;
	}
}

int n, m, num;
bool vis[N];

void cal(int u)
{
	vis[u] = true;
	num++;

	int len = G[u].size();
	for (int i = 0; i < len; i++) {
		edge e = G[u][i];
		if (e.cap > 0 && !vis[e.to])
			cal(e.to);
	}
}

int main() {
	while (~scanf("%d%d", &n, &m)) {
		for (int i = 0; i <= n + 1; i++)
			G[i].clear();
		
		ll xx, ans = 0;
		for (int i = 1; i <= n; i++) {
			scanf("%lld", &xx);
			if (xx > 0) {
				addEdge(0, i, xx);
				ans += xx;
			}
			else if (xx < 0)
				addEdge(i, n + 1, -xx);
		}

		int x, y;
		for (int i = 1; i <= m; i++) {
			scanf("%d%d", &x, &y);
			addEdge(x, y, INF);
		}
		
		ans -= maxFlow(0, n + 1);

		memset(vis, false, sizeof(vis));
		num = 0;
		cal(0);

		printf("%d %lld\n", num - 1, ans);
	}

	return 0;
}