poj2987 Firing 最大权闭合图

题意：某个公司需要裁员，如果要裁u,那么u的下属也要被裁。裁掉一个人可对应获取一个值，求尽可能裁少的人获得最大的值。

思路：最大权闭合图。建立源点S和T，若点i的w权值为正,S向点i连一条容量为w的边，否则，i向T连一条容量为-w的边，对于有向

边i->j，i向j连一条容量为无穷的边。那么最大权和=所有权值为正的值的总和 - 最小割。（证明可以参见Amber大牛的《最小割模型

在信息学竞赛中的应用》，非常详细） 详见代码：

// file name: poj2987.cpp //
// author: kereo //
// create time:  2014年09月04日 星期四 07时05分52秒 //
//***********************************//
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<stack>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MAXN=5000+100;
const int M=60000*2+100;
const int inf=0x3fffffff;
#define L(x) (x<<1)
#define R(x) (x<<1|1)
int n,m,edge_cnt,top,cnt;
int que[MAXN],s[MAXN],vis[MAXN];
int head[MAXN],depth[MAXN],cur[MAXN],gap[MAXN];
struct Edge
{
	int cap,flow,v,next;
}edge[M];
void init(){
	edge_cnt=top=cnt=0;
	memset(head,-1,sizeof(head));
}
void addedge(int u,int v,int cap){
	edge[edge_cnt].v=v; edge[edge_cnt].cap=cap;
	edge[edge_cnt].flow=0; edge[edge_cnt].next=head[u]; head[u]=edge_cnt++;
	edge[edge_cnt].v=u; edge[edge_cnt].cap=0;
	edge[edge_cnt].flow=0; edge[edge_cnt].next=head[v]; head[v]=edge_cnt++;
}
void bfs(int st,int ed){
	memset(depth,-1,sizeof(depth));
	memset(gap,0,sizeof(gap));
	int front=0,rear=0; gap[0]=1; //因为起点,gap[0]=1
	depth[ed]=0; que[rear++]=ed;
	while(front!=rear){
		int u=que[front++];
		for(int i=head[u];i!=-1;i=edge[i].next){
			int v=edge[i].v;
			if(depth[v]!=-1) continue;
			que[rear++]=v;
			depth[v]=depth[u]+1; gap[depth[v]]++;
		}
	}
}
ll isap(int st,int ed){
	bfs(st,ed);
	memcpy(cur,head,sizeof(head));
	top=0;
	int u=st;
	ll ans=0;
	while(depth[u]<n+2) {
		if(u == ed){
			int Min=inf,ins; //下一次从ins之前的边插入
			for(int i=0;i<top;i++)
				if(edge[s[i]].cap-edge[s[i]].flow<Min){
					Min=edge[s[i]].cap-edge[s[i]].flow;
					ins=i;
				}
			for(int i=0;i<top;i++){
				edge[s[i]].flow+=Min; edge[s[i]^1].flow-=Min;
			}
			ans+=Min; top=ins; u=edge[s[top]^1].v;
			continue;
		}
		int v,flag=0; //记录是否能往下走
		for(int i=cur[u];i!=-1;i=edge[i].next){
			v=edge[i].v;
			if(edge[i].cap-edge[i].flow && depth[v]+1 == depth[u]){
				flag=1; cur[u]=i;
				break;
			}
		}
		if(flag){
			s[top++]=cur[u]; u=v;
			continue;
		}
		int d=n+2;
		for(int i=head[u];i!=-1;i=edge[i].next) //找最小的depth
			if(edge[i].cap-edge[i].flow && depth[edge[i].v]<d){
				d=depth[edge[i].v]; cur[u]=i;
			}
		gap[depth[u]]--;
		if(!gap[depth[u]]) return ans; //gap优化
		depth[u]=d+1; gap[depth[u]]++;
		if(u!=st) u=edge[s[--top]^1].v;
	}
	return ans;
}
void dfs(int u){ //dfs确定最小割
	vis[u]=1;
	for(int i=head[u];i!=-1;i=edge[i].next){
		int v=edge[i].v;
		if(!vis[v] && edge[i].cap-edge[i].flow){
			dfs(v); cnt++;
		}
	}
}
int main()
{
	while(~scanf("%d%d",&n,&m)){
		init();
		ll sum=0;
		for(int i=1;i<=n;i++){
			int x;
			scanf("%d",&x);
			if(x>0) 
				addedge(0,i,x),sum+=x;
			else if(x<0)  
				addedge(i,n+1,-x);
		}
		for(int i=0;i<m;i++){
			int u,v;
			scanf("%d%d",&u,&v);
			addedge(u,v,inf);
		}
		sum-=isap(0,n+1);
		memset(vis,0,sizeof(vis));
		dfs(0);
		printf("%d %I64d\n",cnt,sum);
	}
	return 0;
}