【洛谷 P4897】【模板】最小割树（Gomory-Hu Tree）

传送门

题目

题目背景
模板题。做本题之前请确保你会Dinic或ISAP。如果你乱搞过了我请你抽烟

根据相关法律法规，网络流题不允许卡Dinic/ISAP，但可以卡EK，本题数据严格遵循上述条约

题目描述
给定一个nn个点mm条边的无向连通图，多次询问两点之间的最小割

两点间的最小割是这样定义的：原图的每条边有一个割断它的代价，你需要用最小的代价使得这两个点不连通

输入格式
第一行两个数n,mn,m

接下来mm行，每行33个数u,v,wu,v,w，表示有一条连接uu与vv的无向边，割断它的代价为ww

接下来这一行有一个整数QQ，表示询问次数

接下来QQ行，每行两个数u,vu,v，你需要求出uu与vv之间的最小割

注意：因为数据有误，给定图的真实点数应该是n+1n+1个，编号为00到nn

输出格式
输出共QQ行，每行一个整数对应询问的答案

输入输出样例
输入 #1 复制
4 5
1 2 2
2 3 2
4 2 3
4 3 1
1 3 1
3
1 4
2 4
2 3
输出 #1 复制
3
4
4
说明/提示
n\leq 500,\quad m\leq 1500,\quad Q\leq 10^5,\quad 0\leq w\leq 10^4n≤500,m≤1500,Q≤10
5
,0≤w≤10
4

思路

最小割树通常用于求任意两点间的最小割。

首先给出一些结论：

结论一：
mincut(a,b)>=min(mincut(a,c),mincut(c,b))
证明：
对于三点a,b,c，设mincut(a,b),mincut(a,c),mincut(c,b)中最小的为mincut(a,b)。
不妨假设c在残留网络的b端，可以得到mincut(a,c)<=mincut(a,b)
而由假设可得mincut(a,c)>=mincut(a,b)
故mincut(a,c)=mincut(a,b)
事实上这三个值必然是两个相同的较小值和一个较大值。
推论：mincut(a,b)>=min(mincut(a,c1),mincut(c1,c2),…,mincut(ck,b))

通过上述结论可得：对于一棵树，树边(u,v)的边权mincut(u,v)
则对于任意两点a,b，mincut(a,b)=a到b路径上边权的最小值
那么只要建出最小割树，就可以求出任意两点间的最小割了。

咋求？？
这其实是一个分治算法。

步骤一：任意选定一个源节点和一个目的节点。在本例中不失一般性选择节点1为源节点(s)，2为目的节点（t）。则可得最大流，且最小割相应的将点分为两个集合，s集和t集。（s集为残余网络中s能到达的点，t集同理）

步骤二：对于每一个集合，重复步骤一。

步骤三：重复步骤2直到最后每一个点各自在一个集合。

显然这只需跑n-1次最小割

代码

//主要部分
#include<bits/stdc++.h>
using namespace std;
int dfs(int s,int t,int f)
{
	if(s==t||!f) return f;
	int u=s,ret=0;
	for(int i=list[u]; i; i=e[i].next)
	{
		int v=e[i].to;
		if(b[v]!=b[u]+1||!e[i].c) continue;
		int r=dfs(v,t,min(f,e[i].c));
		e[i].c-=r; e[i^1].c+=r; ret+=r; f-=r;
	}
	return ret;
}
bool bfs(int s,int t)
{
	queue<int> q;
	memset(b,0,sizeof(b));
	q.push(s); b[s]=1;
	while(!q.empty())
	{
		int u=q.front(); q.pop();
		for(int i=list[u]; i; i=e[i].next)
		{
			int v=e[i].to;
			if(!e[i].c||b[i]) continue;
			b[v]=b[u]+1;
			q.push(v);
		}
	}
	return b[t];
}
int dinic(int s,int t)
{
	int ret=0;
	while(bfs(s,t))
	{
		ret+=dfs(s,t,0x3f3f3f3f);
	}
	return ret;
}
void build(int *p,int n)
{
	if(n==1) return;
	memset(q1,0,sizeof(q1)); memset(q2,0,sizeof(q2));
	int p1=0,p2=0,val=dinic(p[1],p[2]);
	add2(x,y,val); add2(y,x,val);
	for(int i=1; i<=n; i++)
	{
		if(b[i]) q1[++p1]=p[i];
		else q2[++p2]=p[i];
	}
	for(int i=1; i<=p1; i++) p[i]=q1[i];
	for(int i=1; i<=p2; i++) p[i+p1]=q2[i];
	if(p1) build(p+1,p1);
	if(p2) build(p+p1+1,p2);
}
void dfs2(int u,int f)
{
	for(int i=Last[u]; i; i=E[i].next)
	{
		int v=E[i].to;
		if(v==f) continue;
		f[0][v]=u; dep[v]=dep[u]+1; mn[0][v]=E[i].c;
	}
}
void init()
{
    for (int i = 1; i < K; i++) {
        for (int j = 1; j <= n; j++) {
            f[i][j] = f[i - 1][f[i - 1][j]];
            mn[i][j] = min(mn[i - 1][j], mn[i - 1][f[i - 1][j]]);
        }
    }
}
int solve(int x, int y) {
    int ret = INF;
    if (dep[x] > dep[y]) swap(x, y);
    int sub = dep[y] - dep[x];
    for (int i = K - 1; i >= 0; i--)
        if (sub & (1 << i)) ret = min(ret, mn[i][y]), y = f[i][y];
    if (x == y) return ret;
    for (int i = K - 1; i >= 0; i--) {
        if (f[i][x] != f[i][y]) {
            ret = min(ret, min(mn[i][x], mn[i][y]));
            x = f[i][x], y = f[i][y];
        }
    }
    return min(ret, min(mn[0][x], mn[0][y]));
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1,x,y,c; i<=m; i++) scanf("%d%d",&x,&y,&c),add(x,y,c),add(y,x,c);
	for(int i=1; i<=n; i++) id[i]=i;
	build(id+1,n);
	dep[1]=1;
	dfs2(1,0);
	init();
	scanf("%d",&q);
	while(q--)
	{
		scanf("%d%d",&x,&y);
		printf("%d\n",solve(x,y));
	}
}