[zoj2314]Reactor Cooling——无源汇有上下界可行流

题目大意：

给定一个网络图，每一条有向边有流量的上限和下限，求满足流量守恒下的每一条边的流浪。

思路：

这是一道经典的无源汇有上下界可行流的板子题，也就是所谓的循环流。一个比较主要的思想就是先满足下界之后再进行调整，调整成满足流量守恒和上界的方案。
1. 首先我们先假设每一条边都达到了流量的下界，这样假设后残余网络的值就变成了（上界-下界），流量守恒的条件是每一个点的入流量等于出流量，显然这样来假设流量了之后是不满足流量守恒的，所以我们要在另一个相同的网络里添加一些附加流来满足流量守恒以及不超过残余网络的值。
2. 让我们来假设$a_i$为目前$i$点的（入流量-出流量）的值，相反的在附加网络中$a_i$就变成了（出流量-入流量）的值，我们的目的就是要在附加网络中找到一个合法的流量满足流量差等于$a_i$。
3. 想一想怎么样才可以找到这样一个合法的方案。既然它不满足流量守恒，那我们就在不满足流量守恒的地方添加多余的边使得满足流量守恒，多余的流量可以在图外新建虚拟源点和汇点来满足，如果$a_i>0$，则需要从源点连一条容量为$a_i$的边，如果$a_i<0$，则需要从这个点连一条$-a_i$的边到汇点。
4. 若这个图有合法的方案，那么在这个新图中的表述一定是从原点到汇点的最大流满足把新添加的边给流满，同时还可以合法地在图中添加一些环流，相反的我们就只需要求出从源点到汇点的最大流是否满足流满就好了，在这个基础上添加环流或者不添加环流都是满足上下界的可行的方案。
同时如果有一个办法计算环流之中的流量的话，这个流量一定是最小的。
还是很佩服想出这个做法的人。

/*========================
 * Author : ylsoi
 * Problem : zoj2314
 * Algodithm : MAX_FLOW
 * Time : 2018.7.28
 * =====================*/
#include<bits/stdc++.h>

#define REP(i,a,b) for(int i=a;i<=b;++i)
typedef long long ll;

using namespace std;

void File(){
    freopen("zoj2314.in","r",stdin);
    freopen("zoj2314.out","w",stdout);
}

const int maxn=200+10;
const int maxm=maxn*maxn+10;
const int inf=0x3f3f3f3f;
int T,n,m,ss,tt,ans,a[maxn],sum,low[maxm];
int las[maxm<<1],beg[maxn],to[maxm<<1],flow[maxm<<1],cnte=1;

void add(int u,int v,int f){
    las[++cnte]=beg[u]; beg[u]=cnte; to[cnte]=v; flow[cnte]=f;
    las[++cnte]=beg[v]; beg[v]=cnte; to[cnte]=u; flow[cnte]=0;
}

struct dinic{
    int num[maxn],cur[maxn];
    queue<int>qu;
    bool bfs(){
        memset(num,0,sizeof(num));
        num[ss]=1; qu.push(ss);
        while(qu.size()){
            int u=qu.front(); qu.pop();
            for(int i=beg[u];i;i=las[i]){
                if(!flow[i] || num[to[i]])continue;
                num[to[i]]=num[u]+1;
                qu.push(to[i]);
            }
        }
        return num[tt]!=0;
    }
    int dfs(int u,int res){
        if(u==tt || !res)return res;
        int ret=0,f;
        for(int &i=cur[u];i;i=las[i]){
            if(num[to[i]]!=num[u]+1)continue;
            if((f=dfs(to[i],min(res,flow[i])))){
                res-=f;
                ret+=f;
                flow[i]-=f;
                flow[i^1]+=f;
            }
            if(!res)break;
        }
        return ret;
    }
    void work(){
        while(bfs()){
            REP(i,ss,tt)cur[i]=beg[i];
            ans+=dfs(ss,inf);
        }
    }
}Q;

int main(){
    File();
    scanf("%d",&T);
    while(T--){
        sum=0; ans=0; cnte=1;
        memset(beg,0,sizeof(beg));
        memset(a,0,sizeof(a));
        scanf("%d%d",&n,&m);
        ss=0; tt=n+1;
        int u,v,l,f;
        REP(i,1,m){
            scanf("%d%d%d%d",&u,&v,&l,&f);
            low[i]=l;
            add(u,v,f-l);
            a[v]+=l; a[u]-=l;
        }
        REP(i,1,n){
            if(a[i]>0)sum+=a[i];
            if(a[i]>0)add(ss,i,a[i]);
            else add(i,tt,-a[i]);
        }
        Q.work();
        if(sum!=ans)puts("NO");
        else{
            puts("YES");
            for(int i=2;i<=(m<<1);i+=2)
                printf("%d\n",low[i/2]+flow[i^1]);
        }
    }
    return 0;
}