HihoCoder - 1369 网络流一 最大流 Edmonds-Karp

题目链接

题意：n个点，m条道路，每条道路有最大容量，求从点1到点n最大的流量。

思路：网络流最大流模板题。

网络流最大流：先了解一下网络流中的增广路。

增广路：一条从源点(出发的点)到汇点(到达的点)的路径，该路径上所有前向弧(与链的方向相同)为非饱和弧，后向弧(与链的方向相反)为非零弧。

则根据增广路定义可知，若图中不存在增广路，则当前的可行流即为最大流。因为若存在增广路，则该路上前向弧均为不饱和弧，即仍有可增加的流量空间。

所以我们只要不停求增广路直到不存在即可。常用的算法其中有一个Edmonds-Karp。

Edmonds-Karp：当一个可行流目前已经无法进行扩充，但也有可能并非答案，有可能存在通过某点i有不同的走法而导致有更优解。若是通过回溯来解决的话时间复杂度会大大增加，而该算法重要的就是引入了一个反向边的概念，就是与题目所给的对应边方向相反，这样每次对边修改的时候即可对其反向边进行相反的修改，在下一次修改反向边时就相当于对该边操作的回溯。之后我们就可以不停找增广路，取min(一条边的容量-该边当前流量)，不停的增广，构建新的网络，最后就能求的答案。然而若是无目的寻找增广路时间复杂度也是相当大的。通过优化，我们可以每次选择最短的增广路进行增广，所以我们应采取bfs而不是dfs求增广路。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define PI acos(-1)
#define INF 0x3f3f3f3f
#define NUM 510
#define debug true
#define ll long long
#define lowbit(x) ((-x)&x)
#define ffor(i,d,u) for(int i=d;i<=u;++i)
#define _ffor(i,u,d) for(int i=u;i>=d;--i)
#define mst(array,Num) memset(array,Num,sizeof(array))
const int p = 1e9+7;
int n,m;
int path[NUM][NUM];//path存该边当前剩余容量
int pre[NUM],flow[NUM];
template <typename T>
void read(T& x)
{
    x=0;
    char c;T t=1;
    while(((c=getchar())<'0'||c>'9')&&c!='-');
    if(c=='-'){t=-1;c=getchar();}
    do(x*=10)+=(int)(c-'0');while((c=getchar())>='0'&&c<='9');
    x*=t;
}
template <typename T>
void write(T x)
{
    int len=0;char c[21];
    if(x<0)putchar('-'),x*=(-1);
    do{++len;c[len]=(x%10)+'0';}while(x/=10);
    _ffor(i,len,1)putchar(c[i]);
}
inline int bfs()
{
    queue < int > q;
    int x;
    mst(pre,-1);
    pre[1]=0;
    q.push(1);
    flow[1]=INF;
    while(!q.empty())
    {
        x=q.front();
        q.pop();
        if(x==n)break;//一次求一条增广路
        ffor(i,1,n)
        {
            if(path[x][i]<=0||pre[i]!=-1)continue ;
            pre[i]=x;
            flow[i]=min(flow[x],path[x][i]);
            q.push(i);
        }
    }
    return (pre[n]==-1)?(-1):flow[n];
}
inline int solve()
{
    int change,x,ans=0;
    while((change=bfs())!=-1)
    {
        x=n;
        while(pre[x]!=0)
        {
            path[pre[x]][x]-=change;
            path[x][pre[x]]+=change;//反向边进行相反操作
            x=pre[x];
        }
        ans+=change;
    }
    return ans;
}
inline void AC()
{
    int x,y,w;
    read(n),read(m);
    mst(path,0);
    ffor(i,1,m)
    {
        read(x),read(y),read(w);
        path[x][y]+=w;
    }
    write(solve());
}
int main()
{
    AC();
    return 0;
}