最大流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; }

浙公网安备 33010602011771号