最大流 Dinic算法

dinic算法的主体是利用DFS来求增广路。可以理解成从源节点倒无穷多的水下去，看看最后有多少能够流倒汇。配合当前弧优化，dinic算法的效率还是不错的。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 10006
#define maxe 200006
#define INF 1000000000
using namespace std;
inline char nc(){
    static char buf[100000],*i=buf,*j=buf;
    return i==j&&(j=(i=buf)+fread(buf,1,100000,stdin),i==j)?EOF:*i++;
}
inline int _read(){
    char ch=nc();int sum=0;
    while(!(ch>='0'&&ch<='9'))ch=nc();
    while(ch>='0'&&ch<='9')sum=sum*10+ch-48,ch=nc();
    return sum;
}
int n,e,S,T,tot,ans,cur[maxn],que[maxn],lnk[maxn],dep[maxn],son[maxe],nxt[maxe],f[maxe],c[maxe];
void add(int x,int y,int z){
    nxt[++tot]=lnk[x];son[tot]=y;c[tot]=z;lnk[x]=tot;
}
bool bfs(){
    memset(dep,0,sizeof(dep));
    int hed=0,tal=1;que[1]=S;dep[S]=1;
    while(hed!=tal) for(int j=lnk[que[++hed]];j;j=nxt[j]) if((!dep[son[j]])&&f[j]<c[j])
        dep[son[j]]=dep[que[hed]]+1,que[++tal]=son[j];
    return dep[T]>0;
}
int dfs(int x,int w){
    int sum=w;
    if(x==T)return w;
    for(int j=cur[x];j;j=nxt[j]) if(dep[son[j]]==dep[x]+1&&f[j]<c[j]){
        int t=dfs(son[j],min(c[j]-f[j],sum));
        f[j]+=t;f[j^1]-=t;
        sum-=t;
        if(!sum){cur[x]=j+1;break;}
    }
    return w-sum;
}
int main(){
    freopen("dinic.in","r",stdin);
    freopen("dinic.out","w",stdout);
    n=_read();e=_read();S=_read();T=_read();tot=1;//-1有问题，j等于0就停了 
    for(int i=1,x,y,z;i<=e;i++)x=_read(),y=_read(),z=_read(),add(x,y,z),add(y,x,0);
    while(bfs()){
        memcpy(cur,lnk,sizeof(lnk));
        ans+=dfs(S,INF);
    }
    printf("%d\n",ans);
    return 0;
}