SGU 194 Reactor Cooling（无源汇有上下界最大流）

Description
给n个点，及m根管子，每根管子用来流某种液体，管子是单向的，每时每刻每根管子流进来的物质要等于流出去的物质，要使得m条管子组成一个循环体，里面流液体。并且满足每根管子一定的流量限制，范围为[Li,Ri].即要满足每时刻流进来的不能超过Ri(最大流问题)，同时最小不能低于Li。保证没有自环，两相邻点间最多有一根管子。问这个循环体是否存在，如果存在则输出每个管子里的流量，不存在则输出YES
Input
第一行为两个整数n和m表示点数和管子数，之后m行每行四个整数a,b,l,r表示点a到点b之间有一条流量限制为[l,r]的管子(1<=n<=200,0<=l<=r<=10^5)
Output
如果存在满足条件的循环体，则输出YES以及m个管子中的流量，如果不存在则输出NO
Sample Input
4 6
1 2 1 3
2 3 1 3
3 4 1 3
4 1 1 3
1 3 1 3
4 2 1 3
Sample Output
YES
1
2
3
2
1
1
Solution
无源汇有上下界最大流裸题，构造源点s和汇点e，令l(u,v)和r(u,v)为从u到v的管子的流量下限和上限，令f[v]=sum(l(u,v))-sum(l(v,u))，如果f[v]< 0则从v点到汇点e建流量为-f[v]的边，如果f[v]> 0则从源点s到v建流量为f[v]的边，对有管子连接的两点u,v之间建一条流量为r(u,v)-l(u,v)的边，对这张图跑一边最大流，如果与s相连的所有边都满流说明可行流存在，每条边的流量即为该边流量上限减去残余网络中这条边正向弧中的流量，否则说明可行流不存在
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
#define maxn 222
#define maxm 44444
#define INF 0x3f3f3f3f
int head[maxn],cur[maxn],d[maxn],st[maxm],f[maxn],l[maxm],r[maxm];
int s,e,no,n,m,sum;//s为源点,e为汇点,no为边数 
struct point
{
    int u,v,flow,next;
    point(){};
    point(int x,int y,int z,int w):u(x),v(y),next(z),flow(w){};
}p[maxm];
void add(int x,int y,int z)//从x到y建容量为z的边 
{
    p[no]=point(x,y,head[x],z);//前向弧,标号为偶 
    head[x]=no++;
    p[no]=point(y,x,head[y],0);//后向弧,标号为奇 
    head[y]=no++;
}
void init()//初始化 
{
    memset(f,0,sizeof(f));
    memset(head,-1,sizeof(head));
    no=sum=0;
}
bool bfs()
{
    int i,x,y;
    queue<int>q;
    memset(d,-1,sizeof(d));
    d[s]=0; 
    q.push(s);
    while(!q.empty())
    {
        x=q.front();    
        q.pop();
        for(i=head[x];i!=-1;i=p[i].next)
        {
            if(p[i].flow&& d[y=p[i].v]<0)
            {
                d[y]=d[x]+1;
                if(y==e)    
                    return true;
                q.push(y);
            }
        }
    }
    return false;
}
int dinic()//最大流 
{
    int i,loc,top,x=s,nowflow,maxflow=0;
    while(bfs())
    {
        for(i=s;i<=e;i++)   
            cur[i]=head[i];
        top=0;
        while(true)
        {
            if(x==e)
            {
                nowflow=INF;
                for(i=0;i<top;i++)
                {
                    if(nowflow>p[st[i]].flow)
                    {
                        nowflow=p[st[i]].flow;
                        loc=i;
                    }
                }
                for(i=0;i<top;i++)
                {
                    p[st[i]].flow-=nowflow;
                    p[st[i]^1].flow+=nowflow;
                }
                maxflow+=nowflow;
                top=loc;    
                x=p[st[top]].u;
            }
            for(i=cur[x];i!=-1;i=p[i].next)
                if(p[i].flow&&d[p[i].v]==d[x]+1) 
                    break;
            cur[x]=i;
            if(i!=-1)
            {
                st[top++]=i;
                x=p[i].v;
            }
            else 
            {
                if(!top)    
                    break;
                d[x]=-1;
                x=p[st[--top]].u;
            }
        }
    }
    return maxflow;
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        init();
        s=0;
        e=n+1;
        for(int i=0;i<m;i++)
        {
            int u,v;
            scanf("%d%d%d%d",&u,&v,&l[i],&r[i]);
            add(u,v,r[i]-l[i]);
            f[u]-=l[i],f[v]+=l[i]; 
        }
        for(int i=1;i<=n;i++)
            if(f[i]>0)sum+=f[i],add(s,i,f[i]);
            else add(i,e,-f[i]);
        int ans=dinic();
        if(ans!=sum)printf("NO\n");
        else
        {
            printf("YES\n");
            for(int i=0;i<m;i++)
                printf("%d\n",r[i]-p[2*i].flow);
        }
    }
    return 0;
}