【转】上下界流

建一张新图添加两个新点源点S和汇点T，若原图上存在边I-J，最大流量为A，最小为B，则新图上建边S-J容量为A，I-J容量为B-A，I-E容量为A。再新加入一条边为原图上的汇点到源点容量为正无穷。在新图上求最大流，若源点流出的边全部满流则有解，把新图上的流量加上原图的最低流量就是原图上的流量，把所有原图汇点流出的流全部加起来就是最小流，若求最大流则在原图上保证流量大于最低流量的情况下增广。

#include<iostream> const int maxint=0x7ffffff; int from[5100],max,n,map[510][510],cop[510][510],s,e,map2[510][510],tp[510][510],ll[510],cl,op,list[510*510]; bool used[510]; int calc_min(int a,int b) { return a<b?a:b; } void zengguang() { int t=ll[e],p=op; max+=t; while (from[p]!=0) { int a=list[p]; if (from[p]>0) { int b=list[from[p]]; tp[b][a]+=t; p=from[p]; } else { int b=list[-from[p]]; tp[a][b]-=t; p=-from[p]; } } } int find() { memset(used,0,sizeof(used)); ll[s]=maxint,cl=1,op=2,used[s]=1,list[1]=s; while (cl<op) { int a=list[cl]; for (int i=1;i<=e;i++) { int b=i; if (used[b]) continue; if (map2[a][b]-tp[a][b]) { ll[b]=calc_min(ll[a],map2[a][b]-tp[a][b]); list[op]=b; from[op]=cl; used[b]=1; if (b==e) return 1; op++; continue; } if (tp[b][a]) { ll[b]=calc_min(ll[a],tp[b][a]); list[op]=b; from[op]=-cl; used[b]=1; if (b==e) return 1; op++; continue; } } cl++; } return 0; } int main(void) { int a,b,c,d; scanf("%d",&n); s=n+1,e=s+1; map2[n][1]=maxint; while (scanf("%d%d%d%d",&a,&b,&c,&d)==4) { map[a][b]+=d; cop[a][b]+=c; map2[s][b]+=c; map2[a][b]+=d-c; map2[a][e]+=c; } while (find()) zengguang(); max=0; for (int i=1;i<=n;i++) max+=tp[i][n]+cop[i][n]; printf("%d",max); return 0; } 

 

参考资料：

1. 容量有上下界的网络流_信息学NOI联赛辅导专题 王建德讲稿

2.http://blog.sina.com.cn/s/blog_4e7bbc500100ct4u.html###