【模板】dinic算法

dinic算法用于解决最大流问题。

注意每次BFS之前把dist数组清空，源点的dist设为1。

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define inf 1000000000
using namespace std;
int ans,tot,vert,edg,S,T,fr[100005],to[200005],nxt[200005],f[200005];
int que[100005],dist[100005];
bool inq[100005];
void adde(int p,int q,int F)
{to[tot]=q;f[tot]=F;nxt[tot]=fr[p];fr[p]=tot++;}
bool bfs();
int dfs(int,int);
int main()
{
    memset(fr,-1,sizeof(fr));
    scanf("%d%d%d%d",&vert,&edg,&S,&T);
    int p,q,w,i;
    for(i=1;i<=edg;i++)
    {
        scanf("%d%d%d",&p,&q,&w);
        adde(p,q,w); adde(q,p,0);
    }
    while(bfs()) ans += dfs(S,inf);
    printf("%d",ans);
    return 0;
 } 
 bool bfs()
 {
     memset(dist,0,sizeof(dist));
     int l=0,r=0,i,t,now;
     que[++r] = S; dist[S] = 1;
     while(l<r)
     {
         now = que[++l];
         if(now==T) return true;
         for(i=fr[now];i!=-1;i=nxt[i])
         {
             t = to[i];
             if(!dist[t]&&f[i])
             {
                 dist[t] = dist[now]+1;
                 que[++r] = t;
             }
        }
    }
    return false;
 }
 int dfs(int now,int maxx)
 {
     if(now==T) return maxx;
     int i,t,ret=0,ss;
     for(i=fr[now];i!=-1;i=nxt[i])
     {
         t = to[i];
         if(dist[t]==dist[now]+1&&f[i])
         {
             ss = dfs(t,min(maxx-ret,f[i]));
             ret+=ss;
             f[i] -= ss; f[i^1] += ss;
             if(ret==maxx) return maxx;
        }
    }
    return ret;
 }