P3376 【模板】网络最大流

题目

题目

思路

真是神奇，一个初一生智推最大流（初一巨爷请忽略
其实是这么考虑，最大流的解法是由暴力推出来，暴力肯定会TLE，所以需要优化。

---

本文使用Dinic

---

这时我们要理解最大流算法的精髓：撤销，即建反边

大多数人对于建反边这个操作不甚理解，其实它就是所谓的撤销思路，请看下图：

这是原图，其中蓝点为s，黄点为t。
根据输入顺序不同，有可能第一次搜到的是这样的：

显然这不是最优解，dfs可以求出最优解，但是需要通过回溯取消后面的边来求，费时。
那么怎样才不费时呢，回溯是一种取消，但是我们可以通过建反边来撤销，这样就不再需要回溯了。

---

我们先看看建反边后怎么样：

一开始，所有的反边权值均为0。按第一次跑，跑完后路上的反边全部加上了这次新增流量，而原来路上的边全部减去了这次新增流量。而接下来继续dfs的时候，不为0的反边也可以跑了。
接下来我们跑出一个Z型（使用了中间的反边）然后正反2边一抵消，相当于中间的边没有用过（也就是撤销了中间边的使用），但同时我们又增加了最大流量，而且抵消后的2条路依然是从s到t！

---

但同时会有一个问题：如果跑dfs时，一直在同一条边的正反摩擦，怎么办？
没关系，我们可以在每一趟dfs前用一个bfs先求层（注意这里求层只能使用不为0的边，不论正反），只能从上一层去下一层，如果某个时刻已经到不了t了就可以结束了。

---

我们还可以进一步优化：其实每一趟dfs就是在跑一个残量网络，那么在当前情况下，如果某个点去不了t，那么显然这个点就废了，没必要在下面继续对该点dfs。

---

code:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
using namespace std;
int n,m,s,t,x,y,tot=2,head[201],dep[201],l,r,u[201],r2;
long long w,ans;
long long mn(long long x,long long y)
{
	return (x>y?y:x);
}
struct f{
	int to,net;
	long long w;
} a[10011];
void add(int x,int y,long long w)
{
	a[tot].to=y,a[tot].w=w,a[tot].net=head[x],head[x]=tot++;
	return;
}
bool bfs()
{
	memset(dep,0,sizeof(dep));
	l=0,r=0;
	u[r++]=s;
	dep[s]=1;
	while (l<r)
	{
		r2=r;
		for (int i=l;i<r;i++)
		{
			for (int j=head[u[i]];j;j=a[j].net)
			{
				if (a[j].w!=0&&dep[a[j].to]==0)
				{
					u[r2++]=a[j].to;
					dep[a[j].to]=dep[u[i]]+1;
				}
			}
		}
		l=r,r=r2;
	}
	return dep[t];
}
long long dfs(int d,long long in)
{
	if (d==t) return in;
	long long out=0;
	int uw=dep[d]+1;
	for (int j=head[d];j&&in;j=a[j].net)
	{
		if (a[j].w==0||dep[a[j].to]!=uw) continue;
		long long s=dfs(a[j].to,mn(in,a[j].w));
		out+=s,in-=s;
		a[j].w-=s,a[j^1].w+=s;
	}
	if (out==0)
	{
		dep[d]=0;
	}
	return out;
}
int main()
{
	scanf("%d%d%d%d",&n,&m,&s,&t);
	for (int i=0;i<m;i++)
	{
		scanf("%d%d%lld",&x,&y,&w);
		add(x,y,w);
		add(y,x,0);
	}
	while (bfs())
	{
		ans+=dfs(s,1e18);//s是源点，进入它的总量应该是无限大
		//这里用1e18来代替，因为1e18已经大于所有可能值了
	}
	cout<<ans;
    return 0;
}