有源汇带上下界最大流

有源汇带上下界最大流

　　Loj 模板

　　本题的网络有源有汇，且每条边有上下界，求最大流。

　　其实就算是有源汇，普通的dinic还是处理有上下界的情况。因此，问题仍需要转化。其实，我们可以在t到s连接一条容量为无穷大的边，那么网络就变成了无源汇的?，无源汇的带上下界的可行流我们会求啊，流程仍然一样，先让初流量都等于下界，然后残流网络跑最大流，如果新建的边都跑满，就存在可行流，否则不存在。不过这道题还多了一个要求，就是求最大流。最大流应该是初流量加上我们能再添加的最大量。而此时的最大量来自s到t之间的一些自由流，这些流当中有些并没有跑满，可以尝试寻找增广路。因此我们只需要把t到s的边删去，跑一下s到t的最大流，加上初流就是答案。

#include<bits/stdc++.h>
using namespace std;
const int N=300;
const int M=5e4+10;
const int inf=1<<29;
int ans,n,m,s,t,S,T,tot=1,sout,flow,maxflow,lin[N],in[M],A[N],d[N];
struct node{int y,v,n;}e[M];
void add(int x,int y,int w){
    e[++tot].y=y;e[tot].n=lin[x];lin[x]=tot;e[tot].v=w;
    e[++tot].y=x;e[tot].n=lin[y];lin[y]=tot;e[tot].v=0;
}
bool bfs(int id){
    memset(d,0,sizeof(d));
    queue<int>q;
    if(id==1) q.push(S),d[S]=1;
    else q.push(s),d[s]=1;
    while(q.size()){
        int x=q.front();q.pop();
        for(int i=lin[x];i;i=e[i].n){
            int y=e[i].y;
            if(e[i].v&&!d[y]){
                d[y]=d[x]+1;
                q.push(y);
                if(id==1&&y==T) return 1;
                if(id==2&&y==t) return 1;
            }
        }
    }
    return 0;
}
int dinic(int x,int flow,int id){
    if(id==1&&x==T) return flow;
    if(id==2&&x==t) return flow;
    int rest=flow;
    for(int i=lin[x];rest&&i;i=e[i].n){
        int y=e[i].y;
        if(e[i].v&&d[y]==d[x]+1){
            int k=dinic(y,min(rest,e[i].v),id);
            if(!k) d[y]=0;
            rest-=k;
            e[i].v-=k;
            e[i^1].v+=k;
        }
    }
    return flow-rest;
}
int main(){
    scanf("%d%d%d%d",&n,&m,&s,&t);
    for(int i=1;i<=m;++i){
        int x,y,l,r;scanf("%d%d%d%d",&x,&y,&l,&r);
        add(x,y,r-l);
        A[y]+=l,A[x]-=l;
        //if(x==s) ans+=l;
    }
    S=0,T=n+1;
    for(int i=1;i<=n;++i){
        if(A[i]>0) add(S,i,A[i]),sout+=A[i];
        else if(A[i]<0) add(i,T,-A[i]);
    }
    add(t,s,inf);
    while(bfs(1))
        while(flow=dinic(S,inf,1)) maxflow+=flow;
    if(maxflow!=sout){
        printf("please go home to sleep\n");
    }else{
        ans+=e[tot].v;
        maxflow=0;
        e[tot].v=e[tot^1].v=0;
        while(bfs(2))
            while(flow=dinic(s,inf,2)) maxflow+=flow;
        ans+=maxflow;
        printf("%d\n",ans);
        //printf("%d\n",maxflow);
    }
    return 0;
}

 

 posted on 2019-01-13 21:13 kgxpbqbyt 阅读( ...) 评论( ...) 编辑 收藏