Dinic算法

Dinic算法

Dinic算法用于求解最大网络流，改进了EK算法。

---

EK算法求增广路的时候一次bfs只是求了一个增光路径，然后再在增广路径上进行$flow[x][y]+=w,flow[y][x]-=w$操作。但是一次只寻找一条增广路径速度有点慢，我们考虑一次寻找多条增广路径。从而诞生了Dinic算法！

为什么一次只考虑一条增广路径有些不好呢？我们看看已下这个经典的例子

这个图片中，我们直到一条可能的路径是从$S->V->U->t$，进行$flow[x][y]+=w,flow[y][x]-=w$操作后（此处的$x,y$代表增广路中的每一条边的两个节点），下一次寻找的增广路径是$S->u->v->t$，对于$v$和$u$点的$flow[x][y]+=w,flow[y][x]-=w$操作会抵消上一次的操作，出现这个问题是因为我们每次只是寻找一条增广路径造成的，试想一下如果我们第一次从$S$走到$V$之后,让$V$沿着两条路径走到$t$，也就是$V->t$和$V->u->t$，那么这个问题是不是很好的就解决了？而且一次可以找到多个多条增广路。下面来看下代码！（若代码有任何问题请指出！但思路是对的^ - ^）

#include<bits/stdc++.h>
using namespace std;
const int maxn = 100;
int level[maxn];
int flow[maxn][maxn];
int u,v;
int head[maxn];
int cnt=0;
int vis[maxn];
bool bfs(int u);
int Dinic(int u,int t);
int dfs(int u,int w);
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    int x,y,w;
    for(int i=0;i<m;i++){
        cin>>x>>y>>w;
        flow[x][y]=w;
    }
    return 0;
}
int Dinic(int u,int v){
    int ans=0;
    while(bfs(u)!=-1){
        ans+=dfs(u,v);
    }
    return ans;
}
int dfs(int u, int w){
    if(u==w)
        return w;
    int temp=0;
    for(int i=0;i<maxn;i++){
        if(level[i]=level[u]+1){
            int x = dfs(i,min(w,flow[u][i]));
            flow[u][i]-=x;
            flow[i][u]+=x;
            temp+=x;
            w-=x;
            if(w<=0)
                break;
        }
    }
    return temp;
}
bool bfs(int u){
    memset(level,-1,sizeof(level));
    memset(vis,0,sizeof(vis));
    queue<int>q;
    q.push(u);
    level[u]=0;
    vis[u]=1;
    while(!q.empty()){
        int t = q.front();
        q.pop();
        for(int i=0;i<maxn;i++){
            if(flow[t][i]&&!vis[i]){
                q.push(i);
                vis[i]=1;
                level[i]=level[t]+1;
            }
        }
    }
    if(level[v]==-1)
        return false;
    else
        return true;
}

我们一步步的来看，我们每次都bfs网络流，设置已下每个节点的level值，个人认为这一步不设置level也行，设置level应该只是为了等下dfs的时候防止往回计算（也就是说我们希望总的趋势是从S原点一直计算到t汇点的，但是我们一会儿再更新flow[][]的时候我们不仅仅会让正向减小$flow[x][y]-=x$，也会让反向增加$flow[y][x]+=w$，这样一来我们便不知道了我们图原来拓扑的顺序，（若从s源点放少量的水的话，水将摸不清应该怎么流才能流到汇点），注意的是bfs的时候需要判断下$flow[t][i]$，因为bfs是帮助寻找增广路的，如果$flow[][]$是0了那显然不是增广路径。

再来看看dfs做什么，我们每次都是bfs一下然后dfs一下，bfs刚才说了是标记level的，int dfs(int u,int w)，其中u代表当前处理节点，w代表当前处理节点允许通过的最大流（注意只是允许，能不能通过这么大的流还不一定呢，实际上真正通过的流是return回来的那个值）。dfs要做的事情就是找出u节点相邻的节点能通过的流量的总和，然后然后该流量，再用公式解释下：若$level[i]=level[u]+1$也就是说i是t的下一层节点，对所有的i,执行dfs操作，dfs操作返回的是能经过u的流量的做大值，那怎么求能流过i的流量的最大值呢？继续对i的下一层进行dfs啊。我们现在执行对i节点执行dfs操作，注意，我们传进去的w是min(w,flow[u][i]),表示i节点允许通过的最大流，如果w<flow[u][i]那么显然应该选择w，即使你的这条边是flow[u][i]，但是上一级能通过的最大只有w，所以你只能最大通过w。然后是根据返回的int x = dfs(i,min(w,flow[u][i]))更新flow，最后temp相信不难理解，temp就是我们真正能通过的最大流。还有w-=x，因为我们每一次dfs之后我们还需要继续对下一个节点进行dfs，这时允许通过这个节点的最大流减小了x，因为减小的x被上一个i拿走了！