POJ 2987 最大权闭合图

题意：

 一个公司打算解雇一些员工，已知解雇每个人的收益值（可能为负）。解雇某一个人后，他的下属也将被解雇，求最大收益是多少，在取得最大收益时需要解雇多少人。

题解：

建图不多啰嗦，挺裸的~

求最少的割边数目，可以从源点对残量网络进行一次DFS，每个割都会将源点和汇点隔开，所以从源点DFS下去一定会碰到某个割dfs停止，这时遍历过的点数最少的裁员数目。那种直接找满流边的做法是显然不对的奥~

ps：总是wa盛怒之下，把所有int替换了ac了。。。

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>

#define N 100100
#define M 3000000
#define INF 100000000

using namespace std;
//理解最大权闭合图的最小割的含义 
__int64 head[N],to[M],next[M],len[M],sum,n,m,cnt,S,T,layer[N],q[M<<2],num,ans;
bool vis[N];

inline void add(__int64 u,__int64 v,__int64 wp)
{
    to[cnt]=v; len[cnt]=wp; next[cnt]=head[u]; head[u]=cnt++;
    to[cnt]=u; len[cnt]=0; next[cnt]=head[v]; head[v]=cnt++;
}

void read()
{
    memset(head,-1,sizeof head);cnt=0;
    S=0,T=n+1;sum=0;
    for(__int64 i=1,a;i<=n;i++)
    {
        scanf("%I64d",&a);
        if(a>=0)
        {
            sum+=a;
            add(S,i,a);
        }
        else add(i,T,-a);
    }
    for(__int64 i=1,a,b;i<=m;i++)
    {
        scanf("%I64d%I64d",&a,&b);
        add(a,b,INF);
    }
}

bool bfs()
{
    memset(layer,-1,sizeof layer);
    __int64 h=1,t=2,sta;
    q[1]=S; layer[S]=0;
    while(h<t)
    {
        sta=q[h++];
        for(__int64 i=head[sta];~i;i=next[i])
            if(len[i]&&layer[to[i]]<0)
            {
                layer[to[i]]=layer[sta]+1;
                q[t++]=to[i];
            } 
    }
    return layer[T]!=-1;
}

__int64 find(__int64 u,__int64 cur_flow)
{
    if(u==T) return cur_flow;
    __int64 result=0,tmp;
    for(__int64 i=head[u];~i&&result<cur_flow;i=next[i])
        if(len[i]&&layer[to[i]]==layer[u]+1)
        {
            tmp=find(to[i],min(cur_flow-result,len[i]));
            len[i]-=tmp; len[i^1]+=tmp; result+=tmp;
        }
    if(!result) layer[u]=-1;
    return result;
}

void dinic()
{
    ans=0;
    while(bfs()) ans+=find(S,INF);
    
}

void dfs(__int64 u)
{
    if(u==T) return;
    vis[u]=true;
    for(__int64 i=head[u];~i;i=next[i])
        if(!vis[to[i]]&&len[i]>0)
            dfs(to[i]);
}

int main()
{
    while(scanf("%I64d%I64d",&n,&m)!=EOF)
    {
        read();
        dinic();
        num=0;
        memset(vis,0,sizeof vis);
        dfs(S);
        for(__int64 i=1;i<=n;i++)
            if(vis[i]) num++;
        printf("%I64d %I64d\n",num,sum-ans);
    }
    return 0;
}