【BZOJ4771】七彩树(主席树)

题面

BZOJ

题解

如果没有深度限制,每次只询问子树内的颜色个数,除了树套树\(dfs\)序加前驱或者后继强行二维数点之外,还有这样一种做法:

把所有相同颜色的点按照\(dfs\)序排序,每个点给自己的位置贡献\(1\),相邻的两个点给\(lca\)贡献\(-1\)。然后只要区间内存在这种颜色,则其子树内的权值和必定为\(1\)。那么只需要这样子染好所有颜色之后询问子树和。

然而这题要求的是深度在一个范围内的东西。

如果可以离线,我们可以把所有点按照深度排序,从上往下依次加入,计算贡献,然后把询问在\(dep[x]+d\)的地方计算答案。

现在不能离线,就提前把这个东西的线段树给算好,拿主席树直接提前维护好就行了。

时间复杂度是一个\(\log\)。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<set>
using namespace std;
#define MAX 100100
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
struct Node{int ls,rs,v;}t[MAX<<7];
int rt[MAX],tot;
void Modify(int &x,int l,int r,int p,int w)
{
t[++tot]=t[x];x=tot;t[x].v+=w;if(l==r)return;
int mid=(l+r)>>1;
if(p<=mid)Modify(t[x].ls,l,mid,p,w);
else Modify(t[x].rs,mid+1,r,p,w);
}
int Query(int A,int l,int r,int L,int R)
{
if(L<=l&&r<=R)return t[A].v;
int mid=(l+r)>>1,ret=0;
if(L<=mid)ret+=Query(t[A].ls,l,mid,L,R);
if(R>mid)ret+=Query(t[A].rs,mid+1,r,L,R);
return ret;
}
struct Line{int v,next;}e[MAX];
int h[MAX],cnt=1;
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
int dfn[MAX],tim,dep[MAX],top[MAX],fa[MAX],size[MAX],hson[MAX],ln[MAX],low[MAX];
int n,Q,c[MAX],p[MAX];
set<int> S[MAX];set<int>::iterator t1,t2,t3;
void init()
{
memset(rt,0,sizeof(rt));tot=0;
memset(h,0,sizeof(h));cnt=1;
memset(dfn,0,sizeof(dfn));tim=0;
}
void dfs1(int u,int ff)
{
dep[u]=dep[ff]+1;size[u]=1;hson[u]=0;fa[u]=ff;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(v==ff)continue;
dfs1(e[i].v,u);size[u]+=size[v];
if(size[v]>size[hson[u]])hson[u]=v;
}
}
void dfs2(int u,int tp)
{
top[u]=tp;dfn[u]=++tim;ln[tim]=u;
if(hson[u])dfs2(hson[u],tp);
for(int i=h[u];i;i=e[i].next)
if(e[i].v!=fa[u]&&e[i].v!=hson[u])
dfs2(e[i].v,e[i].v);
low[u]=tim;
}
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;
}
bool cmp(int a,int b){return dep[a]<dep[b];}
int St[MAX],sum;
int main()
{
int T=read();
while(T--)
{
init();n=read();Q=read();
for(int i=1;i<=n;++i)c[i]=read(),St[i]=c[i];
sort(&St[1],&St[n+1]);sum=unique(&St[1],&St[n+1])-St-1;
for(int i=1;i<=n;++i)c[i]=lower_bound(&St[1],&St[sum+1],c[i])-St;
for(int i=2;i<=n;++i)Add(read(),i);
dfs1(1,0);dfs2(1,1);
for(int i=1;i<=n;++i)p[i]=i;
sort(&p[1],&p[n+1],cmp);
for(int i=1;i<=sum;++i)S[i].clear();
for(int i=1,j=1;i<=dep[p[n]];++i)
{
rt[i]=rt[i-1];
while(j<=n&&dep[p[j]]==i)
{
int u=p[j];
S[c[u]].insert(dfn[u]);
t1=t2=t3=S[c[u]].find(dfn[u]);
t2--;t3++;Modify(rt[i],1,n,dfn[u],1);
if(t1!=S[c[u]].begin())Modify(rt[i],1,n,dfn[LCA(ln[*t2],u)],-1);
if(t3!=S[c[u]].end())Modify(rt[i],1,n,dfn[LCA(ln[*t3],u)],-1);
if(t1!=S[c[u]].begin()&&t3!=S[c[u]].end())Modify(rt[i],1,n,dfn[LCA(ln[*t2],ln[*t3])],1);
++j;
}
}
int ans=0;
while(Q--)
{
int x=read()^ans,d=read()^ans;
printf("%d\n",ans=Query(rt[min(dep[x]+d,dep[p[n]])],1,n,dfn[x],low[x]));
}
}
}

