bzoj 4012 [HNOI2015]开店 【树链剖分】

最新推荐文章于 2020-02-21 11:23:53 发布
最新推荐文章于 2020-02-21 11:23:53 发布
阅读量834 收藏
点赞数
分类专栏： 树链剖分 文章标签： 树链剖分 线段树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/heheda_is_an_OIer/article/details/50088385
版权

这道题好像各位神犇都是用动态树分治写的，我这么弱只能用树剖水一水了。

dis(x,y)=dis(1,x)+dis(1,y)-dis(1,lca(x,y))

前两个都是定值，求第三项就行了。

每个点x维护不在重链上的子节点的个数*dis(1,x)，查询显然，再yy一些别的东西就可以了。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>


#define ll long long
#define inf 1e9
#define eps 1e-10
#define md
#define N 150010
#define TR 20000000
using namespace std;
struct yts { int x,t,l,ne;} e[2*N];
int n,td;
int num,dfs_cnt,cnt;
int dat[N],v[N],dep[N],sz[N],son[N],fa[N][20],pos[N],top[N],low[N],right[N],dis[N],root[N],s_sz[N];
ll s_dis[N];
ll sum[TR];
int size[TR],ch[TR][2],qs[TR];
struct AA { int age,pos;} a[N];
bool cmp(AA a,AA b) { return a.age<b.age;}


void lisanhua()
{
<span style="white-space:pre">	</span>sort(a+1,a+n+1,cmp);
<span style="white-space:pre">	</span>a[0].age=-1; int now=-1; dat[0]=-1;
<span style="white-space:pre">	</span>for (int i=1;i<=n;i++)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>if (a[i].age!=now) { dat[++td]=a[i].age; now=a[i].age;}
<span style="white-space:pre">		</span>a[i].age=td;
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>dat[td+1]=inf;
<span style="white-space:pre">	</span>//for (int i=1;i<=td;i++) printf("%d ",dat[i]); printf("\n");
}
void put(int x,int y,int l)
{
<span style="white-space:pre">	</span>num++; e[num].x=x; e[num].t=y; e[num].l=l;
<span style="white-space:pre">	</span>e[num].ne=v[x]; v[x]=num;
}
void dfs1(int x)
{
<span style="white-space:pre">	</span>sz[x]=1; son[x]=0;
<span style="white-space:pre">	</span>for (int i=v[x];i;i=e[i].ne)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>int y=e[i].t;
<span style="white-space:pre">		</span>if (y!=fa[x][0])
<span style="white-space:pre">		</span>{
<span style="white-space:pre">			</span>fa[y][0]=x; dep[y]=dep[x]+1; dis[y]=dis[x]+e[i].l;
<span style="white-space:pre">			</span>dfs1(y);
<span style="white-space:pre">			</span>sz[x]+=sz[y]; if (sz[y]>sz[son[x]]) son[x]=y;
<span style="white-space:pre">		</span>}
<span style="white-space:pre">	</span>}
}


void dfs2(int x,int f)
{
<span style="white-space:pre">	</span>pos[x]=++dfs_cnt; top[x]=f;
<span style="white-space:pre">	</span>if (son[x])
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>dfs2(son[x],f);
<span style="white-space:pre">		</span>low[x]=low[son[x]];
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>else low[x]=pos[x];
<span style="white-space:pre">	</span>for (int i=v[x];i;i=e[i].ne)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>int y=e[i].t;
<span style="white-space:pre">		</span>if (y!=fa[x][0]&&y!=son[x])
<span style="white-space:pre">		</span>  dfs2(y,y);
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>right[x]=dfs_cnt;
}


int find1(int x)//最小的大于等于x的。
{
<span style="white-space:pre">	</span>int l=1,r=td+1;
<span style="white-space:pre">	</span>while (l!=r)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>int mid=(l+r)>>1;
<span style="white-space:pre">		</span>if (dat[mid]>=x) r=mid; else l=mid+1;
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>return l;
}


int find2(int x)//最大的小于等于x的。
{
<span style="white-space:pre">	</span>int l=0,r=td;
<span style="white-space:pre">	</span>while (l!=r)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>int mid=(l+r+1)>>1;
<span style="white-space:pre">		</span>if (dat[mid]<=x) l=mid; else r=mid-1;
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>return l;
}


