网络流、费用流、最大流最小割定理……

 囧，今天第一天电脑竞赛补课，就把最大流的BFS增广、先流预推法、最大流最小割定理、最小费用流讲完了。。汗。。

而我，就只记住了BFS增广和最大流最小割定理。最小费用流ms差不多明白了。

所以先讲讲BFS增广求最大流的算法吧。

简单的来说，就是从S(源)开始BFS，直到到达T(汇)or不存在增广路。

所谓增广路就是从S开始到T的一条路径，而且这条路径上的所有边都是非饱和边，即f(i,j)<c(i,j)。因为都是非饱和边，所以这条路径上的边还可以增大经过它们的流，所以称为增广路。而且，这条增广路可以增加的流量，决定于，这条路径上的剩余容量最小的那一条边，可增加的流量就是这最小的剩余流量。

所以说，求最大流的办法就是，不断地寻找增广路，先找到一条增广路，然后把这条路径上的，所有边，减掉，这条增广路可以增大的流，这样重复，直到找不到增广路为止。

下面看看我这段有注释的代码。。

#include<iostream> using namespace std; const long N=110; const long maxint=(1<<31)-1; long n,m,s,t,x,y,z; long g[N][N],pre[N],low[N]; //g[][]表示原始图，pre[i]表示路径中第i个节点的前一个节点 //low[i]表示到达第i个节点时，路径中可增广的流量， //即到第i个节点时，路径上的边剩余容量的最小值。 //low[t]就是整条增广路里的边的剩余容量的最小值， //就是这条路可以增大的流的大小 long q[N],head=0,tail=0; // q[]就是队列，我们用BFS来扩展节点，扩展到了T就表示找到了一条增广路 long Maxflow=0; // 记录最大流的大小 inline void push(long v){ //入队 q[++tail]=v; } inline long pop(){ //出队 return q[++head]; } inline void clear(){ //队列初始化 head=tail=0; } int main(){ cin>>n>>m>>s>>t; for(long i=1;i<=m;i++){ cin>>x>>y>>z; g[x][y]=g[y][x]=z; } while(true){ //初始化 memset(low,-1,sizeof(low)); memset(pre,-1,sizeof(pre)); clear(); push(s);low[s]=maxint; //源点入队，low[s]赋为无穷大，不可能赋其他小的值啊，自己想想为什么？ //开始广搜 while(head<tail){ long x=pop(); for(long i=1;i<=n;i++){ if(g[x][i]>0&&pre[i]<0){ //g[x][i]>0表示这条边还没满啊，pre[i]<0就表示还没前驱，就是还没被扩展啦 push(i); //入队 low[i]=min(low[x],g[x][i]); //然后就考虑这条边的剩余容量是不是比前面所有找到的都小， //是的话就更新呗，不是的话就等于之前的最小的咯。 pre[i]=x; //更新前缀 } } if(pre[t]>0)break; //pre[t]>0就表示T已经被扩展到了，就是说找到一条增广路了，就不用再搜了 } if(pre[t]>0){ //pre[t]>0就是说能找到增广路 long x=t; Maxflow+=low[t]; //这条增广路可以增大的流就是low[t]那么多了，Maxflow加上它，就是当前的最大流 //下面就是把这条增广路上的边全部更新 while(x!=s){ g[x][pre[x]]+=low[t]; //反向边加上这次增加的流，流过来了，那么有机会再流回去。。。 //这句不是太好理解。。 //leokan的话：“ //...为什么不加,不加的话你一开始增广路找错了不就堵塞了? //” g[pre[x]][x]-=low[t]; //正向边减去这次增加的流，就等于剩下的流 x=pre[x]; } }else break; //否则就是找不到增广路了，就推出了。 } cout<<Maxflow; }

 

转自：http://blog.csdn.net/White_Star/archive/2009/01/15/3792083.aspx