Dinic算法模板

没什么好说的，建议直接背过。

“Dinic的复杂度就是个笑话，跟放P一样”

看似 (O(n^2m)) 实则艹过 (n=10^5,m=10^6)

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

const int N=1e4+10,M=2e5+10,INF=1e8;

int n,m,s,t;
int head[N],ver[M],nxt[M],cc[M],tot=0;
void add(int x,int y,int c)
{
    ver[tot]=y;cc[tot]=c;nxt[tot]=head[x];head[x]=tot++;
    ver[tot]=x;cc[tot]=0;nxt[tot]=head[y];head[y]=tot++;
}
int q[N],d[N],cur[N];//当前路优化

bool bfs()
{
    int hh,tt;
    hh=tt=0;
    memset(d,-1,sizeof d);
    q[0]=s,d[s]=0,cur[s]=head[s];
    while(hh<=tt)
    {
        int x=q[hh++];
        for(int i=head[x];~i;i=nxt[i])
        {
            int y=ver[i];
            if(d[y]==-1 &&cc[i])
            {
                d[y]=d[x]+1;
                cur[y]=head[y];//记录当前路
                if(y==t) return 1; 
                q[++tt]=y;
            }
        }
    }
    return 0;
}

int find(int u,int lim)//从源点流向u点的最大流量是lim的话
{
    if(u==t) return lim;
    int flow=0;
    for(int i=cur[u];~i && flow<lim;i=nxt[i])
    {
        cur[u]=i;//记录当前路
        int y=ver[i];
        if(d[y]==d[u]+1 &&cc[i])
        {
            int tmp=find(y,min(cc[i],lim-flow));
            if(!tmp) d[y]=-1;
            cc[i]-=tmp;cc[i^1]+=tmp;flow+=tmp;
        }
    }
    return flow;
}

int dinic()
{
    int res=0,flow;
    while(bfs())
    {
        while(flow=find(s,INF)) res+=flow;
    }
    return res;
}

int main()
{
    scanf("%d%d%d%d",&n,&m,&s,&t);
    memset(head,-1,sizeof head);
    for(int i=1;i<=m;i++)
    { 
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
    }
    printf("%d",dinic());
    return 0;
}