void build_lca()
{
<span style="white-space:pre">	</span>for (int j=1;j<=18;j++)
<span style="white-space:pre">	</span>  for (int i=1;i<=n;i++)
<span style="white-space:pre">	</span>    if (fa[i][j-1]) fa[i][j]=fa[fa[i][j-1]][j-1];
}


int lca(int x,int y)
{
<span style="white-space:pre">	</span>if (dep[x]<dep[y]) swap(x,y); 
<span style="white-space:pre">	</span>int t=dep[x]-dep[y];
<span style="white-space:pre">	</span>for (int i=18;i>=0;i--)
<span style="white-space:pre">	</span>  if (t&(1<<i)) x=fa[x][i];
<span style="white-space:pre">	</span>if (x==y) return x;
<span style="white-space:pre">	</span>for (int i=18;i>=0;i--)
<span style="white-space:pre">	</span>  if (fa[x][i]!=fa[y][i])
<span style="white-space:pre">	</span>    x=fa[x][i],y=fa[y][i];
<span style="white-space:pre">	</span>return fa[x][0];
}


void insert_sz(int &i,int pre,int l,int r,int x)
{
<span style="white-space:pre">	</span>i=++cnt;
<span style="white-space:pre">	</span>size[i]=size[pre]+1; sum[i]=sum[pre]; qs[i]=qs[pre];
<span style="white-space:pre">	</span>if (l==r) return;
<span style="white-space:pre">	</span>int mid=(l+r)>>1;
<span style="white-space:pre">	</span>if (x<=mid)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>insert_sz(ch[i][0],ch[pre][0],l,mid,x);
<span style="white-space:pre">		</span>ch[i][1]=ch[pre][1];
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>else
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>insert_sz(ch[i][1],ch[pre][1],mid+1,r,x);
<span style="white-space:pre">		</span>ch[i][0]=ch[pre][0];
<span style="white-space:pre">	</span>}
}


void insert(int &i,int pre,int l,int r,int x,int dis)
{
<span style="white-space:pre">	</span>i=++cnt;
<span style="white-space:pre">	</span>sum[i]=sum[pre]+dis; size[i]=size[pre]; qs[i]=qs[pre]+1;
<span style="white-space:pre">	</span>if (l==r) return;
<span style="white-space:pre">	</span>int mid=(l+r)>>1;
<span style="white-space:pre">	</span>if (x<=mid)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>insert(ch[i][0],ch[pre][0],l,mid,x,dis);
<span style="white-space:pre">		</span>ch[i][1]=ch[pre][1];
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>else
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>insert(ch[i][1],ch[pre][1],mid+1,r,x,dis);
<span style="white-space:pre">		</span>ch[i][0]=ch[pre][0];
<span style="white-space:pre">	</span>}
}
<span style="white-space:pre">	</span>
ll query_sz(int L,int R,int l,int r,int ql,int qr)
{
<span style="white-space:pre">	</span>if (ql<=l&&r<=qr) return size[R]-size[L];
<span style="white-space:pre">	</span>int mid=(l+r)>>1;
<span style="white-space:pre">	</span>if (qr<=mid) return query_sz(ch[L][0],ch[R][0],l,mid,ql,qr);
<span style="white-space:pre">	</span>if (mid+1<=ql) return query_sz(ch[L][1],ch[R][1],mid+1,r,ql,qr);
<span style="white-space:pre">	</span>return query_sz(ch[L][0],ch[R][0],l,mid,ql,qr)+query_sz(ch[L][1],ch[R][1],mid+1,r,ql,qr);
}


ll query(int L,int R,int l,int r,int ql,int qr)
{
<span style="white-space:pre">	</span>if (ql<=l&&r<=qr) return sum[R]-sum[L];
<span style="white-space:pre">	</span>int mid=(l+r)>>1;
<span style="white-space:pre">	</span>if (qr<=mid) return query(ch[L][0],ch[R][0],l,mid,ql,qr);
<span style="white-space:pre">	</span>if (mid+1<=ql) return query(ch[L][1],ch[R][1],mid+1,r,ql,qr);
<span style="white-space:pre">	</span>return query(ch[L][0],ch[R][0],l,mid,ql,qr)+query(ch[L][1],ch[R][1],mid+1,r,ql,qr);
}


