<更新提示>

<第一次更新>

<第二次更新>新增一道例题及讲解


<正文>

Exkruscal

\(Exkruscal\)又称\(Kruscal\)重构树,是一种利用经典算法\(Kruscal\)来实现的构造算法,可以将一张无向图重构为一棵具有\(2n-1\)个节点的树,这棵树具有许多特殊的性质,可以用来解决许多问题。

那么我们来了解一下这个新算法。

Kruscal

先回顾一下\(Kruscal\)算法吧。这是一个经典的图论算法,用于求解无向图的最小生成树。

其算法大致思路为将边用边表的形式储存,按照权值从小到大排序,然后枚举每一条边,尝试连接不在同一联通块中的两个节点,这一过程用并查集来维护。

重构树

那么我们怎么用\(Kruscal\)算法来重构树呢?我们还是先将边按边权从小到大排序,然后按照\(Kruscal\)算法的流程来,每一次加入一条连接两个不同联通块的点,并将这一过程用并查集来维护。

不同的是,对于每一条\(Kruscal\)算法选中的边,我们要新建一个节点,将该点的点权赋值为选中边的边权,并将该点连向选中边连接的两个点。

这样,对于\(Kruscal\)算法选中的\(n-1\)条边,我们就构造出了一个具有\(2n-1\)个节点的树,其中,树的\(n\)个叶节点对应原图的\(n\)个节点它们没有点权,而其他的所有点均有一个点权。

性质

我们已经了解如何用\(Kruscal\)算法重构树了,那么,这棵树具有如下的性质,可以帮助我们解题。

\(1.\)这棵树具有堆性质(大根堆)

\(2.\)原图两点之间路径的最大边权最小值是重构树上两点\(lca\)的权值

\(3.\)树中代表原图中的点的节点全是叶子节点,其余节点都代表了一条边的边权

\(Code:\)

inline void Exkruscal(void)
{
sort(l+1,l+m+1,compare);//l是边表
for(int i=1;i<=2*n;i++)
fa[i]=i;//并查集初始化
for(int i=1;i<=m;i++)
{
int x=find(l[i].x),y=find(l[i].y);
if(x==y)continue;
fa[x]=fa[y]=++n;
val[n]=l[i].val;//点权
insert(n,x);//连边
insert(n,y);
}
}

Network(BZOJ3732)

Description

给你N个点的无向图 (1 <= N <= 15,000),记为:1…N。

图中有M条边 (1 <= M <= 30,000) ,第j条边的长度为: d_j ( 1 < = d_j < = 1,000,000,000).

现在有 K个询问 (1 < = K < = 15,000)。 每个询问的格式是:A B,表示询问从A点走到B点的所有路径中,最长的边最小值是多少?

Input Format

第一行: N, M, K。

第2..M+1行: 三个正整数:X, Y, and D (1 <= X <=N; 1 <= Y <= N). 表示X与Y之间有一条长度为D的边。

第M+2..M+K+1行: 每行两个整数A B,表示询问从A点走到B点的所有路径中,最长的边最小值是多少?

Output Format

对每个询问,输出最长的边最小值是多少。

Sample Input

6 6 8
1 2 5
2 3 4
3 4 3
1 4 8
2 5 7
4 6 2
1 2
1 3
1 4
2 3
2 4
5 1
6 2
6 1

Sample Output

5
5
5
4
4
7
4
5

解析

\(Solution\ 1\)

由于最小生成树一定是瓶颈生成树,所以我们可以将原图求最小生成树直接重新建图,保证两点之间路径的最大权值最小。然后树上倍增\(LCA\)顺便求最大值即可。

\(Solution\ 2\)

\(Exkruscal\)模板题,直接\(Kruscal\)重构树,然后求\(LCA\)点权即可。

\(Code:\)

