POJ 2987 Firing 最大流 网络流 dinic 模板

https://www.cnblogs.com/137shoebills/p/9100790.html

http://poj.org/problem?id=2987

之前写过这道题，码一个dinic的最大流板子。

经典问题，选了一个点就有些点必须选，输出使选出的点的权值和最大的最少点数，并输出该权值和。

建图就是s向权值为正的点连流量为val的边，权值为负的点向t连流量为-val的边，所有的权值和-最大流就是答案。

类似于https://www.cnblogs.com/137shoebills/p/7786985.html

细节是注意流量会爆int

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<queue>
using namespace std;
#define LL long long
const int maxn=5010;
const LL minf=(LL)1e11;
int n,m,s,t;
LL val[maxn]={};
struct nod{
    int y,next;LL v; 
}e[140010];
int head[maxn]={},tot=1,dep[maxn]={};
queue<int>q;
bool vis[maxn]={}; int cnt=0;
void init(int x,int y,LL v){
    e[++tot].y=y;e[tot].v=v;e[tot].next=head[x];head[x]=tot;
}
int bfs(){
    q.push(s);
    memset(dep,0,sizeof(dep));
    dep[s]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=head[x];i;i=e[i].next){
            if(e[i].v&&(!dep[e[i].y])){
                dep[e[i].y]=dep[x]+1;
                q.push(e[i].y);
            }
        }
    }
    return dep[t];
}
LL dfs(int x,LL fc){
    if(x==t)return fc;
    LL tsn=0;
    for(int i=head[x];i;i=e[i].next){
        if(e[i].v&&dep[e[i].y]==dep[x]+1){
            LL z=dfs(e[i].y,min(fc-tsn,e[i].v));
            tsn+=z;e[i].v-=z;e[i^1].v+=z;
            if(tsn==fc)break;
        }
    }
    return tsn;
} 
LL dinic(){
    LL tsn=0;
    while(bfs())tsn+=dfs(s,minf);
    return tsn;
}
void getit(int x){
    if(x==t)return;
    if(x!=s)++cnt;
    vis[x]=1;
    for(int i=head[x];i;i=e[i].next){
        if(e[i].v&&!vis[e[i].y]){
            getit(e[i].y);
        }
    }
}
int main(){
    scanf("%d%d",&n,&m);
    s=n+1;t=s+1;
    LL ans=0;
    for(int i=1;i<=n;++i){
        scanf("%lld",&val[i]);
        if(val[i]>=0){init(s,i,val[i]); init(i,s,0); ans+=val[i];}
        else{init(i,t,-val[i]);init(t,i,0);}
    }
    for(int i=1;i<=m;++i){
        int x,y;
        scanf("%d%d",&x,&y);
        init(x,y,minf);init(y,x,0);
    }
    ans-=dinic();
    getit(s);
    printf("%d ",cnt);
    printf("%lld
",ans);
    return 0;
}