【最大流模板】

一、FF算法

1.先调用init（）初始化，用add（from，to，cap）建边。然后max_flow(s，t)最大流，s为源点，t为汇点。
2.N是点数，n是点数，m是边数。
3.cap容量，跑完最大流后cap为零即为割边。
const int N=500;
const int INF=0x3f3f3f3f;
struct edge
{
    int to,cap, rev;
    edge(int x=0,int y=0,int z=0)
    {
        to=x;cap=y;rev=z;
    }
};
vector<edge> g[N];
bool used[N];
int n,m,s,t;
void add(int from,int to,int cap)
{
    g[from].push_back(edge(to,cap,g[to].size()));
    g[to].push_back(edge(from,0,g[from].size()-1));
}
int dfs(int v,int t,int f)
{
    if(v==t)return f;
    used[v]=true;
    for(int i=0;i<g[v].size();i++)
    {
        edge &e=g[v][i];
        if(used[e.to]||e.cap<=0)continue;
        int d=dfs(e.to,t,min(e.cap,f));
        if(d>0)
        {
            e.cap-=d;
            g[e.to][e.rev].cap+=d;
            return d;
        }
    }
    return 0;
}
int max_flow(int s,int t)
{
    int flow=0;
    while(1)
    {
        memset(used,0,sizeof(used));
        int f=dfs(s,t,INF);
        if(f==0)return flow;
        flow+=f;
    }
}
void init()
{
    for(int i=1;i<=n;i++)
        g[i].clear();
}

二、Dinic算法（理论O（m*n^2））

1.先调用init（）初始化，用add（from，to，cap）建边。然后max_flow(s,t)求最大流，s为源点，t为汇点。
2.N是点数，n是实际点数，m是实际边数。
3.cap为容量，cap为0表示割边。
const int N=1000;
const int INF=0x3f3f3f3f;
struct edge
{
    int to,cap,rev;
    edge(int x,int y,int z)
    {
        to=x;cap=y;rev=z;
    }
};
vector<edge> g[N];
int level[N],iter[N];
int n,m,s,t;
bool vis[N];
void add(int from,int to,int cap)
{
    g[from].push_back(edge(to,cap,g[to].size()));
    g[to].push_back(edge(from,0,g[from].size()-1));
}
void bfs(int s)
{
    memset(level,-1,sizeof(level));
    level[s]=0;
    queue<int> que;
    que.push(s);
    while(!que.empty())
    {
        int v=que.front();que.pop();
        for(int i=0;i<g[v].size();i++)
        {
            edge &e=g[v][i];
            if(e.cap>0&&level[e.to]<0)
            {
                level[e.to]=level[v]+1;
                que.push(e.to);
            }
        }
    }
}
int dfs(int v,int t,int f)
{
    if(v==t)return f;
    for(int &i=iter[v];i<g[v].size();i++)
    {
        edge &e=g[v][i];
        if(e.cap>0&&level[e.to]>level[v])
        {
            int d=dfs(e.to,t,min(f,e.cap));
            if(d>0)
            {
                e.cap-=d;
                g[e.to][e.rev].cap+=d;
                return d;
            }
        }
    }
    return 0;
}
int max_flow(int s,int t)
{
    int  flow=0;
    while(1)
    {
        bfs(s);
        if(level[t]<0)return flow;
        memset(iter,0,sizeof(iter));
        int f;
        while(f=dfs(s,t,INF))
        {
            flow+=f;
        }
    }
}
void init()
{
    for(int i=1;i<=n;i++)g[i].clear();
}