[学习笔记]最小割树(Gomory-Hu Tree)

最小割树(\(\mathcal{Gomory-Hu Tree}\))简明指南

对于单源最短路径,我们有\(SPFA\)\(Dijkstra\),对于多源最短路径,我们有\(Floyd\);对于两点间的最小割,我们有\(Dinic\)\(ISAP\),那么对于多组最小割的询问呢?

这是板子题:

P4897 【模板】最小割树(Gomory-Hu Tree)

这是教程:

首先有一个定理,就是一个\(n\)个点的图上,两点之间只有\(n\)种本质不同的最小割。因此一定存在一棵树,满足树上两点的最小割等于原图上两点的最小割。我们把这样的树称之为“最小割树”。 --Ebola

下面考虑如何建出这样的一棵树。

首先我们在图中随意选取两个点\(u,v\),跑出他们之间的最小割\(cut(u,v)\)。那么满流的边就是图中\(u,v\)之间的割。通过这条割我们把图划分成\(u\)所在的部分和\(v\)所在的部分。把\(u\)所在的那一部分的点集记作\(U\)\(v\)所在的那一部分的点集记作\(V\)

  • 引理:对于任意\(x \in U,y \in V\),有\(cut(x,y) \leq cut(u,v)\)
  • 证明:假设有\(cut(x,y) > cut(u,v)\),那么\(cut(u,v)\)就不能把\(u,v\)割开,因为\(x,y\)依然相连。

那么我们在一张你没有玩过的船新图上在\(u,v\)之间建一条值为\(cut(u,v)\)的边。接下来,我们再分别对\(U,V\)集合进行上述的相同操作,也就是在每一集合中找两个点求最小割,再把集合割开成两个...这样,当我们把原图的点的大集合全部割成一个个点集之后,刚好跑了\(n\)次网络流,建了\(n-1\)条边,也就建成了一棵\(Gomory\)-\(Hu \ Tree\)

  • 引理:对于树上的两点\(u,v\),他们的最小割是他们之间简单路径的最小路径权值。
  • 根据上一条引理,显然可得。

那么我们只需要在树上进行查询们就可以得到任意两点的最小割。可以用倍增,也可以用树链剖分。

这是板子:

#include<bits/stdc++.h>
using namespace std;
int n,m,node[505],dep[505],fa[505][10],mn[505][10];
int cnt,top[505],to[1005],len[1005],nex[1005];
int read()
{
    int re=0;
    char ch=getchar();
    while(!isdigit(ch)) ch=getchar();
    while(isdigit(ch)) re=(re<<3)+(re<<1)+ch-'0',ch=getchar();
    return re;
}
void add_edge(int x,int y,int z)
{
    to[++cnt]=y,len[cnt]=z,nex[cnt]=top[x],top[x]=cnt;
    to[++cnt]=x,len[cnt]=z,nex[cnt]=top[y],top[y]=cnt;
}
namespace GHT
{
    int s,t;
    int tot,cur[505],dep[505],col[505],col_bucket[505];
    int cnt=1,top[505],to[3005],cap[3005],flow[3005],nex[3005];
    void add_edge(int x,int y,int z)
    {
        to[++cnt]=y,cap[cnt]=z,flow[cnt]=0,nex[cnt]=top[x],top[x]=cnt;
        to[++cnt]=x,cap[cnt]=z,flow[cnt]=0,nex[cnt]=top[y],top[y]=cnt;
    }
    bool BFS()
    {
        memset(cur,0,sizeof cur);
        memset(dep,0,sizeof dep);
        dep[s]=1,cur[s]=top[s];
        queue<int>Q;
        Q.push(s);
        while(!Q.empty())
        {
            int now=Q.front();Q.pop();
            for(int i=top[now];i;i=nex[i])
                if(!dep[to[i]]&&cap[i]>flow[i])
                {
                    dep[to[i]]=dep[now]+1;
                    cur[to[i]]=top[to[i]];
                    Q.push(to[i]);
                }
        }
        return dep[t]!=0;
    }
    int DFS(int now,int rest)
    {
        if(now==t) return rest;
        int re=0;
        for(int &i=cur[now];i;i=nex[i])
            if(dep[to[i]]==dep[now]+1&&cap[i]>flow[i])
            {
                int lzq=DFS(to[i],min(rest,cap[i]-flow[i]));
                if(lzq)
                {
                    rest-=lzq,re+=lzq;
                    flow[i]+=lzq,flow[i^1]-=lzq;
                    if(!rest) break;
                }
            }
        return re;
    }
    int Dinic(int x,int y)
    {
        int re=0;s=x,t=y;
        for(int i=1;i<=cnt;i++) flow[i]=0;
        while(BFS()) re+=DFS(s,0x3f3f3f3f);
        return re;
    }
    void get_color(int now,int color)
    {
        col[now]=color;
        for(int i=top[now];i;i=nex[i])
            if(cap[i]>flow[i]&&col[to[i]]!=color)
                get_color(to[i],color);
    }
    void build(int l,int r)
    {
        if(l==r) return ;
        int x=node[l],y=node[l+1];
        int cut=Dinic(x,y);
        get_color(x,++tot);
        int L=l,R=r;
        for(int i=l;i<=r;i++)
            if(col[node[i]]==tot) col_bucket[L++]=node[i];
            else col_bucket[R--]=node[i];
        for(int i=l;i<=r;i++) node[i]=col_bucket[i];
        ::add_edge(x,y,cut);
        build(l,L-1);
        build(R+1,r);
    }
}
void dfs(int now)
{
    for(int i=1;i<=9;i++)
    {
        fa[now][i]=fa[fa[now][i-1]][i-1];
        mn[now][i]=min(mn[now][i-1],mn[fa[now][i-1]][i-1]);
    }
    for(int i=top[now];i;i=nex[i])
    {
        if(to[i]==fa[now][0]) continue;
        dep[to[i]]=dep[now]+1,fa[to[i]][0]=now,mn[to[i]][0]=len[i];
        dfs(to[i]);
    }
}
int getcut(int x,int y)
{
    int re=INT_MAX;
    if(dep[x]<dep[y]) swap(x,y);
    for(int i=9;i>=0;i--) if(dep[fa[x][i]]>=dep[y]) re=min(re,mn[x][i]),x=fa[x][i];
    if(x==y) return re;
    for(int i=9;i>=0;i--) if(fa[x][i]!=fa[y][i]) re=min(re,min(mn[x][i],mn[y][i])),x=fa[x][i],y=fa[y][i];
    return min(re,min(mn[x][0],mn[y][0]));
}
int main()
{
    n=read(),m=read();
    while(m--)
    {
        int x=read(),y=read(),z=read();
        GHT::add_edge(x,y,z);
    }
    for(int i=1;i<=n;i++) node[i]=i;
    GHT::build(1,n);
    dep[1]=1;
    dfs(1);
    m=read();
    while(m--)
    {
        int x=read(),y=read();
        printf("%d\n",getcut(x,y));
    }
    return 0;
}
posted @ 2018-10-11 13:34 UranusITS 阅读(...) 评论(...) 编辑 收藏