Dinic算法解最大流问题

Dinic算法解最大流问题

具体的最大流的定义及原理可以参见这位朋友的博客(简单易懂)

最大流问题详解(侵删)

Dinic算法的基本原理

利用BFS对图进行分层处理，接着使用DFS从S开始，每次层次加一寻找下一个点，直至达到汇点T，然后再回溯回去，以此遍历所有的增广路径，这样就可以满足我们同时求出多条增广路的需求，效率更高

因此可给出Dinic算法的进行框架

①在残量网络上通过BFS求出结点的层次，构造分成图。

②在分层图上通过DFS寻找增广路，在回溯时同时更新边权。

更加具体的原理可以参见这位朋友的博客

浅谈网络最大流 (侵删)

代码部分

下面给出带有注释的代码（代码为借鉴上述博客中代码并增添详细解释）

#include <bits/stdc++.h>
using namespace std;
const long long inf=2225020666;
int n,m,s,t,u,v;
long long w,ans,dis[1005000];
int tot = 1, now[1005000], head[1005000]; 
/*dis数组用以设置每个点的层数，head数组用以构造图，now数组用以当前弧优化*/
struct node {
	int to, net;
	long long val;
} e[1005000];

void add(int u,int v,long long w) { //将点和权值加入图中
	e[++tot].to=v;
	e[tot].val=w;
	e[tot].net=head[u];
	head[u]=tot;
}

int bfs() {  //在残量网络中构造分层图 
	for(int i=1;i<=n;i++) 
		dis[i]=inf;
	queue<int> q;
	q.push(s); //源点入队
	dis[s]=0; //源点为第0层
	now[s]=head[s]; //弧优化
	while(!q.empty()) {
		int x=q.front(); //取出队首点
		q.pop();
		for(int i=head[x];i;i=e[i].net) {
			int v=e[i].to; //所连接的下一个点
			if(e[i].val>0 && dis[v]==inf) { //权值大于0且未遍历过
				q.push(v); //入队
				now[v]=head[v];
				dis[v]=dis[x]+1; //层数为上层加一
				if(v==t) 
					return 1; // 到了源点返回1
			}
		}
	}
	return 0; //队空返回0，结束DFS遍历
}

int dfs(int x,long long sum) {  //sum是整条增广路对最大流的贡献
	if(x==t) 
		return sum;
	long long k,res=0;  //k是当前最小的剩余容量 
	for(int i=now[x];i&&sum;i=e[i].net) {
		now[x]=i;  //当前弧优化 
		int v=e[i].to;
		if(e[i].val>0&&(dis[v]==dis[x]+1)) {
			k=dfs(v,min(sum,e[i].val)); // DFS下一个结点
			if(k==0) dis[v]=inf;  //剪枝，去掉增广完毕的点 
			e[i].val-=k;
			e[i^1].val+=k; //i^1为奇数-1，偶数+1
			res+=k;  //res表示经过该点的所有流量和（相当于流出的总量） 
			sum-=k;  //sum表示经过该点的剩余流量 
		}
	}
	return res;
}

int main() {
	scanf("%d%d%d%d",&n,&m,&s,&t);
	for(register int i=1;i<=m;i++) {
		scanf("%d%d%lld",&u,&v,&w);
		add(u,v,w);
		add(v,u,0); 
	}
	while(bfs()) { // 注意这块是每广搜一次进行一次深搜
		ans+=dfs(s,inf);  //流量守恒（流入=流出） 
	}
	printf("%lld",ans);
	return 0;
}