网络流Edmonds-Karp算法入门

今天自习课没事干,看书自学了一下网络流中的EK算法。(求最大流)

设s为源点,t为汇点,C为容量矩阵,F为流量矩阵,f为最大流量。

1.初始化F,f

2.用BFS在残量网络中找到一条从s到t的最短增广路T,如果T不存在,算法结束。

  最短增广路就是最短路径(s→t)

  一边(u,v)被增广完后,在残量网络中就不存在(u,v)这条边,只有(v,u)这条边了。(即不存在(u,v)这条路径)

3.m=min(u,v)∈T{C(u,v)-F(u,v)} (最短增广路径中残量网络中最小的边)

4.f+=m

5.沿着T修改矩阵,对于任意(u,v)∈T,将F(u,v)增加m,F(u,v)减少m

  F(u,v)之所以要减少m是因为有可能会出现以下情况:

  假设求出增广路径1→2→3→4

  但发现其实1→3→4+1→2→4更优

  这就需要反向建边

6.repeat 2.

 

code:

#include <cstdio>
#include <cstring>
#include <algorithm>
#define F(i,x,y) for(int i=x;i<=y;i++)
using namespace std;

const int MAXN=201;
int N,M,a[MAXN][MAXN];
int pre[MAXN],vis[MAXN],l[MAXN],h,t,S,T;

int BFS()
{
    memset(pre,-1,sizeof(pre));
    memset(vis,0,sizeof(vis));
    h=t=0;pre[S]=S,vis[S]=1;l[++t]=S;
        while(h<t){
            int P=l[++h];
            for(int i=1;i<=N;i++){
                if(a[P][i]>0&&!vis[i]){
                    vis[i]=1;
                    pre[i]=P;
                    if(i==T)return true;
                    l[++t]=i;
                }
            }
        }
    return false;
}

int EK()
{
    int F=0,w=2e9;
        while(BFS()){
            for(int i=T;i!=S;i=pre[i])w=min(w,a[pre[i]][i]);
            for(int i=T;i!=S;i=pre[i]){
                a[i][pre[i]]+=w;
                a[pre[i]][i]-=w;
            }
            F+=w;
        }
    return F;
}

int main()
{
    while(scanf("%d%d%d%d",&N,&M,&S,&T)!=EOF){
        memset(a,0,sizeof(a));
        int u,v,w;
            F(i,1,M){
                scanf("%d%d%d",&u,&v,&w);
                a[u][v]+=w;
            }
        printf("%d\n",EK());
    }return 0;
}

转载于:https://www.cnblogs.com/Cptraser/p/7921455.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
### 回答1: Edmonds-Karp算法是一种最大算法,可以在网络问题中找到最大量。Matlab是一款数学软件,可以用于求解各种数学问题,包括网络问题。可以使用Matlab编写Edmonds-Karp算法的程序或者将现有的Edmonds-Karp算法移植到Matlab中进行求解。 ### 回答2: Edmonds-Karp算法是解决最大问题的一种经典算法。它是基于Ford-Fulkerson算法的改进版本,通过在寻找增广路径时使用BFS(广度优先搜索)来保证找到的路径是最短路径。以下是一个使用Matlab实现Edmonds-Karp算法的伪代码: 1. 输入: 邻接矩阵adjMatrix,表示图的连接关系;起始节点source和目标节点target。 2. 初始化残余容量矩阵residualCapacity,将其与adjMatrix相同并全部初始化为0。 3. 初始化最大量maxFlow为0。 4. 循环直到无法找到增广路径: 5. 使用BFS从source到target寻找增广路径,得到路径path。 6. 若无法找到增广路径,则退出循环。 7. 计算增广路径上的最小残余容量minCapacity。 8. 更新残余容量矩阵residualCapacity: - 对于路径上的每一条边(u, v),将residualCapacity[u][v]减去minCapacity。 - 对于路径上的每一条反向边(v, u),将residualCapacity[v][u]增加minCapacity。 9. 将maxFlow增加minCapacity。 10. 输出maxFlow作为最大量。 在Matlab中实现Edmonds-Karp算法,可以根据以下步骤进行编程: - 首先,使用邻接矩阵adjMatrix表示图的连接关系,可以使用Matlab中的二维数组进行表示。 - 然后,根据输入的source和target,使用BFS算法找到增广路径。可以使用队列(Queue)数据结构实现BFS。 - 计算增广路径上的最小残余容量minCapacity,可以使用带权图中的最小值函数进行计算。 - 更新残余容量矩阵residualCapacity,可以根据路径上的边进行更新操作。 - 循环直到无法找到增广路径,可以使用while循环实现。 - 输出最大量maxFlow作为结果。 总之,通过以上步骤进行编程,即可在Matlab中实现Edmonds-Karp算法来解决最大问题。 ### 回答3: Edmonds-Karp算法是一种用于解决最大问题的算法,它是基于Ford-Fulkerson算法的改进版本。该算法通过在残余网络中不断寻找增广路径来找到最大。 在Matlab中实现Edmonds-Karp算法,可以按以下步骤进行: 1. 创建一个图表示残余网络,图中的边包含它们的容量和量。可以使用Matlab中的图论相关函数(如graph函数)创建图。 2. 初始化每条边的量为0。 3. 在残余网络中寻找增广路径,可以使用广度优先搜索(BFS)算法。从源节点开始,通过检查当前节点的邻居节点来找到增广路径。如果在残余网络中找到一条增广路径,就可以通过调整路径上的边的量来增加总量。 4. 不断重复步骤3,直到在残余网络中找不到增广路径为止,即没有可以进一步增加总量的路径。 5. 返回最大。 在Matlab中实现Edmonds-Karp算法需要对图的表示和广度优先搜索算法有一定的了解。可以使用Matlab中的图论包(如Graph and Digraph)来操作图,并使用队列数据结构来实现BFS算法。 需要注意的是,由于Edmonds-Karp算法的复杂度较高,如果图的规模比较大,算法可能需要较长的运行时间。在实际应用中,可以考虑使用其他更高效的最大算法,如Dinic算法或Push-Relabel算法
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值