http://www.lydsy.com/JudgeOnline/problem.php?id=4012 (题目链接)

题意

  一棵树,每条边有正边权,每个点的点权非负。若干组询问,强制在线,每次查询点权在范围${[L,R]}$之间的点到某一点${U}$的距离和。

Solution

  这道题做法很多啊。动态树分治。

  按照套路,我们存下每个重心的子树到这个重心的距离之和,以及每个重心的子树到这个重心的父亲的距离之和。这两个东西按照点权排序后用一个数组存起来并处理成前缀和,然后每次询问区间${[L,R]}$,我们只需要在子树中二分查找即可。

细节

  样例很强啊,过了样例即可AC。当然记得开LL

代码

// bzoj4012
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<vector>
#include<cmath>
#define LL long long
#define inf 1ll<<30
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std; const int maxn=200010;
int n,head[maxn],a[maxn],par[maxn],len[maxn],bin[30];
int L,R,U,Q,A; struct edge {int to,next,w;}e[maxn<<1];
struct data {
int num;LL w;
friend bool operator < (const data a,const data b) {
return a.num<b.num;
}
};
vector<data> c[maxn],t[maxn]; namespace LittleTrick {
int fa[maxn][30],deep[maxn],cnt;
LL d[maxn]; void link(int u,int v,int w) {
e[++cnt]=(edge){v,head[u],w};head[u]=cnt;
e[++cnt]=(edge){u,head[v],w};head[v]=cnt;
}
void dfs(int x) {
for (int i=1;i<=20;i++) fa[x][i]=fa[fa[x][i-1]][i-1];
for (int i=head[x];i;i=e[i].next) if (e[i].to!=fa[x][0]) {
deep[e[i].to]=deep[x]+1;
d[e[i].to]=d[x]+e[i].w;
fa[e[i].to][0]=x;
dfs(e[i].to);
}
}
int lca(int x,int y) {
if (deep[x]<deep[y]) swap(x,y);
int t=deep[x]-deep[y];
for (int i=0;bin[i]<=t;i++) if (bin[i]&t) x=fa[x][i];
for (int i=20;i>=0;i--) if (fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
return x==y ? x : fa[x][0];
}
LL dis(int x,int y) {
return d[x]+d[y]-2*d[lca(x,y)];
}
}
using namespace LittleTrick; namespace NodeDivide {
int size[maxn],f[maxn],vis[maxn],sum,Dargen; void caldargen(int x,int fa) {
size[x]=1,f[x]=0;
for (int i=head[x];i;i=e[i].next) if (e[i].to!=fa && !vis[e[i].to]) {
caldargen(e[i].to,x);
size[x]+=size[e[i].to];
f[x]=max(f[x],size[e[i].to]);
}
f[x]=max(f[x],sum-size[x]);
if (f[x]<f[Dargen]) Dargen=x;
}
void caldeep(int x,int fa,int p) {
c[p].push_back((data){a[x],dis(x,p)});
t[p].push_back((data){a[x],par[p] ? dis(x,par[p]) : 0});
for (int i=head[x];i;i=e[i].next)
if (e[i].to!=fa && !vis[e[i].to]) caldeep(e[i].to,x,p);
}
void work(int x) {
vis[x]=1;
t[x].push_back((data){-1,0});t[x].push_back((data){inf,0}); //方便二分查找,不用特判了
c[x].push_back((data){-1,0});c[x].push_back((data){inf,0});
caldeep(x,0,x);
for (int i=head[x];i;i=e[i].next) if (!vis[e[i].to]) {
Dargen=0;sum=size[e[i].to];
caldargen(e[i].to,x);
par[Dargen]=x;
work(Dargen);
}
sort(t[x].begin(),t[x].end());len[x]=t[x].size();
sort(c[x].begin(),c[x].end());len[x]=c[x].size();
for (int i=1;i<len[x];i++) t[x][i].w+=t[x][i-1].w;
for (int i=1;i<len[x];i++) c[x][i].w+=c[x][i-1].w;
}
void Init() {
f[Dargen=0]=inf;sum=n;
caldargen(1,0);
work(Dargen);
}
}
using namespace NodeDivide; namespace Query {
int Lower_bound(int x) {
int l=0,r=len[x],res;
while (l<=r) {
int mid=(l+r)>>1;
if (c[x][mid].num<L) res=mid,l=mid+1;
else r=mid-1;
}
return res;
}
int Upper_bound(int x) {
int l=0,r=len[x],res;
while (l<=r) {
int mid=(l+r)>>1;
if (c[x][mid].num<=R) res=mid,l=mid+1;
else r=mid-1;
}
return res;
}
LL query(int x,int tag) {
LL d1=dis(x,U),res=0;
LL l=Lower_bound(x),r=Upper_bound(x);
res+=(c[x][r].w-c[x][l].w)+(r-l-tag)*d1;
res-=(t[x][r].w-t[x][l].w);
if (!par[x]) return res;
return res+query(par[x],r-l);
}
}
using namespace Query; int main() {
bin[0]=1;for (int i=1;i<=20;i++) bin[i]=bin[i-1]<<1;
scanf("%d%d%d",&n,&Q,&A);
for (int i=1;i<=n;i++) scanf("%d",&a[i]);
for (int u,v,w,i=1;i<n;i++) {
scanf("%d%d%d",&u,&v,&w);
link(u,v,w);
}
dfs(1);
Init();
LL ans=0;
for (int a,b,i=1;i<=Q;i++) {
scanf("%d%d%d",&U,&a,&b);
if (i==1) L=min(a%A,b%A),R=max(a%A,b%A);
else L=min((a+ans)%A,(b+ans)%A),R=max((a+ans)%A,(b+ans)%A);
printf("%lld\n",ans=query(U,0));
}
return 0;
}

【bzoj4012】 HNOI2015—开店的更多相关文章

  1. [BZOJ4012][HNOI2015]开店(动态点分治,树链剖分)

    4012: [HNOI2015]开店 Time Limit: 70 Sec  Memory Limit: 512 MBSubmit: 2168  Solved: 947[Submit][Status] ...

  2. BZOJ4012 [HNOI2015]开店

    Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现 ...

  3. BZOJ4012[HNOI2015]开店——树链剖分+可持久化线段树/动态点分治+vector

    题目描述 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现她们面临着一个 ...

  4. BZOJ4012 HNOI2015开店(树链剖分+主席树)

    考虑这样一个问题:一棵树初始全是白点,有两种操作:把一个点染黑:询问某点到所有黑点的距离之和. 注意到树上两点x和y的距离为depth[x]+depth[y]-depth[lca(x,y)]*2.要求 ...

  5. BZOJ4012: [HNOI2015]开店【动态点分治】

    Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现 ...

  6. BZOJ4012 [HNOI2015]开店 (动态点分治)

    Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现 ...

  7. BZOJ4012 [HNOI2015]开店 【动态点分治 + splay】

    题目链接 BZOJ4012 题解 Mychael并没有A掉,而是T掉了 讲讲主要思路 在点分树上每个点开两棵\(splay\), 平衡树\(A\)维护子树中各年龄到根的距离 平衡树\(B\)维护子树中 ...

  8. 【BZOJ4012】[HNOI2015]开店 动态树分治+二分

    [BZOJ4012][HNOI2015]开店 Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点 ...

  9. 【BZOJ4012】开店(主席树)

    [BZOJ4012]开店(主席树) 题面 Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱 ...

  10. [HNOI2015]开店 树链剖分,主席树

    [HNOI2015]开店 LG传送门 蒟蒻表示不会动态淀粉质. 先把点按年龄排序, 设\(dis[i]\)表示\(i\)到根的距离. 把我们要算的东西稍微变下形:\(ans\) \[ = \sum \ ...

随机推荐

  1. Android Studio常用快捷键 - 转

    Android Studio常用快捷键 1. Ctrl+D: 集合了复制和粘贴两个操作,如果有选中的部分就复制选中的部分,并在选中部分的后面粘贴出来,如果没有选中的部分,就复制光标所在的行,并在此行的 ...

  2. 带Alpha通道的色彩叠加问题

    css3的rgba色彩模式.png/gif图片的alpha通道.canvas的rgba色彩模式.css3的阴影.css3的opacity属性等等,这些应用在网页中,有意无意间,我们的页面多了许多半透明 ...

  3. Markdown 编辑器

    桌面编辑器 MarkdownPad Pro 版注册邮箱: Soar360@live.com 授权密钥: GBPduHjWfJU1mZqcPM3BikjYKF6xKhlKIys3i1MU2eJHqWGI ...

  4. JNI探秘-----你不知道的FileInputStream的秘密

    作者:zuoxiaolong8810(左潇龙),转载请注明出处,特别说明:本博文来自博主原博客,为保证新博客中博文的完整性,特复制到此留存,如需转载请注明新博客地址即可. 设计模式系列结束,迎来了LZ ...

  5. 批处理修改Hosts文件

    @echo off set stHosts=127.0.0.1 www.XXX.com FOR /F "eol=# tokens=1 delims=" %%i in (%Syste ...

  6. [dx11]利用SpriteFont绘制中文--本地化文本

    1.下载DirectX11 Tool Kit SDK,解压后编译,生成DirectXTK.lib库文件和MakeSpriteFont应用工具; 2.在Dx11环境基础上,用生成的库文件搭建XTK环境; ...

  7. Python3.7 + jupyter安装(CentOS6.5)

    Python3.7 + jupyter安装(CentOS6.5) 方法一(anaconda): anaconda是一个开源的Python发行版本 包含conda,python等大量的科学包以及依赖 优 ...

  8. k8s之使用secret获取私有仓库镜像

    一.前言 其实这次实践算不上特别复杂,只是在实践过程中遇到了一些坑,以及填坑的方法是非常值得在以后的学习过程中参考借鉴的 二.知识准备 1.harbor是一个企业级的镜像仓库,它比起docker re ...

  9. Python进阶量化交易场外篇4——寻找最优化策略参数

    新年伊始,很荣幸笔者的<教你用 Python 进阶量化交易>专栏在慕课专栏板块上线了,欢迎大家订阅!为了能够提供给大家更轻松的学习过程,笔者在专栏内容之外会陆续推出一些手记来辅助同学们学习 ...

  10. C语言版本:循环单链表的实现

    SClist.h #ifndef __SCLIST_H__ #define __SCLIST_H__ #include<cstdio> #include<malloc.h> # ...