最大流EK算法

最新推荐文章于 2022-08-27 15:06:08 发布

x5gx5gx5g

最新推荐文章于 2022-08-27 15:06:08 发布

阅读量212

点赞数

分类专栏：基本算法文章标签：最大流 EK 算法

本文链接：https://blog.csdn.net/x5gx5gx5g/article/details/85290000

版权

基本算法专栏收录该内容

9 篇文章 0 订阅

订阅专栏

本文深入探讨了最大流问题，并重点解析了Edmonds-Karp算法（EK算法）的工作原理和实现步骤。通过实例分析，解释了如何利用EK算法寻找网络中从源点到汇点的最大传输能力，从而解决实际中的流量分配问题。

摘要由CSDN通过智能技术生成

最大流EK算法

http://blog.csdn.net/x_y_q_/article/details/51999466
//最大流EK算法
//5 4    5是边数 4是节点数
//1 2 40
//1 4 20
//2 4 20
//2 3 30
//3 4 10
//最大流50
#include<stdio.h>
#include<queue>
#include<cstring>
#define arraysize 201
using namespace std;
int maxData = 0x7fffffff;
int capacity[arraysize][arraysize]; //记录残留网络的容量
int flow[arraysize];               //标记从源点到当前节点实际还剩多少流量可用
int pre[arraysize];//标记在这条路径上当前节点的前驱,同时标记该节点是否在队列中
int n,m;                            //n表示边数，m表示结点数
queue<int> myqueue;
int BFS(int src,int des){
    int i;
    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;     //残留图中不再存在增广路径
    else  return 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;
}
int main(){
    int i;int start,end,ci;
    while(~scanf("%d%d",&n,&m)){
        memset(capacity,0,sizeof(capacity));memset(flow,0,sizeof(flow));
        for(i=0;i<n;++i){
            scanf("%d%d%d",&start,&end,&ci);
            if(start == end) continue;            //考虑起点终点相同的情况
            capacity[start][end] +=ci; //此处注意可能出现多条同一起点终点的情况
        }
        printf("%d\n",maxFlow(1,m));
    }
    return 0;
}