洛谷P3703 [SDOI2017]树点涂色(LCT,dfn序,线段树,倍增LCA)
闲话
这是所有LCT题目中的一个异类。
之所以认为是LCT题目,是因为本题思路的瓶颈就在于如何去维护同颜色的点的集合。
只不过做着做着,感觉后来的思路(dfn序,线段树,LCA)似乎要喧宾夺主了。。。(至少在代码上看是如此)
思路分析
一个一个操作来(瞎BB中,这种思路模式并不具有普遍性。。。。。。)
1操作
还好我没学树剖233333以至于(直接想到)只好用LCT来维护颜色。
题目透露出的神奇的性质——每一种颜色,无论在任何时刻,肯定是一条链,而且点的深度严格递增!
而且还特意指定根节点!1操作特意修改x到根节点的颜色!
想到了这里,就不难想到本题的关键模型——LCT中每个Splay辅助树维护同颜色点的集合
于是1操作==access。。。。。。
2操作
x到y路径?split?!
I'm too young too simple
LCT维护了集合,就只能维护集合了。随便再乱搞一下集合就被破坏了。
于是就要再外部维护了。
至于维护什么,现在其实还不能产生很好的思路。。。。。
但我们可以先想到一点:没有了LCT,还要资磁树中任意两点之间的询问?
常见的复杂度正确的方法(套路)能想到的就只有树上差分了吧(设F为状态,那么就形如F[x]+F[y]-F[lca])
于是就可以只维护每个点到根节点路径上的颜色种数
转化一下,在LCT中就等价于每个点到根节点路径所要经过的轻边总数,这里又是一个关键点
查询F[x]+F[y]-2F[lca]+1(+1是因为lca所在的颜色被减了两次)
不会树剖,只好写倍增LCA
然后就可以转而思考如何维护这个状态了
首先初始状态就是每个点的深度,然后就接着考虑修改了。
因为状态只与轻边有关,所以在access时更改就可以啦
可以类比一下LCT维护子树信息和(可参考一下Blog的LCT总结)
access中有替换右儿子的操作,等于把原来一条边变轻,新的一条边变重
那么原来那条边所指的子树状态全部要+1(多了一个轻边),新连上的边所指的子树状态全部要-1
于是问题又出现了。。。
众所周知,LCT可以维护子树信息,但不可以修改子树信息
树剖很好维护就不提了,然后我又不会树剖TOT(我太弱了)
在这紧要关头,dfn序救了我。。。。。。
一个子树,所有的点的dfn序一定是连续的区间
所以维护线段树,表示dfn序的区间,修改的时候在对应的区间修改,查询就单点查
3操作
所有的思路难点,在操作2冗长的思路分析中都攻克了
这里就在线段树里维护状态最大值(树剖也一样),区间查询就好啦
至于写法,线段树里的区间加减法可以写懒标记,也可以实现永久化标记(YL巨佬做法,常数暴踩本蒟蒻,目前rank1)
只不过我试了一下,dfn序线段树写永久化标记因为某些无法描述的玄学问题变得更慢了。。。。。。
思路就这样,有些细节在代码里(Debug一晚上带来的惨痛的经验。。。。。。)
算上in,pup,pdn和main,此程序一共有15个函数。。。。。。
#include<cstdio>
#include<iostream>
using namespace std;
#define I inline
#define R register int
#define G ch=getchar()
#define in(z) G;\
while(ch<'-')G;\
z=ch&15;G;\
while(ch>'-')z*=10,z+=ch&15,G
#define lc x<<1
#define rc x<<1|1
#define pup mx[x]=max(mx[lc],mx[rc])
#define pdn if(lz[x])upd(lc,l[x],m[x],lz[x]),upd(rc,m[x]+1,r[x],lz[x]),lz[x]=0
//都是线段树操作,本题的LCT内部没有维护信息
const int N=100009,M=N*20;
int l[M],m[M],r[M],mx[M],lz[M],f[N],c[N][2],st[N][20],o[N];
int p=1,he[N],ne[N<<1],to[N<<1],d[N],dfn[N],at[N],mr[N],now;
//at是dfn的反表示,mr表示每个点的子树在dfn序区间中的右端点(左端点是它自己)
I void dfs(R x,R fa){//建树预处理
d[now=at[dfn[x]=++p]=x]=d[st[x][0]=f[x]=fa]+1;//一堆信息的预处理压进了一行
for(R&i=o[x];(st[x][i+1]=st[st[x][i]][i]);++i);//倍增LCA预处理
for(R i=he[x];i;i=ne[i])
if(fa!=to[i])dfs(to[i],x);
mr[x]=now;
}
I int lca(R x,R y){//求LCA
if(d[x]<d[y])swap(x,y);
for(R i=o[x];i>=0;--i)
if(d[st[x][i]]>=d[y])x=st[x][i];//Debug中的错误1:>=写成了>
if(x==y)return x;
for(R i=o[x];i>=0;--i)
if(st[x][i]!=st[y][i])x=st[x][i],y=st[y][i];
return st[x][0];
}
I void build(R x,R s,R e){//建线段树
l[x]=s;r[x]=e;m[x]=(s+e)>>1;
if(s==e){mx[x]=d[at[s]];return;}//利用反表示找到初始状态
build(lc,s,m[x]);build(rc,m[x]+1,e);pup;
}
I void upd(R x,R s,R e,R v){//区间修改
if(l[x]==s&&r[x]==e){mx[x]+=v;lz[x]+=v;return;}//注意mx也要变
pdn;
if(e<=m[x])upd(lc,s,e,v);
else if(s>m[x])upd(rc,s,e,v);
else upd(lc,s,m[x],v),upd(rc,m[x]+1,e,v);
pup;
}
I int get(R s){//单点查值,与区间查值分开了,为了减小常数
R x=1;
while(l[x]!=r[x]){
pdn;x=(lc)+(s>m[x]);
}
return mx[x];
}
I int ask(R x,R s,R e){//区间查值
if(l[x]==s&&r[x]==e)return mx[x];
pdn;
if(e<=m[x])return ask(lc,s,e);
if(s>m[x])return ask(rc,s,e);
return max(ask(lc,s,m[x]),ask(rc,m[x]+1,e));
}
I bool nrt(R x){//LCT部分
return c[f[x]][0]==x||c[f[x]][1]==x;
}
I void rot(R x){
R y=f[x],z=f[y],k=c[y][1]==x,w=c[x][!k];
if(nrt(y))c[z][c[z][1]==y]=x;c[x][!k]=y;c[y][k]=w;
f[w]=y;f[y]=x;f[x]=z;//Debug中的错误2:y写成了x
}
I void splay(R x){
R y;
while(nrt(x)){
if(nrt(y=f[x]))rot((c[f[y]][0]==y)^(c[y][0]==x)?x:y);
rot(x);
}
}
I int frt(R x){//有别于传统意义下的findroot
while(c[x][0])x=c[x][0];
return x;
}
I void access(R x){
for(R w,y=0;x;x=f[y=x]){
splay(x);
if(c[x][1])w=frt(c[x][1]),upd(1,dfn[w],dfn[mr[w]],1);
if((c[x][1]=y))w=frt(y),upd(1,dfn[w],dfn[mr[w]],-1);
//Debug中的错误3:这里更新要找原子树的根(即深度最小的那个点)
//而不能把辅助树的根当原子树根直接upd(1,dfn[c[x][1]],dfn[mr[c[x][1]]],1)
}
}
int main(){
register char ch;
R n,m,i,a,b,op,x,y;
in(n);in(m);
for(i=1;i<n;++i){
in(a);in(b);
to[++p]=b;ne[p]=he[a];he[a]=p;
to[++p]=a;ne[p]=he[b];he[b]=p;
}
p=0;dfs(1,0);build(1,1,n);
while(m--){
in(op);in(x);
if(op==1)access(x);
else if(op==2){
in(y);
printf("%d\n",get(dfn[x])+get(dfn[y])-get(dfn[lca(x,y)])*2+1);
}
else printf("%d\n",ask(1,dfn[x],dfn[mr[x]]));
}
return 0;
}
洛谷P3703 [SDOI2017]树点涂色(LCT,dfn序,线段树,倍增LCA)的更多相关文章
- [BZOJ4817][SDOI2017]树点涂色(LCT+DFS序线段树)
4817: [Sdoi2017]树点涂色 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 692 Solved: 408[Submit][Status ...
- BZOJ.4817.[SDOI2017]树点涂色(LCT DFS序 线段树)
题目链接 操作\(1.2\)裸树剖,但是操作\(3\)每个点的答案\(val\)很不好维护.. 如果我们把同种颜色的点划分到同一连通块中,那么向根染色的过程就是Access()! 最初所有点间都是虚边 ...
- bzoj4817/luogu3703 树点涂色 (LCT+dfs序+线段树)
我们发现,这个染色的操作他就很像LCT中access的操作(为什么??),然后就自然而然地想到,其实一个某条路径上的颜色数量,就是我们做一个只有access操作的LCT,这条路径经过的splay的数量 ...
- P3703 [SDOI2017]树点涂色 LCT维护颜色+线段树维护dfs序+倍增LCA
\(\color{#0066ff}{ 题目描述 }\) Bob有一棵\(n\)个点的有根树,其中1号点是根节点.Bob在每个点上涂了颜色,并且每个点上的颜色不同. 定义一条路径的权值是:这条路径上的点 ...
- bzoj 3779 重组病毒 好题 LCT+dfn序+线段树分类讨论
题目大意 1.将x到当前根路径上的所有点染成一种新的颜色: 2.将x到当前根路径上的所有点染成一种新的颜色,并且把这个点设为新的根: 3.查询以x为根的子树中所有点权值的平均值. 分析 原题codec ...
- [BZOJ4817][SDOI2017]树点涂色:Link-Cut Tree+线段树
分析 与[BZOJ3779]重组病毒唯一的区别是多了一个链上求实链段数的操作. 因为每条实链的颜色必然不相同且一条实链上不会有两个深度相同的点(好像算法的正确性和第二个条件没什么关系,算了算了),画图 ...
- [Sdoi2017]树点涂色 [lct 线段树]
[Sdoi2017]树点涂色 题意:一棵有根树,支持x到根染成新颜色,求x到y颜色数,求x子树里点到根颜色数最大值 考场发现这个信息是可减的,但是没想到lct 特意设计成lct的形式! 如何求颜色数? ...
- 【BZOJ4817】[Sdoi2017]树点涂色 LCT+线段树
[BZOJ4817][Sdoi2017]树点涂色 Description Bob有一棵n个点的有根树,其中1号点是根节点.Bob在每个点上涂了颜色,并且每个点上的颜色不同.定义一条路径的权值是:这条路 ...
- 【BZOJ4817】【SDOI2017】树点涂色 [LCT][线段树]
树点涂色 Time Limit: 10 Sec Memory Limit: 128 MB[Submit][Status][Discuss] Description Bob有一棵n个点的有根树,其中1 ...
随机推荐
- win10触摸板手势
尴尬的发现,今天才开始使用win10的手势,之前都是单指操作/笑哭 参考:http://www.sohu.com/a/63678223_230077 https://support.microsoft ...
- 2017-2018-2 20155231《网络对抗技术》实验五: MSF基础应用
2017-2018-2 20155231<网络对抗技术>实验五: MSF基础应用 实践目标 掌握信息搜集的最基础技能与常用工具的使用方法. 实验内容 (1)各种搜索技巧的应用 比如IP2L ...
- 20155317王新玮《网络对抗技术》实验9 web安全基础实践
20155317王新玮<网络对抗技术>实验9 web安全基础实践 一.实验准备 1.0 实验目标和内容 Web前端HTML.能正常安装.启停Apache.理解HTML,理解表单,理解GET ...
- WinRT IO相关整理
虽然一般UWP开发还是依赖.Net for UWP,但有时还是需要调用WinRT API.特别是在IO部分,WinRT有着和.Net似曾相识但又不尽相同的接口.在此对经常用到的一些地方进行一下整理. ...
- .NetCore利用BlockingCollection实现简易消息队列
前言 消息队列现今的应用场景越来越大,常用的有RabbmitMQ和KafKa. 我们用BlockingCollection来实现简单的消息队列. 实现消息队列 用Vs2017创建一个控制台应用程序.创 ...
- 浅谈android Service和BroadCastReceiver
1.题记 Android中的服务和windows中的服务是类似的东西,服务一般没有用户操作界面,它运行于系统中不容易被用户发觉,可以使用它开发如监控之类的程序. 广播接收者(BroadcastRece ...
- Salesforce随笔: 解决被指定给Chatter相关用户的RecordType无法被删除的问题
被指定给以下三组用户的RecordType无法在对应的Profile里取消占用: Chatter External User Chatter Free User Chatter Moderator U ...
- 《杜增强讲Unity之Tanks坦克大战》1-准备工作
0.案例介绍 0.1开始界面 点击Play Now 进入游戏界面 左边的坦克使用ws控制前后移动,ad键左右旋转,空格键开火 右边的坦克使用方向键上下控制前后移动,方向键左右键实现左右旋转 ...
- 20135202闫佳歆--week4 两种方式使用同一个系统调用--实验及总结
实验四 使用库函数API和C代码中嵌入汇编代码两种方式使用同一个系统调用 在这里我选择的是第20号系统调用,getpid. 1.使用库函数API: 代码如下: /* getpid.c */ #incl ...
- linux内核分析--计算机是如何工作的
实验部分 使用gcc -S -o main.s main.c -m32命令将源代码编译成汇编代码. 源代码如下: int g(int x) { return x + 9; } int f(int x) ...