网络流（二）

本篇博客主要内容是最大流模板的代码实现

1.Dinic算法的思想

这里说明一下，不写EK的原因是EK求最大流与Dinic相比要慢一些，所以在这里只写了Dinic我不会EKqwq

首先我们建立残留网络，在从中找到所有的增广路径，然后把所有增广路经所能流的流量加起来，即为最大流

可是现在有一个问题，就是如果流网络内存在一个边权均为正数的环，那么就会有无数条增广路径，会造成代码的死循环

所以我们引入分层图的概念

• 分层图

分层图实际上是从源点开始，每扩展一次就把层数加一（源点层数是0）

然后我们用dfs求出所有的增广路径

• 当前弧优化

该优化指的是我们在dfs每一个点时，我们不需要每次都从第一个点开始遍历，标记一下遍历到哪里了，然后为下一次减少时间

贴一下代码

#include <bits/stdc++.h>
using namespace std;
const int N=10010,M=200010,INF=1e8;
int n,m,S,T;
int h[N],e[M],ne[M],f[M],idx;
int q[N],d[N],cur[N];
void add(int a,int b,int c)
{
    e[idx]=b;f[idx]=c;ne[idx]=h[a];h[a]=idx++;
    e[idx]=a;f[idx]=0;ne[idx]=h[b];h[b]=idx++;
}
bool bfs()
{
    int hh=0,tt=-1;
    memset(d,-1,sizeof d);
    q[++tt]=S;d[S]=0;cur[S]=h[S];
    while(hh<=tt)
    {
        int t=q[hh++];
        for(int i=h[t];i!=-1;i=ne[i])
        {
            int ver=e[i];
            if(d[ver]==-1&&f[i])
            {
                d[ver]=d[t]+1;
                cur[ver]=h[ver];
                if(ver==T)return true;
                q[++tt]=ver;
            }
        }
    }
    return false;
}
int find(int u,int limit)
{
    if(u==T)return limit;
    int flow=0;
    for(int i=cur[u];i!=-1&&flow<limit;i=ne[i])
    {
        cur[u]=i;
        int ver=e[i];
        if(d[ver]==d[u]+1&&f[i])
        {
            int t=find(ver,min(f[i],limit-flow));
            if(!t)d[ver]=-1;
            f[i]-=t;f[i^1]+=t;flow+=t;
        }
    }
    return flow;
}
int dinic()
{
    int res=0,sum;
    while(bfs())while(sum=find(S,INF))res+=sum;
    return res;
}
int main()
{
    cin>>n>>m>>S>>T;
    memset(h,-1,sizeof  h);
    while(m--)
    {
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
    }
    cout<<dinic()<<endl;
    return 0;
}