笔记 图论 网络流 EK算法

原题

#include<iostream>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=1010,M=200010,inf=1e8;//边的数量是10w所以残留网络边为20W
int n,m,S,T;//点 边 源点 汇点
int h[N],e[M],f[M],ne[M],idx;
//h表示点 e表示边 f表示容量 ne表示链 idx当前的点
int q[N],d[N],pre[N];
//q是用来宽搜的队列  d表示从起点走到每个点时，当前路径上所有边的容量最小值
//pre 记录路径的前驱点
bool st[N];
//用来bfs判重
void add(int a,int b,int c){
	e[idx]=b,f[idx]=c,ne[idx]=h[a],h[a]=idx++;
	e[idx]=a,f[idx]=0,ne[idx]=h[b],h[b]=idx++;//建立反向边
}
bool bfs(){
	int hh=0,tt=0;//宽搜的队头
	memset(st,false,sizeof st);
	q[0]=S,st[S]=true;d[S]=inf;
	while(hh<=tt){
		int t=q[hh++];
		for(int i=h[t];~i;i=ne[i]){
			int ver=e[i];//当前点的编号
			if(!st[ver]&&f[i]){//当前容量大于0
				st[ver]=true;
				//这里有点像最短路
				d[ver]=min(d[t],f[i]);//上一个点的容量和当前容量的最小值
				pre[ver]=i;
				if(ver==T)return T;//已经找到了终点
				q[++tt]=ver;//双指针维护队列区间
			}
		}
	}
	return false;
}
int EK(){
	int r=0;
	while(bfs()){//找增广路
		r+=d[T];//d[T]表示残余网络中起点到终点路径上的最小值
		for(int i=T;i!=S;i=e[pre[i]^1]){
			f[pre[i]]-=d[T];//计算去掉这次增广后路径上剩余容量
			f[pre[i]^1]+=d[T];//反向边要加上这个量
		}
	}
	return r;
}
int main(){
	cin>>n>>m>>S>>T;
	memset(h,-1,sizeof h);//清空邻接表
	while(m--){
		int a,b,c;//每个边的端点和容量
		cin>>a>>b>>c;
		add(a,b,c);
	}
	cout<<EK()<<endl;
}