【POJ 2987】Firing （最小割-最大权闭合子图）

裁员
【问题描述】
在一个公司里，老板发现，手下的员工很多都不务正业，真正干事员工的没几个，于是老板决定大裁员，每开除一个人，同时要将其下属一并开除，如果该下属还有下属，照斩不误。给出每个人的贡献值和从属关系，求在最大贡献值的前提下最小剩下多少人及最大贡献值。留下多少人无所谓，现在老板想知道留下的人最大的贡献值是多少。
【输入描述】
第一行两个整数n,m，表示有多少个员工与多少个从属关系。
第二行n个整数，表示每个员工的贡献值。
接着m行，每行两个数x,y，表示x是y的下属,一个员工可能有多个下属但不会有多个上司
【输出描述】
包括两个数，表示最大贡献值前提下最小剩下多少人及最大贡献值和。
【其他说明】：
0 < n ≤ 5000
0 ≤ m ≤ 60000
员工价值≤107
样例中留下4,5号员工是最好情况

5 5
8
-9
-20
12
-10
1 2
2 5
1 4
3 4
4 5

2 2

1 #include<cstdio>
2 #include<cstdlib>
3 #include<cstring>
4 #include<iostream>
5 #include<algorithm>
6 #include<queue>
7 using namespace std;
8 #define Maxn 5010
9 #define Maxm 60010
10 #define INF 0xfffffff
11 #define LL long long
12
13 struct node
14 {
15     LL x,y,f,o,next;
16 }t[Maxm*2];LL len;
17 LL first[Maxn];
18
19 void ins(LL x,LL y,LL f)
20 {
21     t[++len].x=x;t[len].y=y;t[len].f=f;
22     t[len].next=first[x];first[x]=len;t[len].o=len+1;
23     t[++len].x=y;t[len].y=x;t[len].f=0;
24     t[len].next=first[y];first[y]=len;t[len].o=len-1;
25 }
26
27 LL mymin(LL x,LL y) {return x<y?x:y;}
28
29 LL w[Maxn],st,ed;
30 LL dis[Maxn];
31
32 queue<LL > q;
33 bool bfs()
34 {
35     while(!q.empty()) q.pop();
36     memset(dis,-1,sizeof(dis));
37     q.push(st);dis[st]=0;
38     while(!q.empty())
39     {
40         LL x=q.front();
41         for(LL i=first[x];i;i=t[i].next) if(t[i].f>0)
42         {
43             LL y=t[i].y;
44             if(dis[y]==-1)
45             {
46                 dis[y]=dis[x]+1;
47                 q.push(y);
48             }
49         }
50         q.pop();
51     }
52     if(dis[ed]==-1) return 0;
53     return 1;
54 }
55
56 LL ffind(LL x,LL flow)
57 {
58     if(x==ed) return flow;
59     LL now=0;
60     for(LL i=first[x];i;i=t[i].next) if(t[i].f>0)
61     {
62         LL y=t[i].y;
63         if(dis[y]==dis[x]+1)
64         {
65             LL a=ffind(y,mymin(flow-now,t[i].f));
66             t[i].f-=a;
67             t[t[i].o].f+=a;
68             now+=a;
69         }
70         if(now==flow) break;
71     }
72     if(now==0) dis[x]=-1;
73     return now;
74 }
75
76 LL max_flow()
77 {
78     LL ans=0;
79     while(bfs())
80     {
81         ans+=ffind(st,INF);
82     }
83     return ans;
84 }
85
86 bool vis[Maxn];
87 void dfs(LL x)
88 {
89     vis[x]=1;
90     for(LL i=first[x];i;i=t[i].next) if(!vis[t[i].y]&&t[i].f>0)
91         dfs(t[i].y);
92 }
93
94 int main()
95 {
96     LL n,m;
97     LL ans=0,sum=0;
98     scanf("%lld%lld",&n,&m);
99     for(LL i=1;i<=n;i++) {scanf("%lld",&w[i]);if(w[i]>0) ans+=w[i];}
100     st=n+1,ed=st+1;
101     for(LL i=1;i<=m;i++)
102     {
103         LL x,y;
104         scanf("%lld%lld",&x,&y);
105         ins(x,y,INF);
106     }
107     for(LL i=1;i<=n;i++) if(w[i]>0) ins(st,i,w[i]);
108     for(LL i=1;i<=n;i++) if(w[i]<0) ins(i,ed,-w[i]);
109     LL fl=max_flow();
110     ans-=fl;
111     memset(vis,0,sizeof(vis));
112     dfs(st);
113     for(LL i=1;i<=n;i++) if(vis[i]) sum++;
114     printf("%lld %lld\n",sum,ans);
115     return 0;
116 }

太长了不放了，自己看吧。。