网络流 dinic

        现在只学了最大流，现在对最大流的认识是可以做一些类似有一个源点和一个汇点，源点和汇点之间会有多条道路，问最后汇入汇点的流量最大能是多少。

        想到这个就会去想直接从源点开始搜索不就行了？搜到汇点后就累加一次答案，但是这样可能不是最优的，比如下面这个图

按上面的来的话最大是2，但是答案是4，也很好看出来；所以为了避免这样的问题，就要引入这个算法最核心的地方----后悔机制，思想应该和反悔贪心差不多(反悔贪心早忘光了，，),实现这个机制我们需要引入反向边，给它赋一层较真实的意义的话就是和退流管道一样，就和下面这样一样，一开始我们先让反向边的流量都为0，之后走了一条和上图一样的路(这样流进汇点的路叫做一条增广路)，我们就需要给路过的边的反向边赋予流量了，赋的流量就是该条边的正向边流过的流量，就比如上面的边，正向流过的是2，本来容量就是2，所以正向边还剩0的容量可用，反向边在加上正向流过的流量2，这意思就是这条边有2个流量可以反悔；

 上面这张图就是建了反边之后的，发现它还多了一条路可以走，答案是4，中间的那条边正向走了2，反向又走了2，其实相当于没走这条边；

        明白了这个反悔机制后算法也就比较好理解了，EK算法就是不断用bfs去找增广路然后不断地更新最大值，时间复杂度是O(n*m*m)，鉴于还没用到过就先略写一下；

        dicnic算法算是EK的优化，EK是每一次都要去遍历一遍剩下的边，只要边还有流量就会去遍历，(这些边和节点构成的网络叫做残留网)但是EK遍历完只会找到一条增广路，而dinic会找到多个，然后用dfs去把这些增广路的答案都计算出来，思路就是用bfs将残留网分层，d[u]表示u与源点的距离，d[]相等的一些点被称作是在一层上，一个点d[u]只能去访问的d[v]=d[u]+1的点，这样有利于dfs同时处理多条增广路，感觉这样就像是一棵树一样，去深搜自己的儿子，代码中还有一些优化，直接看代码就可以了

#include <bits/stdc++.h>
#define endl '\n'
using namespace std;
#define int long long
const int N = 1e6+100;
int head[200005],cnt=1;
struct Edge
{
	int next,to,w;
}e[400005];
void addedge(int from,int to,int w)
{
	e[++cnt].to=to;
	e[cnt].w=w;
	e[cnt].next=head[from];
	head[from]=cnt;
}
int n,m,s,t;
int d[200005],cur[200005];
//d为层数，cur[u]表示u这个点第一个该访问的边是哪一个
bool bfs()
{//找增广路
	memset(d,0,sizeof(d));
	queue<int>q;
	q.push(s);d[s]=1;
	while(!q.empty())
	{
		int u=q.front();q.pop();
		for(int i=head[u];i;i=e[i].next)
		{
			int j=e[i].to;
			if(d[j]==0&&e[i].w)
			{
				d[j]=d[u]+1;
				q.push(j);
				if(j==t) return 1;
			}
		}
	}
	return 0;
}
int dfs(int u,int mf)
{
	if(u==t) return mf;//到了汇点，返回流量
	int sum=0;
	for(int i=cur[u];i;i=e[i].next)
	{
		cur[u]=i;//更新u第一个该访问的边，之前的边已经都被榨干了到不了汇点了，所以再走下去也没有意义
		int j=e[i].to;
		if(d[j]==d[u]+1&&e[i].w)
		{
			int f=dfs(j,min(e[i].w,mf));//返回回来的值是流进汇点的流量
			e[i].w-=f;//正向边减去流进汇点的流量f
			e[i^1].w+=f;//反向边加上流进汇点的流量f,以便后悔
			sum+=f;//累加流进汇点的流量，这就是答案
			mf-=f;//mf就是上一层边更新过来的最大流量值，这一条边分了f之后自然要减去这一条边的f
			if(mf==0) break;//说明上一层的那条边已经被榨干了
		}
	}
	if(sum==0) d[u]=0;//如果走了一遭没有流进汇点的流量，那么这个点u也没什么用了，踢掉以后不会再访问了
	return sum;
}
int dinic()
{
	int flow=0;
	while(bfs())
	{
		memcpy(cur,head,sizeof(head));
		flow+=dfs(s,1e18);
	}
	return flow;
}
signed main()
{
    //cin.tie(0);cout.tie(0);ios::sync_with_stdio(0);
	cin>>n>>m>>s>>t;
	for(int i=1;i<=m;i++)
	{
		int u,v,w;
		cin>>u>>v>>w;
		addedge(u,v,w);addedge(v,u,0);
	}
	int ans=dinic();
	cout<<ans<<endl;
	system("pause");
    return 0;
}