poj 2987 Firing 最大权闭合图 网络流

最大权闭合的求法参看黑书

#include <iostream>#include <cstdio>#include <cstring>using namespace std;const int maxn=5e3+9,maxm=6e4+9;const long long inf=(long long)1<<50;int level[maxn],head[maxn],que[maxn];int lon;struct{    int next,to;    long long w;}e[maxm*10];void edgeini(){    memset(head,-1,sizeof(head));    lon=-1;}void edgemake(int from,int to,long long w){    e[++lon].to=to;    e[lon].w=w;    e[lon].next=head[from];    head[from]=lon;}int makelevel(int s,int t){    memset(level,0,sizeof(level));    level[s]=1;    int top=0;    que[++top]=s;    for(int i=1;i<=top;i++)    {        int u=que[i];        if(u==t) return(1);        for(int k=head[u];k!=-1;k=e[k].next)        if(!level[e[k].to]&&e[k].w)        {            que[++top]=e[k].to;            level[e[k].to]=level[u]+1;        }    }    return(0);}long long dfs(int now,long long maxf,int t){    if(now==t) return(maxf);    long long ret=0;    for(int k=head[now];k!=-1;k=e[k].next)    {        if(e[k].w&&level[e[k].to]==level[now]+1)        {            long long f=dfs(e[k].to,min(maxf-ret,e[k].w),t);            e[k].w-=f;            e[k^1].w+=f;            ret+=f;            if(ret==maxf) return(ret);        }    }    return(ret);}long long dinic(int s,int t){    long long ans=0;    while(makelevel(s,t)) ans+=dfs(s,inf,t);    return(ans);}int flag[maxn];int dfs(int t){    int ret=1;    flag[t]=1;    for(int k=head[t];k!=-1;k=e[k].next)    if(e[k].w>0)    {        int u=e[k].to;        if(!flag[u])        ret+=dfs(u);    }    return(ret);}int main(){    int n,m;    while(scanf("%d %d",&n,&m)!=EOF)    {        edgeini();        long long sum=0;        for(int i=1,tmp;i<=n;i++)        {            scanf("%d",&tmp);            if(tmp>0)            {                sum+=tmp;                edgemake(0,i,tmp);                edgemake(i,0,0);            }            else            {                edgemake(i,n+1,-tmp);                edgemake(n+1,i,0);            }        }        for(int i=1,from,to;i<=m;i++)        {            scanf("%d %d",&from,&to);            edgemake(from,to,inf);            edgemake(to,from,0);        }        memset(flag,0,sizeof(flag));        long long ans=dinic(0,n+1);        printf("%d ",dfs(0)-1);        printf("%lld\n",sum-ans);    }    return 0;}