ll query_line(int L,int R,int l,int r,int ql,int qr)
{
<span style="white-space:pre">	</span>if (ql<=l&&r<=qr) return qs[R]-qs[L];
<span style="white-space:pre">	</span>int mid=(l+r)>>1;
<span style="white-space:pre">	</span>if (qr<=mid) return query_line(ch[L][0],ch[R][0],l,mid,ql,qr);
<span style="white-space:pre">	</span>if (mid+1<=ql) return query_line(ch[L][1],ch[R][1],mid+1,r,ql,qr);
<span style="white-space:pre">	</span>return query_line(ch[L][0],ch[R][0],l,mid,ql,qr)+query_line(ch[L][1],ch[R][1],mid+1,r,ql,qr);
}


int main()
{
#ifndef ONLINE_JUDGE
<span style="white-space:pre">	</span>freopen("data.in","r",stdin); freopen("b.out","w",stdout);
#endif
<span style="white-space:pre">	</span>int Q,A;
<span style="white-space:pre">	</span>scanf("%d%d%d",&n,&Q,&A);
<span style="white-space:pre">	</span>for (int i=1;i<=n;i++)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>scanf("%d",&a[i].age);
<span style="white-space:pre">		</span>a[i].pos=i;
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>for (int i=1;i<n;i++)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>int x,y,l;
<span style="white-space:pre">		</span>scanf("%d%d%d",&x,&y,&l);
<span style="white-space:pre">		</span>put(x,y,l); put(y,x,l);
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>dep[1]=1; dfs1(1); dfs2(1,1); build_lca(); <span style="white-space:pre">	</span>
<span style="white-space:pre">	</span>lisanhua();
<span style="white-space:pre">	</span>int pre=0;
<span style="white-space:pre">	</span>for (int i=1;i<=n;i++)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>int x=a[i].pos;
<span style="white-space:pre">		</span>insert_sz(root[a[i].age],pre,1,n,pos[x]); pre=root[a[i].age];
<span style="white-space:pre">		</span>s_sz[a[i].age]++; s_dis[a[i].age]+=dis[x];
<span style="white-space:pre">		</span>while (x)
<span style="white-space:pre">		</span>{
<span style="white-space:pre">			</span>insert(root[a[i].age],pre,1,n,pos[x],dis[x]); pre=root[a[i].age];
<span style="white-space:pre">			</span>x=fa[top[x]][0];
<span style="white-space:pre">		</span>}
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>for (int i=1;i<=td;i++) s_sz[i]+=s_sz[i-1],s_dis[i]+=s_dis[i-1];
<span style="white-space:pre">	</span>ll ans=0;
<span style="white-space:pre">	</span>for (int i=1;i<=Q;i++)
<span style="white-space:pre">	</span>{
<span style="white-space:pre">		</span>int u,x,y;
<span style="white-space:pre">		</span>scanf("%d%d%d",&u,&x,&y);
<span style="white-space:pre">		</span>int L=(ans+x)%A,R=(ans+y)%A;
<span style="white-space:pre">		</span>if (L>R) swap(L,R);
<span style="white-space:pre">		</span>L=find1(L); R=find2(R);
<span style="white-space:pre">		</span>if (L>td||R<1||L>R)
<span style="white-space:pre">		</span>{
<span style="white-space:pre">			</span>ans=0; printf("0\n");
<span style="white-space:pre">			</span>continue;
<span style="white-space:pre">		</span>}
<span style="white-space:pre">		</span>int rtL=root[L-1],rtR=root[R];
<span style="white-space:pre">		</span>x=u;
<span style="white-space:pre">		</span>ans=((ll)dis[x])*(s_sz[R]-s_sz[L-1])+s_dis[R]-s_dis[L-1];
<span style="white-space:pre">		</span>ll jian=0;
<span style="white-space:pre">		</span>while (x)
<span style="white-space:pre">		</span>{
<span style="white-space:pre">			</span>jian+=query(rtL,rtR,1,n,pos[top[x]],pos[x]);
<span style="white-space:pre">			</span>if (son[x]) jian+=query_line(rtL,rtR,1,n,pos[x]+1,low[x])*dis[x];
<span style="white-space:pre">			</span>x=top[x]; y=fa[x][0];
<span style="white-space:pre">			</span>if (y) jian-=query_sz(rtL,rtR,1,n,pos[x],right[x])*dis[y];
<span style="white-space:pre">			</span>x=y;
<span style="white-space:pre">		</span>}
<span style="white-space:pre">		</span>ans=ans-jian*2;
<span style="white-space:pre">		</span>printf("%lld\n",ans);
<span style="white-space:pre">	</span>}
<span style="white-space:pre">	</span>return 0;
}


