网络流----EK算法求最大流

~~~
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1010,M=20010,INF=0x3f3f3f3f;
int h[N],e[M],ne[M],f[M],idx;//邻接表 头节点 边 next指针 容量
int q[N],pre[N],d[N];//bfs队列 前驱边 当前最大流量
bool st[N];
int n,m,S,T;

void add(int a,int b,int c)
{
    e[idx]=b,ne[idx]=h[a],f[idx]=c,h[a]=idx++;
    e[idx]=a,ne[idx]=h[b],f[idx]=0,h[b]=idx++;
}

bool bfs()
{
    memset(st,0,sizeof st);
    int hh=0,tt=0;
    q[0]=S,st[S]=true,d[S]=INF;
    
    while(hh<=tt)
    {
        auto t=q[hh++];
        for(int i=h[t];i!=-1;i=ne[i])
        {
            int ver=e[i];
            if(!st[ver]&&f[i])
            {
                st[ver]=true;
                pre[ver]=i;
                d[ver]=min(d[t],f[i]);
                if(ver==T) return true;
                q[++tt]=ver;
            }
        }
    }
    return false;
}

int EK()
{
    int res=0;
    while(bfs())
    {
        res+=d[T];
        for(int i=T;i!=S;i=e[pre[i]^1])
        {
            f[pre[i]]-=d[T];
            f[pre[i]^1]+=d[T];
        }
    }
    return res;
}

int main()
{
    cin>>n>>m>>S>>T;
    memset(h,-1,sizeof h);
    
    while(m--)
    {
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
    }
    
    cout<<EK()<<endl;
    
    return 0;
}
~~~