zoj 2314Reactor Cooling

秘制神奇上下界网络流%%%

什么什么有（木）源汇可行流什么的，，看不懂（一下纯属个人sb言论）

看了半天知道点，一个点u，从S连到u的流量是全部流入u的下界，u到T是全部流出u的下界和。（进去出来的约一下）

感觉这个的意思就是保持从进入到出来的下界都符合（强行构造相等？？），并且如果能满流，则上界也符合。那么就是可行的。

看了个有上下界最大流什么的，连一个从T-S的边，然后原图就成了无原汇了，那么再加TT，SS，搞上面的，判断可行的同时可以得出来S-T的流量，是一个可行流量。设为sum1.

那么去掉S-T的边和SS，TT点，在跑了一遍的图上跑最大流，就是使原有的图继续增广，可以得出另一个最大流sum2，那么ans=sum1+sum2。（就看了这么一点，sb）

 #include<bits/stdc++.h>
 #define N 100005
 #define LL long long
 #define inf 0x3f3f3f3f
 #define ls tr[x][0]
 #define rs tr[x][1]
 using namespace std;
 inline int ra()
 {
     int x=,f=; char ch=getchar();
     while (ch<'' || ch>'') {if (ch=='-') f=-; ch=getchar();}
     while (ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}
     return x*f;
 }
 const int S=,T=;
 int n,m,cnt;
 int head[],cur[],h[],q[],in[];
 int low[];
 struct data{int to,next,v;}e[];
 void ine(int x, int y, int v)
 {
     e[++cnt].to=y;
     e[cnt].next=head[x];
     e[cnt].v=v;
     head[x]=cnt;
 }
 void insert(int x, int y, int v)
 {
     ine(x,y,v); ine(y,x,);
 }
 bool bfs()
 {
     for (int i=; i<=T; i++) h[i]=-;
     int l=,r=; q[]=S; h[S]=;
     while (l<r)
     {
         int x=q[l++];
         for (int i=head[x];i;i=e[i].next)
             if (e[i].v && h[e[i].to]==-)
             {
                 h[e[i].to]=h[x]+;
                 q[r++]=e[i].to;
             }
     }
     if (h[T]==-) return ;
     return ;
 }
 int dfs(int x, int f)
 {
     if (x==T) return f;
     int w,ww=;
     for (int i=head[x];i;i=e[i].next)
         if (h[e[i].to]==h[x]+)
         {
             w=dfs(e[i].to,min(e[i].v,f-ww));
             ww+=w; e[i].v-=w; e[i^].v+=w;
             if (ww==f) return f;
         }
     if (!ww) h[x]=-;
     return ww;
 }
 void dinic()
 {
     while (bfs()) dfs(S,inf);
 }
 void build()
 {
     for (int i=; i<=n; i++)
         if (in[i]>) insert(S,i,in[i]);
             else insert(i,T,-in[i]);
 }
 bool jud()
 {
     for (int i=head[S];i;i=e[i].next)
         if (e[i].v) return ;
     return ;
 }
 int main()
 {
     int t=ra();
     while (t--)
     {
         cnt=;
         memset(head,,sizeof(head));
         memset(in,,sizeof(in));
         n=ra(); m=ra();
         for (int i=; i<=m; i++)
         {
             int x=ra(),y=ra(); low[i]=ra(); int w=ra();
             in[x]-=low[i]; in[y]+=low[i];
             insert(x,y,w-low[i]);
         }
         build(); dinic();
         if (!jud()) cout<<"NO"<<endl;
         else{
             cout<<"YES"<<endl;
             for (int i=; i<=m; i++)
                 printf("%d\n",e[(i<<)^].v+low[i]);
         }
         cout<<endl;
     }
     return ;
 }