区间修改的kd-tree,打标记,下传。

每次询问的时候,从询问点向上找到根,然后依次下传下来,再回答询问。

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define N 100001
#define KD 2
int n,root,m,q,qp[2][KD],fa[N],val,delta[N];
bool dn;
struct Node
{
int ch[2],w,minn[KD],maxx[KD],p[KD],id;
void Init()
{
for(int i=0;i<KD;++i)
minn[i]=maxx[i]=p[i];
}
}T[N];
bool operator < (const Node &a,const Node &b){return a.p[dn] < b.p[dn];}
inline void pushup(const int &rt)
{
for(int i=0;i<2;++i)
if(T[rt].ch[i])
for(int j=0;j<KD;++j)
{
T[rt].minn[j]=min(T[rt].minn[j],T[T[rt].ch[i]].minn[j]);
T[rt].maxx[j]=max(T[rt].maxx[j],T[T[rt].ch[i]].maxx[j]);
}
}
inline void pushdown(const int &rt)
{
if(delta[rt])
{
T[rt].w=delta[rt];
for(int i=0;i<2;++i)
delta[T[rt].ch[i]]=delta[rt];
delta[rt]=0;
}
}
int buildtree(int l=1,int r=n,bool d=0)
{
dn=d;
int m=(l+r>>1);
nth_element(T+l,T+m,T+r+1);
T[m].Init();
if(l!=m) T[m].ch[0]=buildtree(l,m-1,d^1);
if(m!=r) T[m].ch[1]=buildtree(m+1,r,d^1);
pushup(m);
return m;
}
void Update(int rt=root)
{
if(qp[0][0] <= T[rt].minn[0] && T[rt].maxx[0] <= qp[0][1]
&& qp[1][0] <= T[rt].minn[1] && T[rt].maxx[1] <= qp[1][1])
{
delta[rt]=val;
return;
}
pushdown(rt);
if(qp[0][0] <= T[rt].p[0] && T[rt].p[0] <= qp[0][1]
&& qp[1][0] <= T[rt].p[1] && T[rt].p[1] <= qp[1][1])
T[rt].w=val;
for(int i=0;i<2;++i)
if(T[rt].ch[i]
&& qp[0][0] <= T[T[rt].ch[i]].maxx[0] && T[T[rt].ch[i]].minn[0] <= qp[0][1]
&& qp[1][0] <= T[T[rt].ch[i]].maxx[1] && T[T[rt].ch[i]].minn[1] <= qp[1][1])
Update(T[rt].ch[i]);
}
void Query(int U)
{
if(fa[U])
Query(fa[U]);
pushdown(U);
}
int zu;
int v[N],next[N],first[N],e;
void AddEdge(const int &U,const int &V)
{
v[++e]=V;
next[e]=first[U];
first[U]=e;
}
int dep[N],dfn[N],dfr[N];
void dfs(int U)
{
dfn[U]=++e;
T[U].w=1;
T[U].p[0]=e;
T[U].p[1]=dep[U];
T[U].id=U;
for(int i=first[U];i;i=next[i])
{
dep[v[i]]=dep[U]+1;
dfs(v[i]);
}
dfr[U]=e;
}
typedef long long ll;
#define MOD 1000000007ll
ll ans;
int ma[N];
int main()
{
// freopen("bzoj4154.in","r",stdin);
scanf("%d",&zu);
for(;zu;--zu)
{
int x,dis;
scanf("%d%d%d",&n,&m,&q);
for(int i=2;i<=n;++i)
{
scanf("%d",&x);
AddEdge(x,i);
}
e=0; dep[1]=1;
dfs(1);
buildtree();
root=(1+n>>1);
for(int i=1;i<=n;++i)
{
ma[T[i].id]=i;
for(int j=0;j<2;++j)
if(T[i].ch[j])
fa[T[i].ch[j]]=i;
}
for(int i=1;i<=q;++i)
{
scanf("%d%d%d",&x,&dis,&val);
if(!val)
{
Query(ma[x]);
ans=(ans+(ll)i*(ll)T[ma[x]].w%MOD)%MOD;
}
else
{
qp[0][0]=dfn[x];
qp[0][1]=dfr[x];
qp[1][0]=dep[x];
qp[1][1]=dep[x]+dis;
Update();
}
}
printf("%lld\n",ans);
for(int i=1;i<=n;++i)
T[i].ch[0]=T[i].ch[1]=0;
memset(delta+1,0,sizeof(int)*n);
memset(fa+1,0,sizeof(int)*n);
ans=e=0;
memset(first+1,0,sizeof(int)*n);
}
return 0;
}

