网络流初步

网络流的概念

网络流里有一些很重要的概念：
可行流
剩余容量
残余网络
增广路

求最大流的算法

所有求最大流的算法的本质都是：当一个可行流的残余网络中不存在可行流（即没有增广路）时，这个可行流是最大流。

EK求最大流

这个算法用的比较少，在此不展示代码。

Dinic求最大流

#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
const int N=1e4+5,M=2e5+5,INF=0x3f3f3f3f;
int n,m,s,t;
int head[N],to[M],ne[M],w[M],idx=1;
void add(int x,int y,int z){
    ne[++idx]=head[x];to[idx]=y;
    w[idx]=z;head[x]=idx;
}
int d[N];//分层图
bool bfs(){
    memset(d,0,sizeof d);
    queue<int> q;
    d[s]=1;q.push(s);
    while(q.size()){
        int u=q.front();
        q.pop();
        for(int i=head[u];i;i=ne[i]){
            int j=to[i];
            if(w[i] && !d[j]){
                d[j]=d[u]+1;
                q.push(j);
                if(j==t) return true;//前提：t必须经过增广路遍历到
            }
        }
    }
    return false;
}
int dinic(int u,int flow){//flow实际上是剩余容量
    if(u==t) return flow;
    int rest=flow;
    for(int i=head[u];i && rest;i=ne[i]){
        int j=to[i];
        if(w[i] && d[j]==d[u]+1){
            int k=dinic(j,min(rest,w[i]));//返回的k是变化的剩余容量
            if(!k) d[j]=0;
            w[i]-=k;
            w[i^1]+=k;
            rest-=k;
        }
    }
    return flow-rest;
}
int main(){
    cin>>n>>m>>s>>t;
    while(m--){
        int x,y,z;
        cin>>x>>y>>z;
        add(x,y,z);add(y,x,0);
    }
    int flow=0,maxflow=0;
    while(bfs()) while(flow=dinic(s,INF)) maxflow+=flow;
    cout<<maxflow<<endl;
    return 0;
}

最大流之二分图匹配