网络最大流Dinic算法

 

一、什么是最大流问题 （摘自网络，原文地址https://www.cnblogs.com/qiucz/p/4601241.html）

假设现在有一个地下水管道网络，有m根管道，n个管道交叉点，现在自来水厂位于其中一个点，向网络中输水，隔壁老王在另外一个点接水，已知由于管道修建的年代不同，有的管道能承受的水流量较大，有的较小，现在求在自来水厂输入的水不限的情况下，隔壁老王能接到的水的最大值？

为解决该问题，可以将输水网络抽象成一个联通的有向图，每根管道是一条边，交叉点为一个结点，从u流向v的管道能承受的最大流量称为容量，设为cap[u][v]，而该管道实际流过的流量设为flow[u][v]，自来水厂称为源点s，隔壁老王家称为汇点t，则该问题求的是最终流入汇点的总流量flow的最大值。

dinic算法就一句话：bfs分层，dfs增广

#include<bits/stdc++.h>
#define N 5005
#define M 10005
#define inf 999999999
using namespace std;

namespace io {
    const int SIZE = 1e7 + 10;
    char inbuff[SIZE];
    char *l, *r;
    inline void init() {
        l = inbuff;
        r = inbuff + fread(inbuff, 1, SIZE, stdin);
    }
    inline char gc() {
        if (l == r) init();
        return (l != r) ? *(l++) : EOF;
    }
    void read(int &x) {
        x = 0; char ch = gc();
        while(!isdigit(ch)) ch = gc();
        while(isdigit(ch)) x = x * 10 + ch - '0', ch = gc();
    }
} using io::read;

int n,m,s,t,num,adj[N],dis[N],q[N];
struct edge{
    int v,w,pre;
}e[M];

void addedge(int u,int v,int w){
    e[num]=(edge){v,w,adj[u]};
    adj[u]=num++;
    e[num]=(edge){u,0,adj[v]};
    adj[v]=num++;
}

int bfs(){
    int i,x,v,tail=0,head=0;
    memset(dis,0,sizeof(dis));
    dis[s]=1;
    q[tail++]=s;
    while(head<tail){
        x=q[head++];
        for(i=adj[x];i!=-1;i=e[i].pre)
        if(e[i].w&&dis[v=e[i].v]==0){
        dis[v]=dis[x]+1;
        if(v==t)
        return 1;
        q[tail++]=v;
        }
    }
    return 0;
}

int dfs(int s,int limit){
    if(s==t)
    return limit;
    int i,v,tmp,cost=0;
    for(i=adj[s];i!=-1;i=e[i].pre)
        if(e[i].w&&dis[s]==dis[v=e[i].v]-1){
            tmp=dfs(v,min(limit-cost,e[i].w));
            if(tmp>0){
                e[i].w-=tmp;
                e[i^1].w+=tmp;
                cost+=tmp;
                if(limit==cost)
                    break;
            }
            else dis[v]=-1;
        }
    return cost;

}

int Dinic(){
    int ans=0;
    while(bfs())
        ans+=dfs(s,inf);
    return ans;
}

int main (){
    read(n);read(m);
    int u,v,w;
    memset(adj,-1,sizeof(adj));
    num=0;
    s=1;
    t=n;
    while(m--){
        //scanf("%d%d%d",&u,&v,&w);
        read(u);read(v);read(w);
        addedge(u,v,w);
    }
    printf("%d\n",Dinic());
}