最大流算法 Dinic算法 ADT

https://www.cnblogs.com/SYCstudio/p/7260613.html

这个博主写的很好完全看得懂就不自己写了

方便自己使用还是把模板抄出来吧。

ADT:

const int max_e = 55000;
const int max_v = 250;
const int INF = 0x3f3f3f3f; //可能有时候要用Long Long

struct Edge{
    int from,to,cap;  //很多时候要用long long int cap，因为可能不够
    int next; //也可以用vector来存 从from出来的下一条边
}edges[max_e];
int head[max_v]; //第一个从这个点出来的边对应的序号
int cnt; //边的数量
int Depth[max_v];
int f[max_e];
int S,T;

void init(){
    cnt = 0;
    memset(head,-1,sizeof(head));
    memset(f,0,sizeof(f));
}

void AddEdge(int f, int t, int cap){
    edges[cnt].from = f;
    edges[cnt].to = t;
    edges[cnt].cap = cap;
    edges[cnt].next = head[f];
    head[f] = cnt;
    cnt++;
    //？？？？？？？？？？
    edges[cnt].from = t;
    edges[cnt].to = f;
    edges[cnt].cap = 0; //反向边
    edges[cnt].next = head[t];
    head[t] = cnt;
    cnt++;
    
}

bool bfs(){
    queue<int>Q;
    while(!Q.empty()) Q.pop();
    memset(Depth,0,sizeof(Depth));
    Depth[S] = 1;
    Q.push(S);
    while(!Q.empty()){
        int u = Q.front();
        Q.pop();
        for(int i=head[u];i!=-1;i=edges[i].next){
            int t = edges[i].to;
            if((Depth[t]==0)&&edges[i].cap>f[i]){
                Depth[t] = Depth[u] + 1;
                Q.push(t);
            }
        }
    }
    if(Depth[T]>0){
        return 1;
    }
    return 0;
}

int dfs(int u, int dist){ //当前节点，当前流量
    if(u==T) return dist;
    for(int i=head[u];i!=-1&&dist;i=edges[i].next){
        int t = edges[i].to;
        if(Depth[t]==(Depth[u]+1)&&edges[i].cap>f[i]){
            int di = dfs(t,min(dist,edges[i].cap-f[i]));
            if(di>0){//增广成功
                f[i] += di;
                f[i^1] -= di;
                return di;
            }
        }
    }
    
    return 0;
}

int dinic(){
    int ans=0;
    while(bfs()){
        while(int d=dfs(S,INF)) {
            ans += d;
        }
    }
    return ans;
}