上篇记录了一下EK,这篇就记录一下我又从头再学的Dinic,(据说是目前最好用的网络流算法之一 (也可能没有之一) )。
推荐博文:https://baijiahao.baidu.com/s?id=1612179096991409044&wfr=spider&for=pc
Dinic
如果你搞明白了EK算法的话,其实Dinic就也不是那么的难以理解了。首先我们可以回想一下EK的思想,一句话就是不断bfs找增广路,直到找不到为止。当然Dinic在最大流中也是不断的找增广路,但是和EK的区别是它先把图分层,然后一下多次增广。
bfs分层
我们先把整个图按照BFS遍历顺序分出不同的层次,分层的原因是(肯定是为了后面DFS增广做准备了)我们分层之后我们可以保证这条路上可以走通并且目前是最短的增广路,(我现在理解就是先走短的顺的再走长的曲折的路,因为有反向边的存在,然后不断更新边权之后,可能会最后会有曲曲折折的通路)。
我们一直分层直到最后到达不了T(终点)结束。
dfs多次增广
这找增广路我们按照分层的层次去找,Dinic比EK优秀的地方就是这,因为我们按层次去找,找不到就回溯,但是这个分层还可以用,我们可以在这个层次上找到这种层次上的所以增广路,也就是一次分层多次增广。因为DFS特点我们不用记录路径,效率也很快。
来个板子试试手:P3376 【模板】网络最大流
完整代码
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<algorithm>
#include<queue>
#include<vector>
#include<math.h>
#include<bits/stdc++.h>
#include<map>
using namespace std;
#define LL long long
#define inf 0x3f3f3f3f
double eps=1e-10;
const LL N=1e5+10;
int n,m,s,t;
struct zxc
{
int to,net,w;
}e[N*2];
int tot=0;
int head[N];
int dep[N];
int vis[N];
int maxflow=0;
void add(int u,int v,int w)
{
e[tot]={v,head[u],w};
head[u]=tot++;
e[tot]={u,head[v],0};
head[v]=tot++;
}
void init()
{
memset(head,-1,sizeof(head));
tot=0;
maxflow=0;
}
int bfs()//分层
{
queue<int>q;
memset(vis,0,sizeof(vis));
memset(dep,0x3f3f3f3f,sizeof(dep));
q.push(s);
vis[s]=1;
dep[s]=0;
while(!q.empty())
{
int u=q.front();
q.pop();
vis[u]=0;
for(int i=head[u];i!=-1;i=e[i].net)
{
int y=e[i].to;
if(dep[y]>dep[u]+1&&e[i].w>0)
{
dep[y]=dep[u]+1;
if(!vis[y])
{
vis[y]=1;
q.push(y);
}
}
}
}
if(dep[t]==0x3f3f3f3f)
{
return 0;
}
else
{
return 1;
}
}
int dfs(int x,int flow)//多次增广
{
int row=0;
if(x==t)return flow;
for(int i=head[x];i!=-1;i=e[i].net)
{
int y=e[i].to;
if(dep[y]==dep[x]+1&&e[i].w>0)
{
if(row=dfs(y,min(flow,e[i].w)))
{
e[i].w-=row;
e[i^1].w+=row;
return row;
}
}
}
return 0;
}
void Dinic()
{
while(bfs())
{
int flow;
while(flow=dfs(s,0x3f3f3f3f))
{
maxflow+=flow;
}
}
}
int main()
{
init();
scanf("%d%d%d%d",&n,&m,&s,&t);//n个点,m条边,起点S,终点T
int u,v,w;
for(int i=1;i<=m;i++)
{
scanf("%d%d%d",&u,&v,&w);
add(u,v,w);
}
Dinic();
printf("%d\n",maxflow);
return 0;
}