[学习笔记]最小割树(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;
}
[学习笔记]最小割树(Gomory-Hu Tree)的更多相关文章
- 最小割树Gomory–Hu tree
fanhq666地址:http://fanhq666.blog.163.com/blog/static/8194342620113495335724/ wiki地址(证明):https://en.wi ...
- bzoj 4519: [Cqoi2016]不同的最小割【最小割树Gomory–Hu tree】
算法详见:http://www.cnblogs.com/lokiii/p/8191573.html 求出点两两之间的最小割之后,把他们扔到map/set里跑即可 可怕的是map和set跑的时间竟然完全 ...
- bzoj 2229: [Zjoi2011]最小割【Gomory–Hu tree最小割树】
这个算法详见http://www.cnblogs.com/lokiii/p/8191573.html 求出两两之间最小割之后暴力统计即可 #include<iostream> #inclu ...
- 【模板】最小割树(Gomory-Hu Tree)
传送门 Description 给定一个\(n\)个点\(m\)条边的无向连通图,多次询问两点之间的最小割 两点间的最小割是这样定义的:原图的每条边有一个割断它的代价,你需要用最小的代价使得这两个点不 ...
- [学习笔记]最小割之最小点权覆盖&&最大点权独立集
最小点权覆盖 给出一个二分图,每个点有一个非负点权 要求选出一些点构成一个覆盖,问点权最小是多少 建模: S到左部点,容量为点权 右部点到T,容量为点权 左部点到右部点的边,容量inf 求最小割即可. ...
- 最小割树(Gomory-Hu Tree)
当我们遇到这样的问题: 给定一个 \(n\) 个点 \(m\) 条边的无向连通图,多次询问两点之间的最小割 我们通常要用到最小割树. 博客 建树 分治.记录当前点集,然后随便找俩点当 \(s\) 和 ...
- LoibreOJ 2042. 「CQOI2016」不同的最小割 最小割树 Gomory-Hu tree
2042. 「CQOI2016」不同的最小割 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据 题目描述 ...
- 最小割树(Gomory-Hu Tree)求无向图最小割详解 附 BZOJ2229,BZOJ4519题解
最小割树(Gomory-Hu Tree) 前置知识 Gomory-Hu Tree是用来解决无向图最小割的问题的,所以我们需要了解无向图最小割的定义 和有向图类似,无向图上两点(x,y)的割定义为一个边 ...
- [模板]最小割树(Gomory-Hu Tree)(luogu4897)
给定一个\(n\)个点\(m\)条边的无向连通图,多次询问两点之间的最小割 两点间的最小割是这样定义的:原图的每条边有一个割断它的代价,你需要用最小的代价使得这两个点不连通 Input 第一行两个数\ ...
随机推荐
- oracle中hint 详解
Hint概述 基于代价的优化器是很聪明的,在绝大多数情况下它会选择正确的优化器,减轻了DBA的负担.但有时它也聪明反被聪明误,选择了很差的执行计划,使某个语句的执行变得奇慢无比. 此时就需要DBA进行 ...
- idea的spring配置
开始使用idea,开始的时候把相关的插件都禁用了,导致在建项目和configuration的时候不出现spring相关字样 到plugins中找到installed的插件,查看spring boot ...
- 2019-8-31-How-to-use-code-to-exit-the-application-in-UWP
title author date CreateTime categories How to use code to exit the application in UWP lindexi 2019- ...
- OK。第一个shell script 脚本
很简单吧.也很容易理解.读取两个字符放到变量firstname和lastname中,然后输出 编辑日期格式 上面是编辑日期格式.下面的例题是按照日期格式来分类创建文件
- Nginx---配置详解(转发)
转载于 http://baijiahao.baidu.com/s?id=1604485941272024493&wfr=spider&for=pc 1.概述 Nginx配置文件的整体 ...
- pandas for python
http://pandas.pydata.org/pandas-docs/stable/user_guide/index.html 不算太难,需要拿一本线性代数看看矩阵原理即可.重点在于考虑如何运用, ...
- Ubuntu下终端命令安装sublime
Ubuntu下终端命令安装sublime出现软件包无法定位 sublime-text-install 且多次换源不成功 建议采用离线安装 安装教程如下 用Ubuntu上的浏览器下载一个 Sublime ...
- 「NOI2018」屠龙勇士 解题报告
「NOI2018」屠龙勇士 首先对于每个龙用哪个剑砍,我们可以用set随便模拟一下得到. 然后求出拿这个剑砍这条龙的答案 \[ atk_ix-p_iy=a_i \] 其中\(atk_i\)是砍第\(i ...
- springboot设置静态资源缓存一年
由于本人所在项目组,所用的项目是前后端分离的,前端是React 的SPA,每次打包都会新版本的静态文件. 然而,在有些时候,这些静态资源是不变的,故可以将资源缓存至用户本地,提升性能. 写法如下,需要 ...
- centos7下命令行配置nginx
本教程中的步骤要求用户拥有root权限 第一步 - 添加Nginx存储库要添加CentOS 7 EPEL仓库,请打开终端并使用以下命令: sudo yum install epel-release第二 ...