[ZJOI2011]最小割【最小割树】

题目描述：

带权无向图，q次询问，每次询问有多少点对的最小割不超过x
n<=150,m<=3000

题目分析：

建出最小割树后每个点dfs一下求出任意两点的最小割，排序之后upper_bound即可。

我写的是论文里的Gusfield(等价流树)。

简单介绍一下最小割树吧。
图截自 2016国家集训队论文《浅谈无向图最小割问题的一些算法及应用》：



我在另一篇博客里面写了证明：



这个证明比较复杂，具体请参阅论文。

Gomory-Hu树的构造也不是很好懂，略过。

递归版：

#include<cstdio>
#include<cctype>
#include<cstring>
#include<queue>
#include<cmath>
#include<algorithm>
#define LL long long
#define maxn 155
#define maxm 6005
using namespace std;
char cb[1<<15],*cs,*ct;
#define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<15,stdin),cs==ct)?0:*cs++)
inline void read(int &a){
	char c;while(!isdigit(c=getchar()));
	for(a=c-'0';isdigit(c=getchar());a=a*10+c-'0');
}
const int inf = 0x3f3f3f3f;
int n,m,S,T;
int fir[maxn],cur[maxn],dis[maxn],nxt[maxm],to[maxm],tot=1,cap[maxm],sum;
inline void line(int x,int y,int z,int rz=0){
	nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y,cap[tot]=z;
	nxt[++tot]=fir[y],fir[y]=tot,to[tot]=x,cap[tot]=rz;
}
queue<int>q;
bool bfs()
{
	memset(dis,0,(n+1)<<2);
	dis[T]=1,q.push(T);
	while(!q.empty()){
		int u=q.front();q.pop();
		for(int i=fir[u];i;i=nxt[i]) if(cap[i^1]&&!dis[to[i]]){
			dis[to[i]]=dis[u]+1;
			q.push(to[i]);
		}
	}
	return dis[S];
}
int dfs(int u,int lim)
{
	if(u==T) return lim;
	int need=lim,delta;
	for(int &i=cur[u];i;i=nxt[i]) if(cap[i]&&dis[u]==dis[to[i]]+1){
		delta=dfs(to[i],min(cap[i],need));
		cap[i]-=delta,cap[i^1]+=delta;
		if(!(need-=delta)) break;
	}
	return lim-need;
}
int Dinic(){
	int flow=0;
	while(bfs()) memcpy(cur,fir,(n+1)<<2),flow+=dfs(S,inf);
	return flow;
}
int Test,Q,x,y,z,s,ans[maxn][maxn],id[maxn],tmp[maxn];
void solve(int l,int r){
	if(l>=r) return;
	for(int i=2;i<=tot;i+=2) cap[i]=cap[i+1]=(cap[i]+cap[i+1])>>1;
	S=id[l],T=id[r];
	int ret=Dinic();
	for(int i=1;i<=n;i++) if(dis[i])
		for(int j=1;j<=n;j++) if(!dis[j])
			ans[i][j]=ans[j][i]=min(ans[i][j],ret);
	int h=l-1,t=r;
	for(int i=l;i<=r;i++) if(dis[id[i]]) tmp[++h]=id[i]; else tmp[t--]=id[i];
	for(int i=l;i<=r;i++) id[i]=tmp[i];
	solve(l,h),solve(h+1,r);
}
int main()
{
	read(Test);
	while(Test--){
		memset(fir,0,sizeof fir),tot=1;
		memset(ans,0x3f,sizeof ans);
		read(n),read(m);
		for(int i=1;i<=m;i++) read(x),read(y),read(z),line(x,y,z,z);
		for(int i=1;i<=n;i++) id[i]=i;
		solve(1,n);
		read(Q);
		while(Q--){
			read(x),s=0;
			for(int i=1;i<=n;i++) for(int j=i+1;j<=n;j++) if(ans[i][j]<=x) s++;
			printf("%d\n",s);
		}
		puts("");
	}
}

非递归版：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define maxn 155
#define maxm 6005
using namespace std;
const int inf = 0x3f3f3f3f;
int n,m,S,T,fa[maxn],dis[maxn][maxn];
bool vis[maxn];
int fir[maxn],cur[maxn],nxt[maxm],to[maxm],c[maxm],tot=1;
vector<int>G[maxn];
inline void line(int x,int y,int z,int rz=0){
	nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y,c[tot]=z;
	nxt[++tot]=fir[y],fir[y]=tot,to[tot]=x,c[tot]=rz;
}
namespace Maxflow{
	int d[maxn],vd[maxn],sz;
	int aug(int u,int augco){
		if(u==T) return augco;
		int need=augco,delta;
		for(int &i=cur[u];i;i=nxt[i]) if(c[i]&&d[u]==d[to[i]]+1){
			delta=aug(to[i],min(need,c[i]));
			c[i]-=delta,c[i^1]+=delta;
			if(!(need-=delta)||d[S]==sz) return augco-need;
		}
		cur[u]=fir[u];
		if(!--vd[d[u]]) d[u]=sz;
		else vd[++d[u]]++;
		return augco-need;
	}
	int SAP(){
		sz=n,memset(d,0,(sz+1)<<2),memset(vd,0,(sz+1)<<2),memset(cur,0,(sz+1)<<2);
		int flow=0;
		while(d[S]<sz) flow+=aug(S,inf);
		return flow;
	}
}
inline void recover(){
	for(int i=2;i<=tot;i+=2) c[i]=c[i^1]=(c[i]+c[i^1])>>1;
}
void dfs(int u){
	vis[u]=1;
	for(int i=fir[u];i;i=nxt[i]) if(c[i]&&!vis[to[i]]) dfs(to[i]);
}
void floyd(int u,int st,int ff,int md){
	dis[st][u]=min(dis[st][u],md);
	for(int i=G[u].size()-1,v;i>=0;i--) if((v=G[u][i])!=ff) floyd(v,st,u,min(md,dis[u][v]));
}
void Gusfield(){
	for(int i=2;i<=n;i++) fa[i]=1;
	for(int u=2,v;u<=n;u++){
		v=fa[u],S=u,T=v,recover();
		dis[u][v]=dis[v][u]=Maxflow::SAP(),G[v].push_back(u),G[u].push_back(v);
		memset(vis,0,(n+1)),dfs(u);
		for(int i=u+1;i<=n;i++) if(fa[i]==v&&vis[i]) fa[i]=u;
	}
	for(int i=1;i<=n;i++) floyd(i,i,0,inf);
}
int tmp[maxn*maxn],cnt;
int main()
{
	int Test,x,y,z,Q;
	scanf("%d",&Test);
	while(Test--){
		scanf("%d%d",&n,&m);
		memset(fir,0,(n+1)<<2),tot=1;
		memset(dis,0x3f,sizeof dis);
		for(int i=1;i<=n;i++) G[i].clear();
		while(m--) scanf("%d%d%d",&x,&y,&z),line(x,y,z,z);
		Gusfield();
		cnt=0;
		for(int i=1;i<=n;i++) for(int j=i+1;j<=n;j++) tmp[++cnt]=dis[i][j];
		sort(tmp+1,tmp+1+cnt);
		scanf("%d",&Q);
		while(Q--) scanf("%d",&x),printf("%d\n",upper_bound(tmp+1,tmp+1+cnt,x)-tmp-1);
		putchar('\n');
	}
}