pku 2987 Firing 最大权闭合图 （dinic做）
http://poj.org/problem?id=2987
题意：
公司打算裁员，裁掉某些员工可以获得正收益，而裁掉某些员工会遭受损失。并且员工之间往往存在一定的关系，当某个员工被裁掉之后，在他的关系之下的所有员工都必须被裁掉。现在要求如何裁员才能获得最大收益。
思路：
胡波涛的论文里http://wenku.baidu.com/view/986baf00b52acfc789ebc9a9.html
各处了解释，证明的最大权闭合图 = 正权点的和 - 最大割（最大流）；
关键是如何建图，如果该点的权值为正则由源点向该点建立权值为该点权值的边，如果为负则由该点向汇点建立权值为该点权值的绝对值的边，然后原图中的边权值为正无穷。然后就是求最小割了，这里还要输出裁掉了几个员工。我们只要搜索最大流之后还能和源点连接点的个数即可，因为如果最大流之后还能建立连接说明，裁掉他之后的利益大于损失。。
ps:用了自己的最大流dinic模板死活tle，换了一下别人的就A了，好吧，我的模板搓了。
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>
#include <set>
#include <map>
#include <string>
#define CL(a,num) memset((a),(num),sizeof(a))
#define iabs(x) ((x) > 0 ? (x) : -(x))
#define Min(a,b) (a) > (b)? (b):(a)
#define Max(a,b) (a) > (b)? (a):(b)
#define ll long long
#define MOD 1073741824
#define lc l,m,rt<<1
#define rc m + 1,r,rt<<1|1
#define pi acos(-1.0)
#define test puts("<------------------->")
#define maxn 100007
#define M 155000
#define N 5007
using namespace std;
//freopen("din.txt","r",stdin);
const __int64 inf=0x3f3f3f3f3f3f3f3fll;
struct node
{
int u,v;
ll w;
int next;
}g[M];
int head[N],ct;
int q[N*100],l,r;
int vt[N];
int out[N];
ll level[N];
ll zw,val[N];
int num;
int n,m;
void add(int u,int v,ll w)
{
g[ct].u = u;
g[ct].v = v;
g[ct].w = w;
g[ct].next = head[u];
head[u] = ct++;
g[ct].u = v;
g[ct].v = u;
g[ct].w = 0;
g[ct].next = head[v];
head[v] = ct++;
}
bool layer(int s,int e)
{
int i;
CL(level,-1);
level[s] = 0;
int l = 0, r= 0;
q[r] = s;
while (l <= r){
int u = q[l++];
for (i = head[u]; i != -1; i = g[i].next)
{
int v = g[i].v;
if (level[v] == -1 && g[i].w > 0)
{
level[v] = level[u] + 1;
q[++r] = v;
if (v == e) return true;
}
}
}
return false;
}
ll dinic(int s,int e){
ll ans = 0;
while (layer(s,e))
{
int top = 0,u = s,i;
for (i = s; i <= e; ++i) out[i] = head[i];
while (out[s] != -1)
{
if (u == e)
{
ll MIN = inf;
for (i = top - 1; i >= 0; --i)
{
MIN = min(MIN,g[q[i]].w);
}
for (i = top - 1; i >= 0; --i)
{
g[q[i]].w -= MIN;
g[q[i]^1].w += MIN;
if (g[q[i]].w == 0) top = i;
}
ans += MIN;
u = g[q[top]].u;
}
else if (out[u] != -1 && g[out[u]].w > 0 && level[u] + 1 == level[g[out[u]].v])
{
q[top++] = out[u];
u = g[out[u]].v;
}
else
{
while (top > 0 && u != s && out[u] == -1)
u = g[q[--top]].u;
out[u] = g[out[u]].next;
}
}
}
return ans;
}
void dfs(int s)
{
vt[s] = 1; num++;
for (int i = head[s]; i != -1; i = g[i].next)
{
int v = g[i].v;
int w = g[i].w;
if (!vt[v] && w > 0) dfs(v);
}
}
void solve(int s,int e)
{
ll tmp = dinic(s,e);//求最小割
CL(vt,0); num = 0;
dfs(s);//求裁员个数
printf("%d %I64d\n",num - 1,zw - tmp);//这里num多加了0
}
int main()
{
//freopen("din.txt","r",stdin);
int i;
int x,y;
while (~scanf("%d%d",&n,&m))
{
CL(head,-1); ct = 0;
//建图
int s = 0, e = n + 1;
for (i = 1; i <= n; ++i)
{
scanf("%I64d",&val[i]);
if (val[i] >= 0)
{
zw += val[i];
add(s,i,val[i]);//点的权值大于0
}
else//小于0
{
add(i,e,-val[i]);
}
}
for (i = 1; i <= m; ++i)
{
scanf("%d%d",&x,&y);
add(x,y,inf);
}
solve(s,e);
}
return 0;
}