SGU 194 Reactor Cooling

http://acm.sgu.ru/problem.php?contest=0&problem=194

题意:m条有向边，有上下界，求最大流。

思路:原图中有u-v low[i],high[i] 建图就是u->v,high[i]-low[i],同时du[u]-=low[i],du[v]+=low[i],然后对于每个点i，如果du[i]>0，那就S->i,du[i],若是du[i]<0那就i->T ,-du[i],然后判断S出去的流是不是满流，不是就无解，否则每条边流量就是基础流量再加上附加流。

啊♂，一开始脑抽了，反向边流量忘记变成0了。。

 #include<algorithm>
 #include<cstdio>
 #include<cmath>
 #include<cstring>
 #include<iostream>
 #define inf 0x7fffffff
 int S,T,nodes;
 int id[];
 int tot,go[],next[],dn[],first[],flow[];
 int dis[],cnt[],op[],ans[],du[],n,m;
 int read(){
     char ch=getchar();int t=,f=;
     while (ch<''||ch>'') {if (ch=='-') f=-;ch=getchar();}
     while (''<=ch&&ch<=''){t=t*+ch-'';ch=getchar();}
     return t*f;
 }
 void insert(int x,int y,int z,int Id){
     tot++;
     go[tot]=y;
     next[tot]=first[x];
     first[x]=tot;
     flow[tot]=z;
     id[tot]=Id;
 }
 void add(int x,int y,int z,int Id){
     insert(x,y,z,Id);op[tot]=tot+;
     insert(y,x,,);op[tot]=tot-;
 }
 int dfs(int x,int f){
     if (x==T) return f;
     int mn=nodes,sum=;
     for (int i=first[x];i;i=next[i]){
         int pur=go[i];
         if (flow[i]&&dis[pur]+==dis[x]){
             int save=dfs(pur,std::min(f-sum,flow[i]));
             sum+=save;
             flow[i]-=save;
             flow[op[i]]+=save;
             if (dis[S]>=nodes||f==sum) return sum;
         }
         if (flow[i]) mn=std::min(mn,dis[pur]);
     }
     if (sum==){
         cnt[dis[x]]--;
         if (cnt[dis[x]]==)
          dis[S]=nodes;
         else
          dis[x]=mn+,cnt[dis[x]]++;
     }
     return sum;
 }
 int main(){
     while (~scanf("%d%d",&n,&m)){
         for (int i=;i<=n+;i++) dis[i]=cnt[i]=first[i]=du[i]=;
         tot=;S=;T=n+;nodes=T+;
         for (int i=;i<=m;i++){
             int u=read(),v=read(),lf=read(),hf=read();
             add(u,v,hf-lf,i);
             du[u]-=lf;
             du[v]+=lf;
             dn[i]=lf;
         }
         for (int i=;i<=n;i++){
             if (du[i]>) add(S,i,du[i],);
             if (du[i]<) add(i,T,-du[i],);
         }
         int Ans;
         while (dis[S]<nodes) Ans+=dfs(S,inf);
         bool flag=true;
         for (int i=first[S];i;i=next[i])
          if (flow[i]!=){
                 flag=false;
                 break;
          }
         if (!flag) {
             puts("NO");
             continue;
         }else{
             puts("YES");
             for (int i=;i<=tot;i++) ans[id[i]]=flow[op[i]];
             for (int i=;i<=m;i++) printf("%d\n",dn[i]+ans[i]);
         }
     }
 }