有上下界的网络流问题学习笔记
刚刚学了有上下界的网络流问题,总结一下。
大概可以分为4种:
1、无源汇的可行流
2、有源汇的可行流
3、有源汇的最大流
4、有源汇的最小流
无源汇的可行流
想像一条水循环系统,无源汇的可行流就是流量在整张图里循环,每一个节点都满足流量守恒,没有源点和汇点这样的特殊点。
但是,在这个问题中,每条边都有一个容量下限Bi和一个容量上限Ci。上限好理解,就是普通最大流中的东西,但下限怎么处理?
想象,如果在这张流网络满足每条边的流量都在【Bi,Ci】内,那么对于每个节点,流经它的流量至少为它的入边的下限和,流出它的流量至少为它出边的下限和。
怎么限制住这个条件呢?于是我们人为地加入一个源点S和汇点T:
S:向每个节点连一条边,容量为它入边的下限和
T:每个节点向T连一条边,容量为它出边的下限和
对每条其它的边,容量变为Ci-Bi,即为正常的最大流的边。
为什么能这样呢?
每条边上的流量可以分为两种:限制流和自由流
限制流就是为了满足下限而必须存在的流,也就是该边的下限
自由流即超出下限的那一部分
我们加入S和T,砍掉其它边的下限,相当于把其它边的限制流抢了过来,由S和T管理,这样子剩下的就是自由流,可以为任意流量。
这样一来,如果该流网络存在可行流,当且仅当求出S到T的最大流后S的出边都满载,T的入边都满载。
ZOJ2314 Reactor Cooling
#include<iostream> #include<cstdio> #include<queue> #include<algorithm> using namespace std; const int maxn=205,maxm=80005,INF=2000000000; int inf[maxn],outf[maxn]; inline int read() { int out=0,flag=1;char c=getchar(); while(c<48||c>57) {if(c=='-') flag=-1;c=getchar();} while(c>=48&&c<=57) {out=out*10+c-48;c=getchar();} return out*flag; } int head[maxn],nedge=0; class EDGE { public: int to,next,f,b; }edge[maxm]; inline void build(int a,int b,int w,int l) { edge[nedge]=(EDGE){b,head[a],w,l}; head[a]=nedge++; edge[nedge]=(EDGE){a,head[b],0,l}; head[b]=nedge++; } int cur[maxn],S,T,d[maxn]; bool vis[maxn]; bool bfs() { fill(vis,vis+maxn,false); queue<int> q; q.push(S); d[S]=0;vis[S]=true; int u,to; while(!q.empty()) { u=q.front(); q.pop(); for(int k=head[u];k!=-1;k=edge[k].next) if(!vis[to=edge[k].to]&&edge[k].f) { d[to]=d[u]+1; vis[to]=true; q.push(to); } } return vis[T]; } int dfs(int u,int minf) { if(u==T||!minf) return minf; int flow=0,f,to; if(cur[u]==-2) cur[u]=head[u]; for(int& k=cur[u];k!=-1;k=edge[k].next) if(d[to=edge[k].to]==d[u]+1&&(f=dfs(to,min(minf,edge[k].f)))) { edge[k].f-=f; edge[k^1].f+=f; flow+=f; minf-=f; if(!minf) break; } return flow; } int maxflow() { int flow=0; while(bfs()) { fill(cur,cur+maxn,-2); flow+=dfs(S,INF); } return flow; } int main() { int CNT=read(); while(CNT--) { fill(head,head+maxn,-1); fill(inf,inf+maxn,0); fill(outf,outf+maxn,0); nedge=0; int N=read(),M=read(),a,b,l,r; S=0;T=N+1; for(int i=1;i<=M;i++) { a=read(); b=read(); l=read(); r=read(); outf[a]+=l; inf[b]+=l; build(a,b,r-l,l); } int cnt=0,flow; for(int i=1;i<=N;i++) build(S,i,inf[i],0),cnt+=inf[i]; for(int i=1;i<=N;i++) build(i,T,outf[i],0); flow=maxflow(); //cout<<cnt<<' '<<flow<<endl; if(cnt==flow) { printf("YES\n"); for(int i=0;(i>>1)<M;i+=2) printf("%d\n",edge[i^1].f+edge[i^1].b); } else printf("NO\n"); } return 0; }
有源汇的可行流
与无源汇的可行流类似,有源汇的可行流就是多了个源点和汇点。这个时候我们可以从另一个角度来看网络流:虽然S和T不满足流量限制,但如果我们人为连一条由T到S的容量为正无穷的边,那么它就转化成了所有点都遵循容量限制的无源汇的流网络。
那么我们再加一个超级源点S‘和超级汇点T’【就相当于上边的人为加入的源汇点】,就可以转化为无缘汇的可行流啦。
有源汇的最大流
有源汇的最大流就是在求出有源汇的可行流之后去掉加入的T到S的无穷大的边,然后再在求完可行流的残量网络中求一遍由S到T的最大流即为最大流
最后结果为可行流中边T->S的流量+新一次最大流中新增的流量。
有源汇的最小流
求法:
1、同样加一个超级源S‘,超级汇T’,【这个时候不加边T->S】求一遍S‘到T’最大流
2、加入边T->S【正无穷】,再求一遍S‘到T’最大流。
3、若S‘出去的边不满流,T’进来的边不满流,则无解,否则解为T->S边的流量。【这些流量为使可行流成立必须加的流】
我很辣鸡,只知道这些。