题目连接:https://www.luogu.org/problemnew/show/U60884

题意:有N个点,标号为1∼N,用N−1条双向带权通道连接,保证任意两个点能互相到达。

Q次询问,问从编号为x的点到达标号L∼R的点其中一个点的最小距离是多少。

说明 :N,Q<1e5,边权<1e4;

思路:不难想到点分树,保存每个点到其“负责”的点的距离,这样的话可以套线段树,线段树保存其他点到点的距离。

但是,点分树上有个需要解决的问题是:如果x顺着点分树向父亲走,那么在父亲保存的线段树中要除去从儿子上来的那一部分(否则的话,不是简单路径,求出来的可能会错)。 想了一下这里很难实现,因为是取min操作,所以线段树上二分估计也不行。 所以卡住了。

然后求问群友,群友提到了虚树,我觉得此题的数据需要没法实现。  后面猛地想通,我们不需要考虑“回走”这种情况,因为我们求的是最小距离,而非简单路径肯定是没有简单路径优的。   那么就是果然如标题所说,是个板子题了。

代码实现:点分树+线段树+树剖求LCA。

点分树:点分治的过程中新建的树,树根是第一次找到是重心,每一层的重心与上一层的重心连边得到点分树。

点分树里保存的是自己作为重心时,会“负责”的点,即此时还没有vis过的,且与自己连通的点。

所有经过“x”到达的点,保存在了两部分里:一是x在点分树里保存的信息。二是x在点分树的祖先里保存的信息(这一部分保留了x向上传递是信息,所以大部分题要把这里除去,此题由于是取min,可以不考虑)。

求LCA:开始用了ST表,但是感觉空间耗费太大,就改为了树剖,跑起来还挺快。

#include<bits/stdc++.h>
#define rep(i,w,v) for(int i=w;i<=v;i++)
#define FOR() for(int i=Laxt[u];i;i=Next[i])
using namespace std;
#define ll long long
#define RG register
#define maxn 100010
#define maxm 200010
#define inf 1e9
int fa[maxn],n,cnt;
int dep[maxn];bool vis[maxn];
int Laxt[maxn],Next[maxm],To[maxm],Len[maxm];
int siz[maxn],fcy[maxn],hson[maxn],Top[maxn];
void add(int u,int v,int len)
{
Next[++cnt]=Laxt[u]; Laxt[u]=cnt; To[cnt]=v; Len[cnt]=len;
}
void dfs1(int u,int ff)
{
fa[u]=ff; siz[u]=;
dep[u]=dep[ff]+;
FOR(){
int v=To[i];if(v==ff)continue;
fcy[v]=fcy[u]+Len[i];
dfs1(v,u); siz[u]+=siz[v];
if(siz[v]>siz[hson[u]]) hson[u]=v;
}
}
void dfs2(int u,int tp)
{
Top[u]=tp;
if(hson[u])dfs2(hson[u],tp);
FOR()
if(To[i]!=fa[u]&&To[i]!=hson[u])
dfs2(To[i],To[i]);
}
int LCA(int u,int v)
{
while(Top[u]^Top[v]) dep[Top[u]]<dep[Top[v]]?v=fa[Top[v]]:u=fa[Top[u]];
return dep[u]<dep[v]?u:v;
}
int Dis(int u,int v){return fcy[u]+fcy[v]-*fcy[LCA(u,v)];}
int Fa[maxn],Size,root,mx,rt[maxn];
void Getroot(int u,int ff)
{
siz[u]=;int ret=;
FOR(){
int v=To[i];if(v==ff||vis[v])continue;
Getroot(v,u);siz[u]+=siz[v];
ret=max(ret,siz[v]);
}
ret=max(ret,Size-siz[u]);
if(ret<mx) mx=ret,root=u;
}
void DFS(int u,int ff)
{
vis[u]=true;Fa[u]=ff;
FOR(){
int v=To[i];if(vis[v])continue;
mx=Size=siz[v];
Getroot(v,u);
DFS(root,u);
}
}
struct in{ int l,r,mn; }s[maxn<<]; int tot;
void modify(int &Now,int L,int R,int pos,int val)
{
if(!Now){Now=++tot; s[tot].mn=inf;}
s[Now].mn=min(s[Now].mn,val);
if(L==R) return; int Mid=(L+R)>>;
if(pos<=Mid) modify(s[Now].l,L,Mid,pos,val);
else modify(s[Now].r,Mid+,R,pos,val);
}
int query(int Now,int L,int R,int l,int r)
{
if(!Now) return inf;
if(l<=L&&r>=R) return s[Now].mn;
int Mid=(L+R)>>,res=inf;
if(l<=Mid) res=min(res,query(s[Now].l,L,Mid,l,r));
if(r>Mid) res=min(res,query(s[Now].r,Mid+,R,l,r));
return res;
}
void Modify(int x)
{
modify(rt[x],,n,x,);
for(int i=x;Fa[i];i=Fa[i]){
int dis=Dis(x,Fa[i]);
modify(rt[Fa[i]],,n,x,dis);
}
}
int Query(int x,int L,int R)
{
int res=query(rt[x],,n,L,R);
for(int i=x;Fa[i];i=Fa[i])
{
int dis=Dis(x,Fa[i]);
res=min(res,dis+query(rt[Fa[i]],,n,L,R));
}
return res;
}
int main()
{
int Q,u,v,len;
scanf("%d",&n);
rep(i,,n-){
scanf("%d%d%d",&u,&v,&len);
add(u,v,len); add(v,u,len);
}
dfs1(,); dfs2(,);
Size=mx=n;
Getroot(,); DFS(root,);
rep(i,,n) Modify(i);
int ans=,l,r,x;
scanf("%d",&Q);
while(Q--){
scanf("%d%d%d",&l,&r,&x);
printf("%d\n",Query(x,l,r));
}
return ;
}

