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&∑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;
}