P3376 【模板】网络最大流（dinic算法）

dinic 算法思路：

①：在残留网络中 bfs 出层次图。
②：在层次图中找出最大的增广流，并入答案。
③：重复①、②步骤，直到在残留网络中找不到 $s$ 到 $t$ 的路径。

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=10010;

int read(){
	int x=0,f=1;char ch=getchar();
	while(!isdigit(ch)){
		if(ch=='-') f=-1;
		ch=getchar();
	}
	while(isdigit(ch)) x=x*10+ch-'0',ch=getchar();
	return x*f;
}

int n,m,S,T,head[N],tot=1,cur[N],dep[N];
struct node{
	int to,nxt,c;
}edge[N];
void add(int x,int y,int c){
	edge[++tot].to=y;
	edge[tot].c=c;
	edge[tot].nxt=head[x];
	head[x]=tot;
}

bool bfs(){
	memset(dep,0,sizeof(dep));
	queue<int> q;q.push(S),dep[S]=1;
	while(!q.empty()){
		int t=q.front();q.pop();
		for(int i=head[t];i;i=edge[i].nxt){
			int y=edge[i].to;
			if(!dep[y]&&edge[i].c){
				dep[y]=dep[t]+1;
				q.push(y);
				if(y==T) return true;
			}
		}
	}
	return false;
}

int dfs(int x,int mf){//多路增广优化
	if(x==T) return mf;
	int sum=0;
	for(int i=cur[x];i;i=edge[i].nxt){//当前弧优化必不可少
		int y=edge[i].to;cur[x]=i;
		if(dep[y]==dep[x]+1&&edge[i].c){
			int f=dfs(y,min(mf,edge[i].c));
			edge[i].c-=f,edge[i^1].c+=f;
			sum+=f,mf-=f;
			if(mf==0) break;
		}
	}
	if(sum==0) dep[x]=0;
	return sum;
}

int dinic(){
	int flow=0;
	while(bfs()){
		memcpy(cur,head,sizeof(head));
		flow+=dfs(S,1e9);
	}
	return flow;
}

signed main(){
	
	n=read(),m=read(),S=read(),T=read();
	for(int i=1;i<=m;i++){int x=read(),y=read(),c=read();add(x,y,c),add(y,x,0);}
	cout<<dinic()<<endl;
	
	return 0;
}```