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

Problem

Solution

• 做模板使用，不证明（不会证 。
• 原图中两点间的最小割为最小割树上两点间最小割，可通过树上倍增询问。
• 最小割树的建立：任意两点间的最小割将当前点集分成两部分，以此为基础分治，每次从当前点集中任取两点求最小割，将最小割作为Gomory-Hu Tree上两点间的边权，并根据最小割将当前点集分成两部分，s那边的与t那边的，边界就是当前点集中只有一个点。
• 相当于跑了n-1次dinic。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn = 505,maxm = 1505,inf = 1e9;
int cnt = 1,last[maxn],ceng[maxn],cnt0=1,belong[maxn];
int cur[maxn],n,m,cnt1,f[maxn][20],dep[maxn],front[maxn],Q,gm[maxn][20];
struct edge{
	int v,next,w;
}e[maxm<<2];
edge te[maxn<<1];
inline void add(int u,int v,int w)
{
	e[++cnt] = {v,last[u],w};
	last[u] = cnt;
	e[++cnt] = {u,last[v],0};
	last[v] = cnt;
}
inline void add_tree(int u,int v,int w)
{
	te[++cnt1] = {v,front[u],w};
	front[u] = cnt1;
	te[++cnt1] = {u,front[v],w};
	front[v] = cnt1;
}
bool bfs(int s,int t)
{
	queue<int>q;
	memset(ceng,0,sizeof(ceng));
	ceng[s]=1;
	q.push(s);
	while(!q.empty())
	{
		int u=q.front();q.pop();
		for(int i=last[u];i;i=e[i].next)
		{
			int v=e[i].v,w=e[i].w;
			if(ceng[v]==0&&w>0)
			{
				ceng[v]=ceng[u]+1;
				q.push(v);
			}
		}
	}
	return ceng[t]; 
}
int dfs(int u,int dis,int t)
{
	if(u==t)
	return dis;
	int minf = dis;
	for(int i=cur[u];i;i=e[i].next)
	{
		cur[u]=e[i].next;
		int v=e[i].v,w=e[i].w;
		if(ceng[v]==ceng[u]+1&&w!=0)
		{
			int di=dfs(v,min(dis,w),t);
			if(di>0)
			{
				e[i].w-=di;
				e[i^1].w+=di;
				dis-=di;
			}
		}
	}
	return minf-dis;
}
int dinic(int s,int t)
{
	int ans=0;
	while(bfs(s,t))
	{
		for(int i=1;i<=n;i++)
		cur[i]=last[i];
		ans+=dfs(s,inf,t);
	}
	return ans;
}
void init()
{
	for(int i=2;i<=cnt;i+=2)
	{
		e[i].w+=e[i^1].w;
		e[i^1].w = 0; 
	}
}
void dfs0(int u,int fa,int from)
{
	if(belong[u]==from)
	belong[u] = cnt0;
	for(int i=last[u];i;i=e[i].next)
	{
		int v=e[i].v,w=e[i].w;
		if(v!=fa&&w>0&&e[i^1].w>0)
		dfs0(v,u,from);
	}
}
int num[maxn],tmp1[maxn],tmp2[maxn];
void build(int l,int r)
{
	if(l==r) return ;
	
	init();
	int s = num[l],t = num[l+1];
	int flow = dinic(s,t),num1=0,num2=0;
	add_tree(s,t,flow);
	for(int i=l;i<=r;i++)
	{
		if(!ceng[num[i]]) tmp1[++num1]=num[i];
		else tmp2[++num2]=num[i];
	}
	for(int i=l;i<=l+num1-1;i++) num[i] = tmp1[i-l+1];
	for(int i=l+num1;i<=r;i++) num[i] = tmp2[i-l-num1+1];
	build(l,l+num1-1);
	build(l+num1,r);
}
void dfs_tree(int u,int fa)
{
	for(int i=1;i<=15;i++)
	{
		f[u][i] = f[f[u][i-1]][i-1];
		gm[u][i] = min(gm[u][i-1],gm[f[u][i-1]][i-1]);
	}
	for(int i=front[u];i;i=te[i].next)
	{
		int v=te[i].v,w=te[i].w;
		if(v!=fa)
		{
			dep[v] = dep[u]+1;
			f[v][0] = u;
			gm[v][0] = w;
			dfs_tree(v,u);
		}
	}
}
int query(int u,int v)
{
	int ans = inf;
	if(dep[u]<dep[v]) swap(u,v);
	for(int i=15;i>=0;i--)
	{
		if(dep[f[u][i]]>=dep[v])
		{
			ans = min(ans,gm[u][i]);
			u = f[u][i];
		}
	}
	for(int i=15;i>=0;i--)
	{
		if(f[u][i]!=f[v][i])
		{
			ans=min(ans,min(gm[u][i],gm[v][i]));
			u = f[u][i];v = f[v][i];
		}
	}
	if(u!=v)
	ans=min(ans,min(gm[u][0],gm[v][0]));
	return ans;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1,u,v,w;i<=m;i++)
	{
		scanf("%d%d%d",&u,&v,&w);
		add(u,v,w);
		add(v,u,w);
	}
	for(int i=1;i<=n;i++)
	num[i] = i;
	build(1,n);
	dep[1] = 1;
	dfs_tree(1,0);
	scanf("%d",&Q);
	for(int i=1,u,v;i<=Q;i++)
	{
		scanf("%d%d",&u,&v);
		cout<<query(u,v)<<endl;
	}
	return 0;
}