heheda_is_an_OIer
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
[BZOJ]4012: [HNOI2015]开店 树链剖分+主席树（线段树合并）
TYB的博客
11-23 280
树链剖分+主席树（线段树合并）
BZOJ-4012 开店（动态树分治）
十分残念的博客
08-14 313
传送门：BZOJ-4012 题意：给定一棵树，每个点有一个颜色，多次询问颜色在[l,r]区间内的所有点与某个点之间的距离之和 又是一个数组开小T到怀疑人生系列。。（我怎么总是管不住这双手root数组只开MX呢） 一个经典的动态树分治题目，用树状数组保存每个节点到该层子树的根节点的距离，然后一层一层计算就行。 = =本来应该是这样的，后来看了下发现颜色有n种，n^2logn的树状数组肯定爆内
【BZOJ 4538】[Hnoi2016]网络 树链剖分+堆
pbihao的博客
02-13 240
很巧妙的一道题，因为是树上的操作，所以我们可以联想到树链剖分，但是存在这样的几个问题。 1.每一次查询都是查询不受影响的重要程度 2.每一个操作还可以撤回，当前的最大值在后面的操作中可能会撤回 3.懒惰标记的处理 主要就是这两个问题。 1.为了解决第一个问题，我们可以在每一次修改的时候修改当前需要修改的区间的补集，也就是说线段树维护的就是不经过改点的最大重要度，也就是该点被破坏以后不
bzoj 4012: [HNOI2015]开店 树链剖分+可持久化线段树
beginend
03-20 264
题意 给出一棵数，点有点权。每次询问所有点权在[L,R]之间的点到点u的距离和。强制在线。 n&lt;=150000,Q&lt;=200000 分析 如果没有点权限制的话，我们可以通过树剖线段树，来维护每个点轻儿子子树的节点数量，节点深度和以及节点数量*该点深度的和。 现在有了权值的限制，那么只要把线段树可持久化一下，每次在对应的线段树里面查询即可。 代码 #includ...
bzoj 4012. [HNOI2015]开店（动态点分治）
qq_41997978的博客
02-21 119
先考虑在点分树上如何处理一个弱化的问题：查询其它点到 uuu 点的距离之和。 在点分树上每个节点维护子树到它的距离和，为了方便向上合并其他节点的贡献，再维护每个节点到其父亲节点的距离和，每到一个父亲减一下就得到除这棵子树外其它点到父亲节点的距离和，再加上其到父亲节点的这条边的贡献即可。 维护三个数组：suma[v],sumb[v],siz[v]suma[v],sumb[v],siz[v]sum...
树链剖分_王天懿.ppt
09-14
4. **BZOJ 3589**：动态树的子树权值修改和链上的点权值求和，需要用到树链剖分的链上修改和查询能力。 5. **BZOJ 3626**：求最近公共祖先（LCA）与深度的累加，可以结合LCA算法和树链剖分优化查询。 总之，树链...
BZOJ3230相似子串测试数据
08-19
标题"BZOJ3230相似子串测试数据"指的是一个与在线判题平台BZOJ（BestCoder ZOJ，中国大学算法竞赛在线评测系统）相关的编程问题，具体问题编号为3230，涉及的主要知识点是"相似子串"的搜索和处理。在编程竞赛中，...
BZOJ上面的题目以及数据
10-13
【BZOJ】全称是“北京邮电大学在线评测系统”，它是一个为编程爱好者和参赛者提供在线编程练习和比赛的平台。该平台涵盖了多种算法和编程问题，旨在提升用户的编程技能和解决复杂问题的能力。BZOJ支持C、C++、Java等...
bzoj 4012: [HNOI2015]开店 （树链剖分+主席树）
clover_hxy的博客
04-19 606
题目描述传送门题解这道题维护和求解的方法和bzoj 3924: [Zjoi2015]幻想乡战略游戏是类似的。 但是这道题有一个[L,R]的区间限制，所以我们用主席树来维护，外层是按照离散化后的xix_i从小到大，内层是dfs序。代码#include<iostream> #include<cstring> #include<algorithm> #include<cstdio> #include<c
【BZOJ4012】【HNOI2015】开店（树链剖分，主席树）
დ♂Hany01's Blog Space♂ღ
09-28 240
Description 风见幽香有一个好朋友叫八云紫，她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学。最近她们灵机一动，打算在幻想乡开一家小店来做生意赚点钱。这样的 想法当然非常好啦，但是她们也发现她们面临着一个问题，那就是店开在哪里，面 向什么样的人群。很神奇的是，幻想乡的地图是一个树形结构，幻想乡一共有 n 个地方，编号为 1 到 n，被 n-1 条带权的边连接起来。每个地方都住着...
[HNOI2016]网络 树链剖分，堆
weixin_30445169的博客
02-28 82
[HNOI2016]网络 LG传送门 表示乱搞比正解难想。 整体二分很好想吧。 但是为了好写快乐，我们选择三个\(\log\)的乱搞。 先树剖，线段树套堆维护区间最大值。对于一次修改，如果是插入，就把树上除了这条链的地方加上这个重要度，如果是删除则反之。注意线段树可以标记永久化，这里用的堆是一种（可能）叫懒惰堆的东西，直接看代码都能理解。 //written by newbiechd #inclu...
bzoj4012 [HNOI2015]开店（动态点分治+二分+STL/树链剖分+主席树）
Icefox的博客
02-23 346
如何求所有点到一个点的距离和呢？类似bzoj3924我们的处理方法。 对于每个节点x维护 s1[x]–x的子树中的点到x的距离 s2[x]–x的子树中的点到fa[x]的距离 然后logn的在重心树上跳一跳，减去重复的，补上少的即可。 这题还有一个年龄的限制，怎么办呢？很容易想到线段树，但是我不会写gg 我们还可以利用vector+二分来做。 把所有点按年龄从小到大排序，做前缀和。每次...
BZOJ 4012 HNOI 2015 开店(shop) 一道简单的点剖题
YxuanwKeith的博客
03-25 1493
HNOI 2015 开店(shop) 一道简单的点剖题。 可以用c++自带函数优化代码量
[动态树分治] BZOJ4012 [HNOI2015]开店
雯舞
03-15 2455
膜PoPoQQQ：http://blog.csdn.net/popoqqq/article/details/45365043 省选前学动态树分治  跑起来飞快 45s+ 看着提交记录 只能说惨不忍睹 爆int呆滞了一早上 还有因为没删文件 没强制在线 RE的 动态树分治 印象中还有 幻想乡战略游戏和捉迷藏  代码注释一大片 丑我也不管了 #include #include #inclu
[BZOJ4012][HNOI2015]开店（动态点分治）
女装的你，如此好看！
04-01 443
题目大意就是：一个 nnn 个节点的树，边和点都带权。每次询问点权在 [L,R][L,R][L,R] 范围内的所有点到 uuu 的距离之和，强制在线。 （注意下面所说的「父亲」「子节点」「子树」「祖先」等概念都是在分治树上的） 考虑动态点分治时对于每个点 uuu ，储存 uuu 的子树内所有点到 uuu 的父亲的距离。 但这样显然还不能满足我们的要求，因为我们每次查询的是点权在 [L,R][...
bzoj4012: [HNOI2015]开店
weixin_30797199的博客
12-18 110
传送门：http://www.lydsy.com/JudgeOnline/problem.php?id=4012 思路：首先我们考虑一个简化的问题： 给定一棵树，每次询问所有点到一个点的距离和。 画个图就能知道：距离和=所有点到根的距离和+点数*u到根的距离-每个点与u的lca到根的距离*2 于是问题就成了求lca的dis和 那么我们先对每个点，把...
BZOJ 4012 HNOI 2015 开店 动态点分治
とある黄鱼のOIブログ
04-03 1038
多次询问给定点到点权在[l,r][l,r]间的所有点距离之和。询问某个点x到其他点的距离之和，可以考虑其到根的路径上的点，发现其和可以表示为（为了简单表示，令c[x]表示x所在子树内的所有点） c[x]到x的距离和+ {c[fa[x]]-c[x]}到fa[x]的距离和+|{c[fa[x]]-c[x]}|*dis(x,fa[x])+ …. c[fa[x]]-c[x]表示x的父亲的子树的所有点除
bzoj 1040 zjoi2008 骑士 树形dp
最新发布
08-08
我们可以使用树形dp来解决这个问题，具体来说，我们可以将这个图看作是一棵树，每个点的父节点是它的前驱或者后继，然后我们从根节点开始，依次向下遍历，对于每个节点，我们可以考虑它的两个儿子，如果它的两个儿子...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值