最大流EK算法
最大流模板:
#include <iostream> #include <queue> //#include <conio.h> using namespace std; #define arraysize 201 int maxData = 0x7fffffff; int capacity[arraysize][arraysize]; //記錄殘留網絡的容量 int flow[arraysize]; //標記從源點到當前節點實際還剩多少流量可用 int pre[arraysize]; //標記在這條路徑上當前節點的前驅,同時標記該節點是否在隊列中 int n,m; queue<int> myqueue; int BFS(int src,int des) //EK算法使用BFS尋找增廣路徑 {int i,j;while(!myqueue.empty()) //隊列清空 myqueue.pop();for(i=1;i<m+1;++i) //初始化前驅 {pre[i]=-1;}pre[src]=0;flow[src]= maxData; //初始化源點的流量為無窮大 myqueue.push(src);while(!myqueue.empty()){int index = myqueue.front();myqueue.pop();if(index == des) //找到了增廣路徑 break;for(i=1;i<m+1;++i) //遍歷所有的結點 {if(i!=src && capacity[index][i]>0 && pre[i]==-1){pre[i] = index; //記錄前驅 flow[i] = min(capacity[index][i],flow[index]); //關鍵:迭代的找到增量 myqueue.push(i); }}} if(pre[des]==-1) //殘留圖中不再存在增廣路徑 return -1;elsereturn flow[des]; } int maxFlow(int src,int des) {int increasement= 0;int sumflow = 0;while((increasement=BFS(src,des))!=-1){int k = des; //利用前驅尋找路徑 while(k!=src){int last = pre[k];capacity[last][k] -= increasement; //改變正向邊的容量 capacity[k][last] += increasement; //改變反向邊的容量 k = last;}sumflow += increasement;} return sumflow; }總結
- 上一篇: NYOJ 113
- 下一篇: Edmonds_Karp 算法 (转)