最大流算法
最大流算法(EK算法)
今天学习了最大流,最大流算法主要有两个,一个是EK算法,一个是dinic算法。EK算法的时间复杂度较大,但容易理解,相对的,dinic算法的时间复杂度较小,但却难以理解。由于我是个蒟蒻,dinic算法的实现我看不懂,所以我只学习了最大流的EK算法(doge)。
EK算法的主要思想是在当前图中找增广路,然后再沿着这条增广路去更新流量,直到找不到新的增广路为止。经过这样的操作后,最大流就是每次更新的流量的总和。在这里要引入一个东西,叫做反向边,反向边的作用详见下面大佬的博客。
大佬的博客
#include <iostream>
#include <cstring>
#include <queue>
#include <algorithm>
using namespace std;
#define N 205
#define INF 0x7fffffff
int n,m;
int pre[N],flag[N],mp[N][N],flow[N][N];
bool bfs(int st,int ed)
{
memset(pre,0,sizeof(pre));
memset(flag,0,sizeof(flag));
queue <int> q;
flag[st]=1;
q.push(st);
while(!q.empty())
{
int now=q.front();
q.pop();
for(int i=1;i<=m;i++)
{
if(i!=now&&!flag[i]&&mp[now][i]-flow[now][i]>0)
{
pre[i]=now;
flag[i]=1;
q.push(i);
if(i==ed)
return true;
}
}
}
return false;
}
int EK(int st,int ed)
{
int ans=0;
while(bfs(st,ed))
{
int delta=INF;
for(int i=ed;i!=st;i=pre[i])
delta=min(delta,mp[pre[i]][i]-flow[pre[i]][i]);
for(int i=ed;i!=st;i=pre[i])
{
flow[pre[i]][i]+=delta;
flow[i][pre[i]]-=delta;
}
ans+=delta;
}
return ans;
}
int main()
{
cin >> n >> m;
while(n--)
{
int s,e,c;
cin >> s >> e >> c;
mp[s][e]+=c;
mp[e][s]+=c;
flow[s][e]=0;
flow[e][s]+=c;
}
cout << EK(1,m) << endl;
return 0;
}
真就板子题,没啥好说的……