最大流dinic算法

该算法的主要思想就是每次bfs判断是否可以到达目的地,并把路径记录下来,然后通过dfs进行多次增广,直到找不到增广路径就停止。

按道理来说时间复杂度应该降低了,为什么我的时间复杂度反而变高了。不太明白,,

具体代码:

#include <iostream>
#include <cstring>
#include <string>
#include <queue>
using namespace std;

const int N =10000+10;
const int M = 2e5+10;

int h[N],to[M],ne[M],cap[M],idx=0;
int dis[N],cur[N];

void add(int a,int b,int c){
    ne[idx]=h[a];
    cap[idx]=c;
    to[idx]=b;
    h[a]=idx++;
}

int bfs(int st,int ed){
    memset(dis,-1,sizeof dis);
    queue<int> q;
    q.push(st);
    dis[st]=0;
    cur[st]=h[st];
    while(q.size()){
        int t = q.front();
        q.pop();
        for(int i=h[t];~i;i=ne[i]){
            int v = to[i];
            if(cap[i]&&dis[v]==-1){
                dis[v]=dis[t]+1;
                cur[v]=h[v];
                if(v==ed)return 1;
                q.push(v);
            }
        }
    }
    return 0;
}


int dfs(int idx,int ed,int f){
    if(idx==ed||f==0)return f;
    int f_out=0;
    for(int i=cur[idx];~i;i=ne[i]){
        cur[idx]=i;
        int v = to[i];
        if(cap[i]&&dis[v]==dis[idx]+1){
            int d = dfs(v,ed,min(f-f_out,cap[i]));
            if(!d)dis[v]=-1;
            cap[i]-=d;
            cap[i^1]+=d;
            f_out+=d;
            if(f_out==f)return f;
        }
    }
    return f_out;
}

int mcmf(int st,int ed){
    int res=0,f;
    while(bfs(st,ed)){
        while(f=dfs(st,ed,0x3f3f3f3f)){
            res+=f;
        }
    }
    return res;
}

int main(){
    int n,m,s,t;
    cin>>n>>m>>s>>t;
    memset(h,-1,sizeof h);
    for(int i=0;i<m;i++){
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
        add(b,a,0);
    }
    cout<<mcmf(s,t)<<endl;
}

 

posted @ 2020-07-20 12:27  kstranger  阅读(90)  评论(0)    收藏  举报