LibreOJ 116 有源汇有上下界最大流 【模板】【上下界最大流】

调完样例后居然一次ac，感动啊...

这是一个模板题有两个点要注意

一是建图的时候要用+=而不能直接=，因为两个点之间的管道可能会多次叠加（特别是处理x和y的时候）

二是在最后算s的流出流量之和时要先把【t到s】和【s到t】的边删掉

像素好低啊。。

#include<iostream>
#include<cstring>
#include<deque>
#define INF 1000000000
using namespace std;

int n,G[250][250],G2[250][250],max_flow,T,sum;

int layer[250],vis[250];
bool countlayer(int s){//还有没有增广路 
    memset(layer,0,sizeof(layer));
    deque<int> q;
    q.push_back(s); layer[s]=1;
    while(!q.empty()){
        int nd = q.front(); q.pop_front();
        if(nd==T) return true;
        for(int i=1;i<=n+2;i++){
            if(G[nd][i] && layer[i]==0){
                layer[i]=layer[nd]+1;
                q.push_back(i);
            }
        }
    }
    return false;
}

int dinic(int s){
    //cout<<"?"<<endl;
    while( countlayer(s) ){
    //    cout<<max_flow<<endl;
    //    cout<<"suc"<<endl;
        memset(vis,0,sizeof(vis));
        deque<int> q;
        q.push_back(s); vis[s]=1;
        while(!q.empty()){
            int nd=q.back();
            if(nd==T){//找到一条增广路 
                int minflow=INF,u;
                for(int i=0;i<q.size()-1;i++){//从上到下的最小管道 
                    int from=q[i],to=q[i+1];
                    if(G[from][to]<minflow){
                        minflow=G[from][to];
                        u=from;
                    }
                }
                //做增广 
                for(int i=0;i<q.size()-1;i++){
                    int from=q[i],to=q[i+1];
                    G[from][to]-=minflow;
                    G[to][from]+=minflow;
                }
                max_flow+=minflow;
                //回溯
                while(!q.empty() && q.back()!=u ){
                    vis[q.back()]=0;
                    q.pop_back();
                }    
            }
            else{//继续往下搜
                int i;
                for(i=1;i<=n+2;i++){
                    if( layer[nd]+1==layer[i] && vis[i]==0 && G[nd][i] ){
                        vis[i]=1;
                        q.push_back(i);
                        break;
                    }
                }
                if(i==n+3) q.pop_back();
            }
        }
    }
}

int main(){
    int m,s,t; cin>>n>>m>>s>>t;
    //x在n+1 ; y在n+2 
    for(int i=1;i<=m;i++){
        int u,v,lower,upper; cin>>u>>v>>lower>>upper;
        sum+=lower;
        G[u][v]+=upper-lower;//非必要弧
        G[u][n+1] += lower;
        G[n+2][v] += lower;//必要弧 
    }
    G[n+1][n+2] = G[t][s] = INF;
    
    for(int i=1;i<=n+2;i++){
        for(int j=1;j<=n+2;j++) G2[i][j]=G[i][j];
    }
    T=n+1;
    dinic(n+2);
    //cout<<max_flow<<" "<<sum<<endl;
    if(max_flow!=sum) cout<<"please go home to sleep"; 
    else{
        int sum1=0;
        //cout<<"!!! "<<G2[s][t]<<" "<<G[s][t]<<endl;
        G[s][t] = G[t][s] = 0;
        for(int i=1;i<=n+2;i++){
        //    cout<<G2[s][i]<<" "<<G[s][i]<<endl;
            sum1+=G2[s][i]-G[s][i];
        }
    //    cout<<sum1<<endl;
        max_flow=0;
        
        T=t;
        dinic(s);
        cout<<sum1+max_flow;
    }
    
    return 0;    
}