BZOJ2055 80人环游世界 网络流 费用流 有源汇有上下界的费用流

https://darkbzoj.cf/problem/2055 

https://blog.csdn.net/Clove_unique/article/details/54864211 ←对有上下界费用流的处理方法 

首先建立附加源汇ss,tt 
对于原图里有的一条边x->y,[l,r],cost，变成x->y,r-l,cost 
每一个点的权di定义为所有流入这个点的边的下界和-所有流出这个点的边的下界和 
对于一个点i，若di>0，ss->i,di,0；若di<0，i->tt,-di,0 
连边t->s,inf,0 
然后对ss,tt做最小费用最大流 
最终的费用为（网络流中计算的费用+原图中有费用的边的下界*这条边的费用）

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
const int maxn=250;
const int minf=10000;
int n,m,S,T,SS,mx,ans;
struct nod{
    int x,y,v,co,rev,next;
}e[maxn*maxn];
int head[maxn]={},tot=0;
int dis[maxn]={},fa[maxn]={},vis[maxn]={};
inline void init(int x,int y,int v,int co){
    e[++tot].x=x;e[tot].y=y;e[tot].v=v;e[tot].co=co;e[tot].rev=tot+1;e[tot].next=head[x];head[x]=tot;
    e[++tot].x=y;e[tot].y=x;e[tot].v=0;e[tot].co=-co;e[tot].rev=tot-1;e[tot].next=head[y];head[y]=tot;
}
queue<int>q;
inline bool SPFA(){
    memset(dis,31,sizeof(dis));
    memset(fa,0,sizeof(fa));
    mx=dis[0];
    q.push(S);dis[S]=0;
    while(!q.empty()){
        int x=q.front();q.pop();vis[x]=0;
        for(int i=head[x];i;i=e[i].next){
            int y=e[i].y;
            if((!e[i].v)||dis[y]<=e[i].co+dis[x])continue;
            dis[y]=e[i].co+dis[x];fa[y]=i;
            if(!vis[y]){
                q.push(y);vis[y]=1;
            }
        }
    }
    return dis[T]!=mx;
}
inline void doit(){
    int val=mx;
    for(int i=fa[T];i;i=fa[e[i].x])val=min(val,e[i].v);
    for(int i=fa[T];i;i=fa[e[i].x]){e[i].v-=val; e[e[i].rev].v+=val; ans+=e[i].co*val;}
}
int main(){
    int x;
    scanf("%d%d",&n,&m);S=n*2+1;T=S+1;SS=T+1;
    for(int i=1;i<=n;i++){
        scanf("%d",&x);
        if(x<=0)continue;
        init(i,i+n,0,0);
        init(S,i+n,x,0); init(i,T,x,0);
    }
    for(int i=1;i<=n;i++){
        for(int j=i+1;j<=n;j++){
            scanf("%d",&x);if(x!=-1)init(i+n,j,minf,x);
        }
    }
    init(S,SS,m,0);
    for(int i=1;i<=n;i++) init(SS,i,minf,0);
    ans=0;
    while(SPFA())doit();
    printf("%d
",ans);
    return 0;
}