网络流/最大流算法与题目总结

最大流模型为有N个湖（点），他们有小溪相连（边），每个小溪都有它们单位时间的最大流量，最大流求的是一个湖（源点）到另一个湖（汇点）单位时间内能流入的最大流量。

求最大流的一个算法为 ford_fulkerson算法，算法思想：

1.从源点开始寻找一条路径到汇点，记录这条路径流量最小的边的流量为m,然后最大流量+=m，然后每条边的流量都减去m，然后反向更新流量（每条边反方向的流量加上m）；

2.重复步骤1直到找不到路径到汇点为止，最终的最大流量即为所求的最大流。

ford_fulkerson C++模板：

void ford_fulkerson(int s,int e) { queue<int>Q; int u,v; while(1) { Maxflow=0;//最大流初始化 memset(maxflow,0,sizeof(maxflow));//每次寻找增广路径都将每个点的流入容量置为0 memset(visit,0,sizeof(visit));//标记一个点是否已经压入队列 maxflow[s]=INT_MAX;//源点的容量置为正无穷 Q.push(s); // 将源点压入队列 while(!Q.empty()) //当队列不为空 { u=Q.front(); Q.pop(); for(v=1;v<=N;v++) { if(!visit[v]&&flow[u][v]>0) { visit[v]=1; father[v]=u;//记录下他的父亲方便往后的正反向更新 Q.push(v); maxflow[v]=(maxflow[u]<flow[u][v]?maxflow[u]:flow[u][v]);//当前点的容量为父亲点容量与边流量的较小者 } } if(maxflow[e]>0) //如果找到了汇点并且汇点容量不为0则清空队列。 { while(!Q.empty()) Q.pop(); break; } } if(maxflow[e]==0)//已经找不到到汇点的增光路经了，就退出整个循环 break; for(i=e;i!=s;i=father[i]) { flow[father[i]][i]-=maxflow[e];//正向更新 flow[i][father[i]]+=maxflow[e];//反向更新 } Maxflow+=maxflow[e];//更新最大流 } }

一个证明需要反向更新的例子：

 

这个图从1到6明显最的流量为3+3=6，但如果不用反向更新，由于用BFS首先找到一条路径1-2-3-6，然后1-2,2-3,3-6这几条边被更新成了0，就再也找不到增广路径到6了，而1-5-3的剩余的5个流量就无路可走了，用反向更新会再在3-2这里加一条流量为3的边，这时候1-5-3的5流量就可以通过这条边到2再到4再到6了，这就是反向更新的需要。

最大流题目：

POJ1273

最大流入门题，代码：

#include<iostream> #include<queue> using namespace std; const int MAXN = 205; const int maxint = 0x3f3f3f3f; int flow[MAXN][MAXN],maxflow[MAXN],father[MAXN]; int i,Maxflow,M,visit[MAXN]; void ford_fulkerson() { queue<int>Q; int u,v; while(1) { memset(maxflow,0,sizeof(maxflow)); memset(visit,0,sizeof(visit)); maxflow[1]=INT_MAX; Q.push(1); while(!Q.empty()) { u=Q.front(); Q.pop(); for(v=1;v<=M;v++) { if(!visit[v]&&flow[u][v]>0) { visit[v]=1; father[v]=u; Q.push(v); maxflow[v]=(maxflow[u]<flow[u][v]?maxflow[u]:flow[u][v]); } } if(maxflow[M]>0) { while(!Q.empty()) Q.pop(); break; } } if(maxflow[M]==0) break; for(i=M;i!=1;i=father[i]) { flow[father[i]][i]-=maxflow[M]; flow[i][father[i]]+=maxflow[M]; } Maxflow+=maxflow[M]; } } int main() { int N,si,ei,ci; while(cin>>N>>M&&N!=EOF) { Maxflow=0; memset(flow,0,sizeof(flow)); for(i=0;i<N;i++) { cin>>si>>ei>>ci; flow[si][ei]+=ci; } ford_fulkerson(); cout<<Maxflow<<endl; } return 0; }

POJ 1459

这题关键是建图，为该图自建一个源点和汇点就行了。。

POJ 1149

 灰常经典，无敌建图题，最大流必做题。

1.增加源点0与汇点t;

2.如果某顾客有某个猪圈的钥匙，如果之前有顾客有这个猪圈的钥匙，那么加一条之前那顾客到当前这个顾客的边，权值为无限大（方便剩下的猪能传递给当前的顾客），否则就在源点0与这个顾客加一条边，权值为这个猪圈的猪数目。还需要加一条该顾客到汇点t的边，权值为该顾客的需要猪的数量。

3.求源点0到汇点t的最大流