Dinic——模板整理

Dinic：最大流算法
时间复杂度：$O(n^2m)$

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn=10005,maxe=200005;
struct Edge { int from,to,cap,flow; } Es[maxe];
int n,m,d[maxn],now[maxn],nxt[maxe],S,T,fir[maxn],pos[maxn],tot=1;
queue <int> que;
void add(int x,int y,int w){
    Es[++tot]=(Edge){x,y,w,0};
    nxt[tot]=fir[x]; fir[x]=tot;
    Es[++tot]=(Edge){y,x,0,0};
    nxt[tot]=fir[y]; fir[y]=tot;
}
bool BFS(){
    memset(d,63,sizeof(d)); int INF=d[0];
    while(!que.empty()) que.pop(); 
    que.push(S); d[S]=0; 
    while(!que.empty()){
        int x=que.front(); que.pop();
        for(int j=fir[x];j;j=nxt[j]) if(d[Es[j].to]==INF&&Es[j].cap>Es[j].flow){
            d[Es[j].to]=d[x]+1; que.push(Es[j].to);
        }
    }
    return d[T]!=INF;
}
int find_DFS(int x,int flow){
    if(x==T||flow==0) return flow;
    int res=0,f;
    for(int &j=pos[x];j;j=nxt[j]){
        if(d[x]+1==d[Es[j].to]&&(f=find_DFS(Es[j].to,min(flow,Es[j].cap-Es[j].flow)))>0){
            Es[j].flow+=f; Es[j^1].flow-=f; 
            res+=f; flow-=f;
            if(flow==0) break;
        }
    }
    return res;
}
int GetMaxFlow() {
    int MaxFlow=0,now;
    while(BFS()){
        for(int i=1;i<=n;i++) pos[i]=fir[i]; 
        MaxFlow+=find_DFS(S,1e+9);
    }
    return MaxFlow;
}
int main(){
    freopen("Dinic.in","r",stdin);
    freopen("Dinic.out","w",stdout); 
    scanf("%d%d%d%d",&n,&m,&S,&T);
    for(int i=1;i<=m;i++){
        int x,y,z; scanf("%d%d%d",&x,&y,&z);
        add(x,y,z);
    }
    printf("%d\n",GetMaxFlow());
    return 0;
}