#include<bits/stdc++.h>
using namespace std;
#define mset(name,val) memset(name,val,sizeof name)
#define filein(str) freopen(str".in","r",stdin)
#define fileout(str) freopen(str".out","w",stdout)
const int N=15020,M=30020,K=20020,MaxlogN=25;
struct Link{int x,y,val;}l[M];
struct edge{int ver,next;}e[2*N];
int n,m,q,t,fa[2*N],val[2*N],depth[2*N],f[2*N][MaxlogN],Last[2*N];
inline void input(void)
{
scanf("%d%d%d",&n,&m,&q);
for(int i=1;i<=m;i++)
scanf("%d%d%d",&l[i].x,&l[i].y,&l[i].val);
}
inline bool compare(Link a,Link b)
{
return a.val<b.val;
}
inline void insert(int x,int y)
{
e[++t].ver=y;e[t].next=Last[x];Last[x]=t;
}
inline int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
inline void Exkruscal(void)
{
sort(l+1,l+m+1,compare);
for(int i=1;i<=2*n;i++)
fa[i]=i;
for(int i=1;i<=m;i++)
{
int x=find(l[i].x),y=find(l[i].y);
if(x==y)continue;
fa[x]=fa[y]=++n;
val[n]=l[i].val;
insert(n,x);
insert(n,y);
}
}
inline void dfs(int x,int dep)
{
depth[x]=dep;
for(int i=Last[x];i;i=e[i].next)
{
int y=e[i].ver;
if(f[x][0]==y)continue;
f[y][0]=x;
dfs(y,dep+1);
}
}
inline void dp(void)
{
f[n][0]=-1;
for(int k=1;(1<<k)<n;k++)
for(int i=1;i<=n;i++)
if(f[i][k-1]<0)f[i][k]=-1;
else f[i][k]=f[f[i][k-1]][k-1];
}
inline int LCA(int x,int y)
{
if(depth[x]>depth[y])
swap(x,y);
for(int d=depth[y]-depth[x],i=0;d;d>>=1,i++)
if(1&d)y=f[y][i];
if(x==y)return x;
for(int i=MaxlogN-1;i>=0;i--)
{
if(f[x][i]^f[y][i])
{
x=f[x][i];
y=f[y][i];
}
}
return f[x][0];
}
inline void solve(void)
{
for(int i=1;i<=q;i++)
{
int a,b;
scanf("%d%d",&a,&b);
printf("%d\n",val[LCA(a,b)]);
}
}
int main(void)
{
input();
Exkruscal();
dfs(n,1);
dp();
solve();
return 0;
}

路径权值

Description

给定一个带权树,树上任意两点间的路径权值d(x,y)定义为x,y这两个点之间路径上的最小值,树上任意一点x的权值定义为这个点到树上其他所有点的路径权值和,即∑d(x,i),现求树上每一个点的路径权值和。

Input Format

首先输入一个整数 n(1≤n≤100000),表示树的点的个数。 

接下来n−1行,每行三个整数 x,y,s(1≤x,y≤n,1≤s≤1000),表示编号为x的节点和编号为y的节点之间存在一条权值为s的边,树上每个点的编号为1~n

Output Format

n行,每行输出每个节点对应的路径权值和

Sample Input

4
1 2 2
2 4 1
2 3 1

Sample Output

4
4
3
3

解析

看到两点路径间的最小值就想到是可以用\(Exkruscal\)的。

先将边按边权从大到小排序,重构树一下,两点\(LCA\)的权值就是树上两点路径的最小值。对于重构树中的每一个非叶子节点\(x\),它的左子树中的每一个叶节点到右子树中的每一个叶节点的真实路径中必然经过\(x\)所代表的边。

我们设\(p_i\)代表原树中节点\(i\)的路径权值和,那么对于刚才提到的一个非叶子节点\(x\),\(x\)的左子树中所有叶子结点的\(p\)值都应该加上\(x\)右子树大小乘\(val_x\),右子树也是同理(加上\(x\)左子树大小乘\(val_x\))。

最后,我们要输出每个节点的\(p\)值。这样的话,问题就转换成了我们需要一个数据结构,能够实现区间加法和单点查询,可以用树状数组配合差分实现。

