BZOJ4154:[IPSC2015]Generating Synergy
浅谈\(K-D\) \(Tree\):https://www.cnblogs.com/AKMer/p/10387266.html
题目传送门:https://lydsy.com/JudgeOnline/problem.php?id=4154
每个节点看做是点\((dfn_i,dep_i)\),然后对于距离\(a\)的不超过\(l\)的在\(a\)子树中的点就是正方形\([dfn_a,dfn_a+siz_a-1][dep_a,dep_a+l]\),直接范围\(cover\)和单点查询即可。
时间复杂度:\(O(n\sqrt{n})\)
空间复杂度:\(O(n)\)
代码如下:
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn=1e5+5,inf=2e9,mo=1e9+7;
int n,m,c,tot,cnt,pps,ans,Case;
int now[maxn],pre[maxn<<1],son[maxn<<1];
int dfn[maxn],node[maxn],siz[maxn],dep[maxn];
int read() {
int x=0,f=1;char ch=getchar();
for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+ch-'0';
return x*f;
}
void add(int a,int b) {
pre[++tot]=now[a];
now[a]=tot,son[tot]=b;
}
void dfs(int u,int D) {
dfn[u]=++cnt,dep[u]=D,siz[u]=1;
for(int p=now[u],v=son[p];p;p=pre[p],v=son[p])
dfs(v,D+1),siz[u]+=siz[v];
}
struct kd_tree {
int root,top;
int stk[maxn],fa[maxn];
struct point {
int cov,col,ls,rs;
int c[2],mn[2],mx[2];
point() {}
point(int _x,int _y) {
mn[0]=mx[0]=c[0]=_x;
mn[1]=mx[1]=c[1]=_y;
col=1,cov=-1,ls=rs=0;
}
bool operator<(const point &a)const {
return c[pps]<a.c[pps];
}
}p[maxn];
int build(int l,int r,int d) {
int mid=(l+r)>>1,u=mid;pps=d;
nth_element(p+l,p+mid,p+r+1);
node[p[u].c[0]]=u;
if(l<mid)p[u].ls=build(l,mid-1,d^1);
if(r>mid)p[u].rs=build(mid+1,r,d^1);
int ls=p[u].ls,rs=p[u].rs;fa[ls]=fa[rs]=u;
for(int i=0;i<2;i++) {
int mn=min(p[ls].mn[i],p[rs].mn[i]);
p[u].mn[i]=min(p[u].mn[i],mn);
int mx=max(p[ls].mx[i],p[rs].mx[i]);
p[u].mx[i]=max(p[u].mx[i],mx);
}
return u;
}
void prepare() {
p[0].mn[0]=p[0].mn[1]=inf;
p[0].mx[0]=p[0].mx[1]=-inf;
for(int i=1;i<=n;i++)
p[i]=point(dfn[i],dep[i]);
root=build(1,n,0);
}
void make_tag(int u,int c) {
p[u].col=p[u].cov=c;
}
void push_down(int u) {
if(p[u].cov==-1)return;
if(p[u].ls)make_tag(p[u].ls,p[u].cov);
if(p[u].rs)make_tag(p[u].rs,p[u].cov);
p[u].cov=-1;
}
void query(int u,int id) {
int tmp=u;
while(fa[u])stk[++top]=fa[u],u=fa[u];
while(top)push_down(stk[top--]);
ans=(ans+1ll*id*p[tmp].col%mo)%mo;
}
void change(int u,int x1,int x2,int y1,int y2,int color) {
if(x2<p[u].mn[0]||x1>p[u].mx[0])return;
if(y2<p[u].mn[1]||y1>p[u].mx[1])return;
bool bo1=(x1<=p[u].mn[0]&&p[u].mx[0]<=x2);
bool bo2=(y1<=p[u].mn[1]&&p[u].mx[1]<=y2);
if(bo1&&bo2) {make_tag(u,color);return;}
push_down(u);
bo1=(x1<=p[u].c[0]&&p[u].c[0]<=x2);
bo2=(y1<=p[u].c[1]&&p[u].c[1]<=y2);
if(bo1&&bo2)p[u].col=color;
if(p[u].ls)change(p[u].ls,x1,x2,y1,y2,color);
if(p[u].rs)change(p[u].rs,x1,x2,y1,y2,color);
}
}T;
void clear() {
ans=tot=cnt=pps=0;
memset(now,0,sizeof(now));
memset(T.fa,0,sizeof(T.fa));
}
int main() {
Case=read();
while(Case--) {
n=read(),c=read(),m=read(),clear();
for(int i=2,x;i<=n;i++)
x=read(),add(x,i);
dfs(1,1);T.prepare();
for(int i=1;i<=m;i++) {
int a=read(),l=read(),opt=read();
if(!opt)T.query(node[dfn[a]],i);
else T.change(T.root,dfn[a],dfn[a]+siz[a]-1,dep[a],dep[a]+l,opt);
}
printf("%d\n",ans);
}
return 0;
}
BZOJ4154:[IPSC2015]Generating Synergy的更多相关文章
- 【BZOJ4154】[Ipsc2015]Generating Synergy KDtree
[BZOJ4154][Ipsc2015]Generating Synergy Description 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问 ...
- 【bzoj4154】[Ipsc2015]Generating Synergy KD-tree
题目描述 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问点a的颜色 输入 第一行一个数T,表示数据组数 接下来每组数据的第一行三个数n,c,q表示结 ...
- 【kd-tree】bzoj4154 [Ipsc2015]Generating Synergy
区间修改的kd-tree,打标记,下传. 每次询问的时候,从询问点向上找到根,然后依次下传下来,再回答询问. #include<cstdio> #include<algorithm& ...
- BZOJ4154: [Ipsc2015]Generating Synergy
Description 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问点a的颜色 Input 第一行一个数T,表示数据组数 接下来每组数据的第一 ...
- BZOJ4154:[Ipsc2015]Generating Synergy(K-D Tree)
Description 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问点a的颜色 Input 第一行一个数T,表示数据组数 接下来每组数据的第一行三 ...
- 【bzoj 4154】[Ipsc2015]Generating Synergy
题目 大概已经掌握熟练码出\(kdt\)的技能了 发现距离子树根节点\(x\)不超过\(l\)的点可以用两种方式来限制,首先\(dfs\)序在\([dfn_x,dfn_x+sum_x)\)中,深度自然 ...
- BZOJ 4154: [Ipsc2015]Generating Synergy KDtree+dfs序
多组数据真tm恶心~ 把 $dfs$序和深度分别看作横纵坐标,然后用 $KDtree$ 数点就可以了~ #include <cstdio> #include <cstring> ...
- 4154: [Ipsc2015]Generating Synergy
Description 给定一棵以1为根的有根树,初始所有节点颜色为1,每次将距离节点a不超过l的a的子节点染成c,或询问点a的颜色 区间修改单点查询kdtree #include<iostre ...
- [bzoj4154][Ipsc2015]Generating Synergy_KD-Tree_dfs序
Generating Synergy bzoj-4154 Ipsc-2015 题目大意:给定一棵n个节点树,m个操作,支持:将一个点周围所有距该点距离不超过l的子结点的颜色改成另一种颜色:查询单点颜色 ...
随机推荐
- 添加code到github上
第一步:github上新建远程仓库 1. 在 https://github.com/ 注册账号 2. new 一个新仓库 (1) 点击加号下的`New repository` (2)在Reposit ...
- 初识python---简介,简单的for,while&if
一编程语言:编程语言是程序员与计算机沟通的介质: 编程语言的分类: 1机器语言:是用二进制代码表示的计算机能直接识别和执行的一种机器指令的集合. 优点:灵活,直接执行和速度快 ...
- FANCO工程机械云平台
此系统专门为工程机械设备使用单位定制合理的生产要求,监控生产状况,快速体现产值,通过算法计算绩效等
- 快乐学习 Ionic Framework+PhoneGap 手册1-1{创建APP项目}
快乐学习 Ionic Framework+PhoneGap 手册1-1 * 前提必须安装 Node.js,安装PhoneGap,搭建Android开发环境,建议使用真机调试 {1.1}= 创建APP项 ...
- jQuery滑动杆打分插件
在线演示 本地下载
- ORA-01034和ORA-27101的错误
我本机安装的数据库版本是ORACLE 11G R2,用plsql连接时候,报ora-12514如下错误: 但是在cmd里用sqlplus连接已经创建的用户时候,报如下错误: ORA-01034 - O ...
- BufferingForwardingAppender in log4net
https://blog.csdn.net/szx1999/article/details/50073857 7. 写日志会影响系统性能吗? 写日志必然是会消耗一定资源的,而RollingFileAp ...
- eclipse添加删除插件-eclipse marketplace
源文地址:http://jingyan.baidu.com/article/cdddd41c5c883353cb00e19e.html 在有些版本的eclips上并没有eclipse marketpl ...
- HDU 4348 To the moon (主席树区间更新)
题意:首先给你n个数,开始时间为0,最后按照操作输出 给你四种操作: 1. C l r d : 在(l,r)区间都加上d,时间加一2. Q l r : 询问现在(l,r)的区间和3. H l r ...
- http 常见的错误码
翻译自:https://www.w3.org/Protocols/rfc2616/rfc2616-sec10.html 常见错误码 一.信息 1XX (Information 1xx) ——这一类的状 ...