【bzoj4012】 HNOI2015—开店
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—开店的更多相关文章
- [BZOJ4012][HNOI2015]开店(动态点分治,树链剖分)
4012: [HNOI2015]开店 Time Limit: 70 Sec Memory Limit: 512 MBSubmit: 2168 Solved: 947[Submit][Status] ...
- BZOJ4012 [HNOI2015]开店
Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现 ...
- BZOJ4012[HNOI2015]开店——树链剖分+可持久化线段树/动态点分治+vector
题目描述 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现她们面临着一个 ...
- BZOJ4012 HNOI2015开店(树链剖分+主席树)
考虑这样一个问题:一棵树初始全是白点,有两种操作:把一个点染黑:询问某点到所有黑点的距离之和. 注意到树上两点x和y的距离为depth[x]+depth[y]-depth[lca(x,y)]*2.要求 ...
- BZOJ4012: [HNOI2015]开店【动态点分治】
Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现 ...
- BZOJ4012 [HNOI2015]开店 (动态点分治)
Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现 ...
- BZOJ4012 [HNOI2015]开店 【动态点分治 + splay】
题目链接 BZOJ4012 题解 Mychael并没有A掉,而是T掉了 讲讲主要思路 在点分树上每个点开两棵\(splay\), 平衡树\(A\)维护子树中各年龄到根的距离 平衡树\(B\)维护子树中 ...
- 【BZOJ4012】[HNOI2015]开店 动态树分治+二分
[BZOJ4012][HNOI2015]开店 Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点 ...
- 【BZOJ4012】开店(主席树)
[BZOJ4012]开店(主席树) 题面 Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱 ...
- [HNOI2015]开店 树链剖分,主席树
[HNOI2015]开店 LG传送门 蒟蒻表示不会动态淀粉质. 先把点按年龄排序, 设\(dis[i]\)表示\(i\)到根的距离. 把我们要算的东西稍微变下形:\(ans\) \[ = \sum \ ...
随机推荐
- 20155301 Exp7 网络欺诈防范
20155301 Exp7 网络欺诈防范 1.基础问题回答 (1)通常在什么场景下容易受到DNS spoof攻击 (2)在日常生活工作中如何防范以上两攻击方法 2.实践过程记录 简单应用SET工具建立 ...
- 【来龙去脉系列】AutoMapper一款自动映射框架
前言 通常在一个应用程序中,我们开发人员会在两个不同的类型对象之间传输数据,通常我们会用DTOs(数据传输对象),View Models(视图模型),或者直接是一些从一个service或者Web AP ...
- HTTPUTILS
maven依赖 <dependency> <groupId>org.apache.httpcomponents</groupId> <artifactId&g ...
- lnmp如何实现伪静态,默认目录伪静态
LNMP一键安装包一直是我小内存VPS的首选安装环境,因为它占用资源少,性能高.最新发布的lnmp0.7还增加了对LNMPA的支持,LNMPA使用Nginx作为前端服务能够更快更及时的静态页面.js. ...
- OpenCV实战:人脸关键点检测(FaceMark)
Summary:利用OpenCV中的LBF算法进行人脸关键点检测(Facial Landmark Detection) Author: Amusi Date: 2018-03-20 ...
- dubbo见解
调用关系说明 服务容器负责启动,加载,运行服务提供者. 服务提供者在启动时,向注册中心注册自己提供的服务. 服务消费者在启动时,向注册中心订阅自己所需的服务. 注册中心返回服务提供者地址列表给消费者, ...
- FastDFS教程Ⅰ-文件服务器安装与Nginx配置
1.简介 FastDFS是一个开源的轻量级分布式文件系统,它对文件进行管理,功能包括:文件存储.文件同步.文件访问(文件上传.文件下载)等,解决了大容量存储和负载均衡的问题.特别适合以文件为载 ...
- Unity插件-NGUI学习笔记
Anchors 的作用 类似Android里面的.9格式图片的功能, 边框可以随着文字变大而变大. 实现方法: 用NGUI 创建一个Sprite, 命名为TextBg, 一个Label, Label的 ...
- mac 安装 tomcat 配置
前面的话:记录下 Mac 安装配置 Tomcat 过程 1. 下载安装 Tomcat 下载 Tomcat 地址(官方地址):https://tomcat.apache.org/download-80. ...
- Alpha Version Release Of Teamwork: Appendix 1 BUG BASH
在为期一周的发布周中,我们将app本身最后的细节完善,功能代码到位,UI不断改进和优化,团队在开始准备发布之前,对整个APP进行了一次BUG检查,每个人都部署了app在自己的android设备上进行测 ...