还有一个问题,如果做区间加法时,子树中叶节点所对应的编号不连续怎么办,我们可以前序遍历给每一个叶子节点重新标一个编号,用于树状数组,最后输出时输出每个节点对应编号的值即可。

\(Code:\)

#include<bits/stdc++.h>
using namespace std;
const int N=100000+20;
int dif[N],n,n_,Last[N*2],t,f[N*2],val[N*2],size[N*2],cnt;
struct LINK{int x,y,val;}l[N];
struct EDGE{int ver,next;}e[N*2];
struct TREE{int l,r;}tree[N*2];
inline bool compare(LINK a,LINK b)
{
return a.val>b.val;
}
inline int find(int x)
{
return f[x]==x?x:f[x]=find(f[x]);
}
inline void insert(int x,int y)
{
e[++t].ver=y;e[t].next=Last[x];Last[x]=t;
}
inline int lowbit(int x){return x&(-x);}
inline void modify(int pos,int val)
{
for(;pos<=n_;pos+=lowbit(pos))dif[pos]+=val;
}
inline int query(int pos)
{
int res=0;
for(;pos;pos-=lowbit(pos))res+=dif[pos];
return res;
}
inline void input(void)
{
scanf("%d",&n);
for(int i=1;i<n;i++)
scanf("%d%d%d",&l[i].x,&l[i].y,&l[i].val);
n_=n;
}
inline void Exkruscal(void)
{
sort(l+1,l+n,compare);
for(int i=1;i<=2*n-1;i++)
f[i]=i;
for(int i=1;i<n;i++)
{
int x=find(l[i].x),y=find(l[i].y);
f[x]=f[y]=++n_;
val[n_]=l[i].val;
insert(n_,x);insert(n_,y);
}
}
inline void dfs(int x)
{
size[x]=(x<=n);
tree[x].l=++cnt;
for(int i=Last[x];i;i=e[i].next)
{
int y=e[i].ver;
dfs(y);
size[x]+=size[y];
}
tree[x].r=cnt;
}
inline void solve(void)
{
for(int i=n+1;i<=n_;i++)
{
int lson=0,rson=0;
for(int j=Last[i];j;j=e[j].next)
{
if(lson)rson=e[j].ver;
else lson=e[j].ver;
}
modify(tree[lson].l,size[rson]*val[i]);
modify(tree[lson].r+1,-size[rson]*val[i]);
modify(tree[rson].l,size[lson]*val[i]);
modify(tree[rson].r+1,-size[lson]*val[i]);
}
}
int main(void)
{
input();
Exkruscal();
dfs(n_);
solve();
for(int i=1;i<=n;i++)
printf("%d\n",query(tree[i].l));
return 0;
}

<后记>