【BZOJ4771】七彩树(主席树)的更多相关文章

  1. 线段树简单入门 (含普通线段树, zkw线段树, 主席树)

    线段树简单入门 递归版线段树 线段树的定义 线段树, 顾名思义, 就是每个节点表示一个区间. 线段树通常维护一些区间的值, 例如区间和. 比如, 上图 \([2, 5]\) 区间的和, 为以下区间的和 ...

  2. 【BZOJ4771】七彩树 主席树+树链的并

    [BZOJ4771]七彩树 Description 给定一棵n个点的有根树,编号依次为1到n,其中1号点是根节点.每个节点都被染上了某一种颜色,其中第i个节点的颜色为c[i].如果c[i]=c[j], ...

  3. [BZOJ4771]七彩树(主席树)

    https://blog.csdn.net/KsCla/article/details/78249148 用类似经典的链上区间颜色计数问题的做法,这个题可以看成是询问DFS在[L[x],R[x]]中, ...

  4. 【BZOJ3439】Kpm的MC密码 trie树+主席树

    Description 背景 想Kpm当年为了防止别人随便进入他的MC,给他的PC设了各种奇怪的密码和验证问题(不要问我他是怎么设的...),于是乎,他现在理所当然地忘记了密码,只能来解答那些神奇的身 ...

  5. 学习笔记--函数式线段树(主席树)(动态维护第K极值(树状数组套主席树))

    函数式线段树..资瓷 区间第K极值查询 似乎不过似乎划分树的效率更优于它,但是如果主席树套树状数组后,可以处理动态的第K极值.即资瓷插入删除,划分树则不同- 那么原理也比较易懂: 建造一棵线段树(权值 ...

  6. HDU5790 Prefix 字典树+主席树

    分析:这个题和spoj的d_query是一个题,那个是求一段区间里有多少个不同的数字,这里是统计有多少个不同的前缀 用字典树进行判重,(和查询不同的数字一样)对于每个不同的前缀,只保留它最后一次出现的 ...

  7. bzoj 3545&&3551: [ONTAK2010]Peaks &&加强版 平衡树&&并查集合并树&&主席树

    3545: [ONTAK2010]Peaks Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 635  Solved: 177[Submit][Stat ...

  8. BZOJ 4539: [Hnoi2016]树 [主席树 lca]

    4539: [Hnoi2016]树 题意:不想写.复制模板树的子树,查询两点间距离. *** 终于有一道会做的题了...... 画一画发现可以把每次复制的子树看成一个大点来建一棵树,两点的lca一定在 ...

  9. 线段树(单标记+离散化+扫描线+双标记)+zkw线段树+权值线段树+主席树及一些例题

    “队列进出图上的方向 线段树区间修改求出总量 可持久留下的迹象 我们 俯身欣赏” ----<膜你抄>     线段树很早就会写了,但一直没有总结,所以偶尔重写又会懵逼,所以还是要总结一下. ...

  10. UOJ#218. 【UNR #1】火车管理 线段树 主席树

    原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ218.html 题解 如果我们可以知道每次弹出栈之后新的栈顶是什么,那么我们就可以在一棵区间覆盖.区间求和 ...

随机推荐

  1. charles如何设置弱网

  2. (Git 学习)Git SSH Key 创建步骤

    首先感谢segmentfalut上的朋友对我帮助. 首先:查看你是否有../ssh 这个文件:怎么查看:找到你的git安装目录,在安装目录下查看是否./ssh,以我的为例: 在C盘/Users/11/ ...

  3. Es5中的类和静态方法 继承

    Es5中的类和静态方法 继承(原型链继承.对象冒充继承.原型链+对象冒充组合继承) // es5里面的类 //1.最简单的类 // function Person(){ // this.name='张 ...

  4. Migrate MySQL database using dump and restore

    kaorimatz/mysqldump-loader: Load a MySQL dump file using LOAD DATA INFILEhttps://github.com/kaorimat ...

  5. 父级div宽度100%,子级一个div宽度固定,另一个宽度自适应

    <!DOCTYPE html> <html> <head> <title>布局测试</title> <style type=" ...

  6. python3 常见的两种文件上传方法

    1.上传页面带input type格式send_keys传值方式上传不能大于60k(具体看开发设置的value)文件大小 fx.find_element_by_id('xx').send_keys(r ...

  7. JMeter压测分布式部署

    监控JMeter压力机的性能

  8. [转帖]Linux的标准输入 标准输出和错误输出

    Linux标准输入.输出和错误和文件重定向 专题 https://www.cnblogs.com/softidea/p/3965093.html 感觉自己对 这一块的理解一直不好 昨天同事给了一个 b ...

  9. Mysql 5.7 Windows 版本(zip)的安装简单过程

    1. 下载zip包 https://cdn.mysql.com//Downloads/MySQL-5.7/mysql-5.7.25-winx64.zip 2. 找一个目录解压缩 3. 简单进行安装: ...

  10. .Net中EF通用数据层小结

    增删改查: using System; using System.Collections.Generic; using System.Data; using System.Data.Entity; u ...