luoguU60884 【模板】动态点分治套线段树的更多相关文章

  1. BZOJ4317Atm的树&BZOJ2051A Problem For Fun&BZOJ2117[2010国家集训队]Crash的旅游计划——二分答案+动态点分治(点分树套线段树/点分树+vector)

    题目描述 Atm有一段时间在虐qtree的题目,于是,他满脑子都是tree,tree,tree…… 于是,一天晚上他梦到自己被关在了一个有根树中,每条路径都有边权,一个神秘的声音告诉他,每个点到其他的 ...

  2. [APIO2019] [LOJ 3146] 路灯 (cdq分治或树状数组套线段树)

    [APIO2019] [LOJ 3146] 路灯 (cdq分治或树状数组套线段树) 题面 略 分析 首先把一组询问(x,y)看成二维平面上的一个点,我们想办法用数据结构维护这个二维平面(注意根据题意这 ...

  3. BZOJ.4553.[HEOI2016&TJOI2016]序列(DP 树状数组套线段树/二维线段树(MLE) 动态开点)

    题目链接:BZOJ 洛谷 \(O(n^2)\)DP很好写,对于当前的i从之前满足条件的j中选一个最大值,\(dp[i]=d[j]+1\) for(int j=1; j<i; ++j) if(a[ ...

  4. P3157 [CQOI2011]动态逆序对(树状数组套线段树)

    P3157 [CQOI2011]动态逆序对 树状数组套线段树 静态逆序对咋做?树状数组(别管归并QWQ) 然鹅动态的咋做? 我们考虑每次删除一个元素. 减去的就是与这个元素有关的逆序对数,介个可以预处 ...

  5. [TJOI2017][bzoj4889] 不勤劳的图书管理员 [线段树套线段树]

    题面 传送门 思路 考虑两本书的位置交换对答案的贡献: (为了方便描述,用"左边那本"和"右边那本"称呼两本我们要交换的书,"中间那本"是我 ...

  6. 【BZOJ3295】动态逆序对(线段树,树状数组)

    [BZOJ3295]动态逆序对(线段树,树状数组) 题面 Description 对于序列A,它的逆序对数定义为满足iAj的数对(i,j)的个数.给1到n的一个排列,按照某种顺序依次删除m个元素,你的 ...

  7. 【BZOJ2001】[HNOI2010]城市建设(CDQ分治,线段树分治)

    [BZOJ2001][HNOI2010]城市建设(CDQ分治,线段树分治) 题面 BZOJ 洛谷 题解 好神仙啊这题.原来想做一直不会做(然而YCB神仙早就切了),今天来怒写一发. 很明显这个玩意换种 ...

  8. BZOJ 1901 Zju2112 Dynamic Rankings 树状数组套线段树

    题意概述:带修改求区间第k大. 分析: 我们知道不带修改的时候直接上主席树就可以了对吧?两个版本号里面的节点一起走在线段树上二分,复杂度是O((N+M)logN). 然而这里可以修改,主席树显然是凉了 ...

  9. bzoj3196 二逼平衡树 树状数组套线段树

    题目传送门 思路:树状数组套线段树模板题. 什么是树状数组套线段树,普通的树状数组每个点都是一个权值,而这里的树状数组每个点都是一颗权值线段树,我们用前缀差分的方法求得每个区间的各种信息, 其实关键就 ...

随机推荐

  1. 1 RAID技术入门

    序   RAID一页通整理所有RAID技术.原理并配合相应RAID图解,给所有存储新人提供一个迅速学习.理解RAID技术的网上资源库,本文将持续更新,欢迎大家补充及投稿.中国存储网一如既往为广大存储界 ...

  2. .NET Core 之 Nancy 基本使用

    Nancy简介 Nancy是一个轻量级的独立的框架,下面是官网的一些介绍: Nancy 是一个轻量级用于构建基于 HTTP 的 Web 服务,基于 .NET 和 Mono 平台,框架的目标是保持尽可能 ...

  3. JVM的基本结构和JVM的内存结构

    这里概要介绍一下JVM在启动后,作为操作系统的一个进程的基本结构,以及从操作系统角度看,JVM如何管理它从操作系统里申请来的内存的,也就是JVM的内存结构或者叫JVM内存模型. 1.JVM的基本结构 ...

  4. charles 4.2.1 Ubuntu破解版安装

    charles 4.2.1 Ubuntu破解版安装 下载 charles-proxy-4.2.1_amd64.tar.gz 破解版 charles.jar 破解包 解压 sudo tar -zxvf ...

  5. 【转帖】Intel AMD 龙芯2019年12月份最新产品线

    Intel未来三代U集体曝光:14nm退回去了! https://news.cnblogs.com/n/651244/ 不过没搞懂 为啥中芯国际已经开始量产14nm了 龙芯为什么不用.. 3A4000 ...

  6. Go语言【数据结构】字符串

    字符串 简介 一个字符串是一个不可改变的字节序列,字符串通常是用来包含人类可读的文本数据.和数组不同的是,字符串的元素不可修改,是一个只读的字节数组.每个字符串的长度虽然也是固定的,但是字符串的长度并 ...

  7. C#使用CSS选择器抓取页面内容

    最近在查wpf绘图资料时,偶然看到Python使用CSS选择器抓取网页的功能.觉得很强,这里用C#也实现一下. 先介绍一下CSS选择器 在 CSS 中,选择器是一种模式,用于选择需要添加样式的元素. ...

  8. 如何定位 golang 进程 hang 死的 bug

    之前在 golang 群里有人问过为什么程序会莫名其妙的 hang 死然后不再响应任何请求.单核 cpu 打满. 这个特征和我们公司的某个系统曾经遇到的情况很相似,内部经过了很长时间的定位分析总结,期 ...

  9. 网页调用文件另存为js

    查看引用是否正常,页面添加html代码. <a id="downLoad" onclick="oDownLoad('downLoad')">下载&l ...

  10. JBOOT使用总结

    @Override public SwAdmin findById(long id) { return DAO.findFirst("SELECT * FROM sw_admin WHERE ...