『Kruscal重构树 Exkruscal』的更多相关文章

  1. loj2876 水壶 [JOISC 2014 Day2] kruscal重构树

    正解:kruscal重构树+bfs 解题报告: 我永远喜欢loj! 感觉这题和这题挺像的,,,预处理和解题方法都是,,,所以大概整体二分能过去? 但因为做这题主要是入门一下kruscal重构树,,,所 ...

  2. kruscal重构树略解

    我们先看一道题:Luogu P4197 Peaks 这道题珂以用启发式合并+主席树来做 那么强制在线呢?(bzoj 3551 [ONTAK2010]Peaks加强版) 离线做法就不行了 我们就要用一个 ...

  3. BZOJ_3545_[ONTAK2010]Peaks_主席树+倍增+kruscal重构树+dfs序

    BZOJ_3545_[ONTAK2010]Peaks_主席树+倍增+kruscal重构树 Description 在Bytemountains有N座山峰,每座山峰有他的高度h_i.有些山峰之间有双向道 ...

  4. [NOI2018]归程(kruscal重构树)

    [NOI2018]归程 题面太长辣,戳这里 模拟赛上写了一个spfa (关于spfa,它已经死了),然后一个st表水完暴力跑路.考后说是Kruscal重构树或者可持久化并查集???这都是些什么东西.不 ...

  5. $ CometOJ-Contest\#11\ D$ $Kruscal$重构树

    正解:$Kruscal$重构树 解题报告: 传送门$QwQ$ 发现一个图上搞就很麻烦,考虑变为生成树达到原有效果. 因为在询问的时候是要求走到的点编号尽量小,发现这个时候点的编号就成为限制了,于是不难 ...

  6. 洛谷$P4768\ [NOI2018]$归程 $kruscal$重构树

    正解:$kruscal$重构树 解题报告: 传送门$QwQ$ 语文不好选手没有人权$TT$连题目都看不懂真的要哭了$kk$ 所以先放个题目大意?就说给定一个$n$个点,$m$条边的图,每条边有长度和海 ...

  7. 学习笔记:Kruscal 重构树

    网上感觉没有什么很详细 + 证明的讲解啊) 前置:Kruskal 求最小生成树. 这个算法可以将一棵树 / 无向连通图重构成一颗有性质的新树. 算法可以解决一些树上瓶颈边权之类的问题,可以把需要持久化 ...

  8. 【题解】洛谷P1967 [NOIP2013TG] 货车运输(LCA+kruscal重构树)

    洛谷P1967:https://www.luogu.org/problemnew/show/P1967 思路 感觉2013年D1T3并不是非常难 但是蒟蒻还是WA了一次 从题目描述中看出每个点之间有许 ...

  9. BZOJ_5415_[Noi2018]归程_kruscal重构树+倍增+最短路

    BZOJ_5415_[Noi2018]归程_kruscal重构树+倍增 Description www.lydsy.com/JudgeOnline/upload/noi2018day1.pdf 好久不 ...

随机推荐

  1. 类中被final修饰的成员变量需要初始化

    类中被final修饰的成员变量需要初始化,否则编译不通过,因为final修饰后不能再赋值,因此必须初始化.

  2. 昂贵的聘礼 POJ - 1062(最短路)

    年轻的探险家来到了一个印第安部落里.在那里他和酋长的女儿相爱了,于是便向酋长去求亲.酋长要他用10000个金币作为聘礼才答应把女儿嫁给他.探险家拿不出这么多金币,便请求酋长降低要求.酋长说:" ...

  3. 基本排序算法(Java)

    基本排序算法 (Java) 经过几次笔试,发现自己的Java基础有些薄弱,基本的排序算法掌握的还不够熟练,需要多加学习总结. 1. 选择排序 思想: 给定一个整数数组,例 int[] a ={38,6 ...

  4. 「雅礼集训 2017 Day5」矩阵

    填坑填坑.. 感谢wwt耐心讲解啊.. 如果要看这篇题解建议从上往下读不要跳哦.. 30pts 把$A$和$C$看成$n$个$n$维向量,那$A_i$是否加入到$C_j$中就可以用$B_{i,j}$表 ...

  5. JS for循环 if判断、white循环。小练习

    1----输入正整数n,求1-n的和. var n=prompt("请输入一个正整数"); var sum=0; for (var i=1;i<=n;i++) { sum=s ...

  6. idea设置调用方法时提示方法注释

    如图所示:打开file-->setting-->Editor-->General,搜索show,然后勾选上Show quick documentation on mouse move ...

  7. 使用mongo获取文章

    先在命令行中输入mongo,进入mongo,然后 show dbs # 从结果中发现有cmb_demo_23_hacker use cmb_demo_23_hacker db.all_in_one.f ...

  8. H5 _拖放使用

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  9. 马昕璐 201771010118《面向对象程序设计(java)》第十八周学习总结

    实验十八  总复习 实验时间 2018-12-30 1.实验目的与要求 (1) 综合掌握java基本程序结构: (2) 综合掌握java面向对象程序设计特点: (3) 综合掌握java GUI 程序设 ...

  10. Qt5+MSVC2015环境将VS2015编写的控制台程序转化为GUI程序

    如题所述,如何将VS2015编写的控制台程序转化为Qt5+MSVC2015环境编译的GUI程序? 最近想到这个操作,类似于Linux下使用的命令行操作转到Windows下使用GUI操作,看了控制台的命 ...