【kd-tree】bzoj4154 [Ipsc2015]Generating Synergy的更多相关文章

  1. BZOJ4154:[Ipsc2015]Generating Synergy(K-D Tree)

    Description 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问点a的颜色 Input 第一行一个数T,表示数据组数 接下来每组数据的第一行三 ...

  2. BZOJ4154: [Ipsc2015]Generating Synergy

    Description 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问点a的颜色   Input 第一行一个数T,表示数据组数 接下来每组数据的第一 ...

  3. 【BZOJ4154】[Ipsc2015]Generating Synergy KDtree

    [BZOJ4154][Ipsc2015]Generating Synergy Description 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问 ...

  4. BZOJ4154:[IPSC2015]Generating Synergy

    浅谈\(K-D\) \(Tree\):https://www.cnblogs.com/AKMer/p/10387266.html 题目传送门:https://lydsy.com/JudgeOnline ...

  5. [bzoj4154][Ipsc2015]Generating Synergy_KD-Tree_dfs序

    Generating Synergy bzoj-4154 Ipsc-2015 题目大意:给定一棵n个节点树,m个操作,支持:将一个点周围所有距该点距离不超过l的子结点的颜色改成另一种颜色:查询单点颜色 ...

  6. 【BZOJ4154】Generating Synergy【kd树】

    题意 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问点a的颜色 分析 我们以dfs序为横坐标,深度为纵坐标,建kd树.我们每次更新,都是在kd树中更 ...

  7. 【bzoj4154】[Ipsc2015]Generating Synergy KD-tree

    题目描述 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问点a的颜色 输入 第一行一个数T,表示数据组数 接下来每组数据的第一行三个数n,c,q表示结 ...

  8. 【bzoj 4154】[Ipsc2015]Generating Synergy

    题目 大概已经掌握熟练码出\(kdt\)的技能了 发现距离子树根节点\(x\)不超过\(l\)的点可以用两种方式来限制,首先\(dfs\)序在\([dfn_x,dfn_x+sum_x)\)中,深度自然 ...

  9. 【Symmetric Tree】cpp

    题目: Given a binary tree, check whether it is a mirror of itself (ie, symmetric around its center). F ...

随机推荐

  1. No.1__C#

    这是第一篇C#的日记,到现在为止已经学习了一个礼拜的C#了.由于是实习中才开始学习,所以这次不准备像在大学学习那样,拿着课本划重点,背概念.这应当是一门实践的课程,应该一边编程,一边学.这是到公司第一 ...

  2. Flex 加载shp

    至于gis格式比较常见的shp是开源的,网上开源的as代码也很多 这个支持的shp算比较好的 源码在这边http://files.cnblogs.com/files/haibalai/shp.rar, ...

  3. google快捷键,通过浏览器本身来查看

    今天玩google浏览器时发现一个不需要访问google就能查看其快捷键的方式,再此记录一下,以备后用. 1:打开Google浏览器 2:打开开发者工具,有三种方法 2-1:按F12快捷键 2-2:按 ...

  4. window.parent

    今天出现了一个严重的bug,页面在载入时居然导致了页面崩溃,找了近5个小时的代码,才找到了原因(擦,我只能说我水平有限,浪费了太长时间了~). 源代码: var curr = window; whil ...

  5. 使用celery的backend异步获取结果

    惯例先贴出相关参考的文档: http://docs.celeryproject.org/en/stable/getting-started/next-steps.html http://docs.ce ...

  6. SQL笔记 [长期更新] (-2015.4)

    [遍历所有表,复制表结构,复制表数据] --插入语句SELECT * INTO A FROM B 是在还没有A表的情况下,直接通过B表创建并把B表数据复制到A表里面,之后A,B表的结构和数据完全一样. ...

  7. SQL递归

    递归一般出现在树形结构中 1:根据孩子节点查找所有父节点 With T As ( Select * From U_Companies TB Where CompanyID=80047 Union Al ...

  8. 使用jna调用dll,jdk位数和dll位数的关系

    最近在学习jna,发现dll文件能能否成功调用取决于jdk位数. 32位jdk只能使用32位的dll,64位jdk只能使用64位的dll,否则位数不对应的话报的错是 "Exception i ...

  9. srum.3

    1. 确保product backlog井然有序.(参考示例图1) 2. Sprint周期,一个冲刺周期,长度定为两周,本学期还有三个冲刺周期. 3. 确定Sprint目标. 3.1 产品负责人概括产 ...

  10. Mybatis按SQL查询字段的顺序返回查询结果

    在SpringMVC+Mybatis的开发过程中,可以通过指定resultType="hashmap"来获得查询结果,但其输出是没有顺序的.如果要按照SQL查询字段的顺序返回查询结 ...