最小割树（Gomory-Hu Tree），任意两点最小割

正题

Portal
可以得到这个定理：用一组最小割将S集和T集分开，则S集的一个点到T集的一个点的最小割权值都等于这个最小割权值。
具体证明大概是考虑三个点之间的割：设三个点$a,b,c$
若割$(a,b)$后，$c$在a集，那么存在割$(a,b)\ge$割$(c,b)$
很好证明，因为$c$最多走回$a$之后然后再开始割，这样得出来的割一定不会比原来大，如果$c$在$b$集是相同的，因为这个定理对于任三点都适用，所以可以得到割$(a,b)=$割$(c,b)$。
那么我们每次选两个点出来跑最小割，剩下的往两边继续分治就可以了，注意，分治在做的时候也是按全图来跑。
建出一棵最小割树在上面做LCA查询就好了。
时间复杂度$O(n^{3}m)$

#include<bits/stdc++.h>
using namespace std;

const int N=510,M=1510;
struct edge{
	int y,nex,c;
}s[M<<2];
vector<pair<int,int> > nex[N];
int first[N],len=1,n,m,a[N],S,T,tot,b[N],head[N],fa[N][9],mmin[N][9],q;
int qs[N],st,ed,d[N],dep[N];

void ins(int x,int y,int c){
	s[++len]=(edge){y,first[x],c};first[x]=len;
	s[++len]=(edge){x,first[y],0};first[y]=len;
}

bool bfs(){
	memset(d,0,sizeof(d));
	qs[st=ed=1]=S;d[S]=1;
	for(int i=1;i<=n;i++) head[i]=first[i];
	while(st<=ed){
		int x=qs[st++];
		for(int i=first[x];i!=0;i=s[i].nex) if(s[i].c && d[s[i].y]==0)
			d[s[i].y]=d[x]+1,qs[++ed]=s[i].y;
	}
	return d[T]!=0;
}

int dfs(int x,int t){
	if(x==T) return t;
	int tot=0;
	for(int&i=head[x];i!=0;i=s[i].nex) if(s[i].c && d[s[i].y]==d[x]+1){
		int my=dfs(s[i].y,min(t-tot,s[i].c));
		tot+=my;s[i].c-=my;s[i^1].c+=my;
		if(t==tot) break;
	}
	return tot;
}

int Dinic(){
	int tot=0,dx=0;
	while(bfs()){
		dx=dfs(S,1e9);
		while(dx) tot+=dx,dx=dfs(S,1e9);
	}
	return tot;
}

void bt(int l,int r){
	if(l==r) return ;
	S=a[l];T=a[l+1];
	int tmp=Dinic();
	nex[a[l]].push_back(make_pair(a[l+1],tmp));
	nex[a[l+1]].push_back(make_pair(a[l],tmp));
	b[0]=0;for(int i=l;i<=r;i++) if(d[a[i]]) b[++b[0]]=a[i];
	int mid=b[0]+l-1;
	
	for(int i=l;i<=r;i++) if(!d[a[i]]) b[++b[0]]=a[i];
	for(int i=l;i<=r;i++) a[i]=b[i-l+1];
	for(int i=2;i<=len;i+=2) s[i].c+=s[i^1].c,s[i^1].c=0;
	bt(l,mid);bt(mid+1,r);
}

void gt(int x){
	for(int i=0;i<nex[x].size();i++) if(nex[x][i].first!=fa[x][0]){
		int y=nex[x][i].first,c=nex[x][i].second;
		fa[y][0]=x;mmin[y][0]=c;dep[y]=dep[x]+1;
		for(int k=1;k<=8;k++) fa[y][k]=fa[fa[y][k-1]][k-1],mmin[y][k]=min(mmin[y][k-1],mmin[fa[y][k-1]][k-1]);
		gt(y);
	}
}

int solve(int x,int y){
	if(dep[x]>dep[y]) swap(x,y);
	int ans=1e9;
	for(int k=8;k>=0;k--) if(dep[fa[y][k]]>=dep[x]) ans=min(ans,mmin[y][k]),y=fa[y][k];
	if(x==y) return ans;
	for(int k=8;k>=0;k--) if(fa[x][k]!=fa[y][k]) ans=min(ans,min(mmin[x][k],mmin[y][k])),x=fa[x][k],y=fa[y][k];
	return min(ans,min(mmin[x][0],mmin[y][0]));
}

int main(){
	scanf("%d %d",&n,&m);
	int x,y,c;
	for(int i=1;i<=m;i++) scanf("%d %d %d",&x,&y,&c),ins(x,y,c),ins(y,x,c);
	for(int i=1;i<=n;i++) a[i]=i;
	bt(1,n);dep[1]=1;gt(1);
	scanf("%d",&q);
	while(q--){
		scanf("%d %d",&x,&y);
		printf("%d\n",solve(x,y));
	}
}