hdu3061 最大权闭合集

闭合图定义:一个有向图G=(V,E)的闭合图是该有向图的一个点集，且该点集的所有出边都还指向该点集

另外 闭合图可以是几个联通图的合体的部分子图 

如果给每一个点一个权值的话 就可以求最大权闭合图

解法：另外建立一个回答你和源点

按照权值分类：把权值为正的点和源点相连接 然后权值就是该点的权值 把权值为负的点和汇点相连接 权值就是该点权值的相反数 

最后是托普关系 反拓扑序列建图 i->j  图上就是 j－>i 跑最大流就可以了

#include<algorithm>
#include<cstdio>
#include<vector>
#include<iostream>
using namespace std;
const int MAXn=500+9;
int n,answer,d[MAXn],gap[MAXn],st=1,ed;
struct edge{int to,cap,rev;};
vector<edge>g[MAXn];
void add(int u,int v,int cap){
    g[u].push_back((edge){v,cap,(int)(g[v].size())});
    g[v].push_back((edge){u,0,(int)(g[u].size()-1)});
}
int sap(int u,int flow)
{
    if(u==ed)
        return flow;
    int res=flow,t;
    for(int i=0;i<g[u].size();i++){
        edge &e=g[u][i];
        if(e.cap && d[u]==d[e.to]+1)
        {
            t=sap(e.to,min(res,e.cap));
            e.cap-=t,g[e.to][e.rev].cap+=t;
            if(!(res-=t))
                return flow;
        }
    }
    if(!--gap[d[u]])
        d[st]=n;
    ++gap[++d[u]];
    return flow-res;
}
int main()
{
    int m,i,t1,t2,t3;
    while(~scanf("%d%d",&n,&m))
    {
        int x;
        n+=2;
        for(int i=0;i<=n;++i)g[i].clear();
        ed=n;
        st=n-1;
        answer=0;
        int sum=0;
        for(int i=1;i<=n-2;++i)
        {
            scanf("%d",&x);
            if(x>0){sum+=x;add(st,i,x);}
            else add(i,ed,-x);
        }
        memset(d,0,sizeof d),memset(gap,0,sizeof gap);
        for(i=1;i<=m;++i)
        {
            scanf("%d%d",&t1,&t2);
            add(t1,t2,~0U>>1);
        }
        for(gap[0]=n;d[st]<n;)
            answer+=sap(st,~0U>>1);
        
        printf("%d\n",sum-answer);
    }
    return 0;
}