初学最小割树

前言

证明什么的一概不会，问就是感性理解，总之先坑了。。。

最小割树的构造

任取当前点集中的两个点\(s,t\)，求出\(s\)与\(t\)在原图上的任意一组最小割（设大小为\(v\)），然后就可以在最小割树上给\(s,t\)两点间连一条边权为\(v\)的边。

接着，把点集按照最小割分成两部分（即与\(s\)连通的一部分和与\(t\)连通的一部分），对两部分的点分别递归构造。

显然这样一定能构造出一棵树，这棵树就叫最小割树。

最小割树的性质

任意两点\(x,y\)间的最小割就是最小割树上\(x,y\)间路径中的最小边权。

因此只要预处理一下就可以倍增查询了。

代码：\(O(n^3m)\)

#include<bits/stdc++.h>
#define Tp template<typename Ty>
#define Ts template<typename Ty,typename... Ar>
#define Reg register
#define RI Reg int
#define Con const
#define CI Con int&
#define I inline
#define W while
#define N 500
#define M 1500
#define LN 9
#define INF (int)1e9
using namespace std;
int n,m;
namespace FastIO
{
	#define FS 100000
	#define tc() (FA==FB&&(FB=(FA=FI)+fread(FI,1,FS,stdin),FA==FB)?EOF:*FA++)
	#define pc(c) (FC==FE&&(clear(),0),*FC++=c)
	int OT;char oc,FI[FS],FO[FS],OS[FS],*FA=FI,*FB=FI,*FC=FO,*FE=FO+FS;
	I void clear() {fwrite(FO,1,FC-FO,stdout),FC=FO;}
	Tp I void read(Ty& x) {x=0;W(!isdigit(oc=tc()));W(x=(x<<3)+(x<<1)+(oc&15),isdigit(oc=tc()));}
	Ts I void read(Ty& x,Ar&... y) {read(x),read(y...);}
	Tp I void writeln(Ty x) {W(OS[++OT]=x%10+48,x/=10);W(OT) pc(OS[OT--]);pc('\n');}
}using namespace FastIO;
namespace D//网络流
{
	#define add(x,y,f) (_[++ee].nxt=lnk[x],_[lnk[x]=ee].to=y,_[ee].F=f)
	int s,t,ee=1,lnk[N+5],cur[N+5];struct edge {int to,nxt,F;}e[2*M+5],_[2*M+5];
	I void Add(CI x,CI y,CI f) {add(x,y,f),add(y,x,f);}//加一条无向边
	I void Init() {for(RI i=1;i<=ee;++i) e[i]=_[i];}//把原图的边复制一份
	int q[N+5],d[N+5];I bool BFS() {RI i,k,H,T;for(i=1;i<=n;++i) d[i]=0;d[q[H=T=1]=s]=1;
		W(H<=T&&!d[t]) for(i=lnk[k=q[H++]];i;i=e[i].nxt) !d[e[i].to]&&e[i].F&&(d[q[++T]=e[i].to]=d[k]+1);return d[t];}
	I int DFS(CI x=s,RI f=INF) {if(x==t||!f) return f;RI i,o,g=0;for(i=cur[x];i&&f;i=e[i].nxt)
		(d[x]+1)==d[e[i].to]&&(o=DFS(e[i].to,min(f,e[i].F)),e[i].F-=o,e[i^1].F+=o,g+=o,f-=o);return cur[x]=i,g;}
	I int MaxFlow() {RI g=0;W(BFS()) memcpy(cur,lnk,sizeof(lnk)),g+=DFS();return g;}//最小割=最大流
}
class MinCutTree
{
	private:
		#define add(x,y,z) (e[++ee].nxt=lnk[x],e[lnk[x]=ee].to=y,e[ee].v=z)
		int ee,lnk[N+5];struct edge {int to,nxt,v;}e[2*N+5];
		int id[N+5],idl[N+5],idr[N+5];I void Solve(CI l,CI r)//建最小割树
		{
			RI i,v,tl=0,tr=0;if(l>=r) return;D::Init(),D::s=id[l],D::t=id[r],v=D::MaxFlow();//任取两点求一组最小割
			for(add(id[l],id[r],v),add(id[r],id[l],v),i=l;i<=r;++i) (D::d[id[i]]?idl[++tl]:idr[++tr])=id[i];//最小割树上连边；根据最小割划分点集
			for(i=1;i<=tl;++i) id[l+i-1]=idl[i];for(i=1;i<=tr;++i) id[l+tl+i-1]=idr[i];Solve(l,l+tl-1),Solve(l+tl,r);//递归构造
		}
		int d[N+5],f[N+5][LN+1],g[N+5][LN+1];I void dfs(CI x)//预处理
		{
			RI i;for(i=1;i<=LN;++i) f[x][i]=f[f[x][i-1]][i-1],g[x][i]=min(g[x][i-1],g[f[x][i-1]][i-1]);//倍增预处理
			for(i=lnk[x];i;i=e[i].nxt) e[i].to^f[x][0]&&(d[e[i].to]=d[f[e[i].to][0]=x]+1,g[e[i].to][0]=e[i].v,dfs(e[i].to),0);
		}
	public:
		I void Build() {for(RI i=1;i<=n;++i) id[i]=i;Solve(1,n),dfs(1);}
		I int Q(RI x,RI y)//树上倍增
		{
			#define U(x,i) (t=min(t,g[x][i]),x=f[x][i])
			RI i,t=INF;for(d[x]<d[y]&&(swap(x,y),0),i=0;d[x]^d[y];++i) (d[x]^d[y])>>i&1&&U(x,i);
			if(x==y) return t;for(i=LN;~i;--i) f[x][i]^f[y][i]&&(U(x,i),U(y,i));return U(x,0),U(y,0),t;
		}
}T;
int main()
{
	RI i,x,y,z;for(read(n,m),i=1;i<=m;++i) read(x,y,z),D::Add(x,y,z);
	T.Build();RI Qt;read(Qt);W(Qt--) read(x,y),writeln(T.Q(x,y));return clear(),0;
}

标签：FS,割树,int,最小,初学,FO,id,define
来源： https://www.cnblogs.com/chenxiaoran666/p/MinCutTree.html