网络流的概念
网络流里有一些很重要的概念:
可行流
剩余容量
残余网络
增广路
求最大流的算法
所有求最大流的算法的本质都是:当一个可行流的残余网络中不存在可行流(即没有增广路)时,这个可行流是最大流。
EK求最大流
这个算法用的比较少,在此不展示代码。
Dinic求最大流
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
const int N=1e4+5,M=2e5+5,INF=0x3f3f3f3f;
int n,m,s,t;
int head[N],to[M],ne[M],w[M],idx=1;
void add(int x,int y,int z){
ne[++idx]=head[x];to[idx]=y;
w[idx]=z;head[x]=idx;
}
int d[N];//分层图
bool bfs(){
memset(d,0,sizeof d);
queue<int> q;
d[s]=1;q.push(s);
while(q.size()){
int u=q.front();
q.pop();
for(int i=head[u];i;i=ne[i]){
int j=to[i];
if(w[i] && !d[j]){
d[j]=d[u]+1;
q.push(j);
if(j==t) return true;//前提:t必须经过增广路遍历到
}
}
}
return false;
}
int dinic(int u,int flow){//flow实际上是剩余容量
if(u==t) return flow;
int rest=flow;
for(int i=head[u];i && rest;i=ne[i]){
int j=to[i];
if(w[i] && d[j]==d[u]+1){
int k=dinic(j,min(rest,w[i]));//返回的k是变化的剩余容量
if(!k) d[j]=0;
w[i]-=k;
w[i^1]+=k;
rest-=k;
}
}
return flow-rest;
}
int main(){
cin>>n>>m>>s>>t;
while(m--){
int x,y,z;
cin>>x>>y>>z;
add(x,y,z);add(y,x,0);
}
int flow=0,maxflow=0;
while(bfs()) while(flow=dinic(s,INF)) maxflow+=flow;
cout<<maxflow<<endl;
return 0;
}