首页 > 代码库 > 最大流之dinic
最大流之dinic
先用bfs预处理出层次图,然后在层次图上用dfs找增广路径,理论复杂度O(n*n*m)
const int INF=0xfffffff ; struct node{ int s,t,cap,nxt ; }e[200005] ; int m,n,head[10005],level[10005],cnt ; void add(int s,int t,int cap) { e[cnt].s=s ;e[cnt].t=t ;e[cnt].cap=cap ;e[cnt].nxt=head[s] ;head[s]=cnt++ ; e[cnt].s=t ;e[cnt].t=s ;e[cnt].cap=0 ;e[cnt].nxt=head[t] ;head[t]=cnt++ ; } bool build(int s,int t) { memset(level,-1,sizeof(level)) ; queue <int> q ; level[s]=0 ; q.push(s) ; while(!q.empty()) { int u=q.front() ; q.pop() ; for(int i=head[u] ;i!=-1 ;i=e[i].nxt) { int tt=e[i].t ; if(level[tt]==-1 && e[i].cap>0) { level[tt]=level[u]+1 ; if(tt==t)return true ; q.push(tt) ; } } } return false ; } int find(int s,int t,int flow) { if(s==t)return flow ; int ret=0,a ; for(int i=head[s] ;i!=-1 ;i=e[i].nxt) { int tt=e[i].t ; if(level[tt]==level[s]+1 && e[i].cap>0 && (a=find(tt,t,min(flow,e[i].cap)))) { e[i].cap-=a ; e[i^1].cap+=a ; return a ; } } return 0 ; } int dinic(int s,int t) { int flow,ret=0 ; while(build(s,t)) while(flow=find(s,t,INF)) ret+=flow ; return ret ; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。