P3703-[SDOI2017]树点涂色【LCT,线段树】

正题

题目链接:https://www.luogu.com.cn/problem/P3703


题目大意

nnn个点的一棵树开始所有点有不同的颜色,mmm次操作

  1. 将根节点到xxx节点的路径上染上一种新的颜色
  2. 询问一条路径的不同颜色个数
  3. 询问一个节点的子树中的一个xxx使得xxx到根节点的颜色最多。

解题思路

操作111LCTLCTLCTaccessaccessaccess操作很相似。相同颜色之间就是实边,不同颜色之间就是虚边。

操作222就是之间px+py−2pLCA+1p_x+p_y-2p_{LCA}+1px+py2pLCA+1就好了,但是考虑到操作333,所以维护一个dfndfndfn序和线段树就可以查询子树最大值了。

之后维护一个LCTLCTLCT,在accessaccessaccess操作切换虚实边的时候修改一下线段树就好了,并且需要注意我们不能直接拿SplaySplaySplay的根的子树,要找到实际的树中的根,所以SplaySplaySplay一直往左就好了。

好像还有树链剖分的做法,线段树查询的时候维护一下末尾颜色好像就可以了,这里不多讲(我也不会)

时间复杂度O(nlog⁡2n)O(n\log^2 n)O(nlog2n)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e5+10;
struct node{int to,next;
}a[N<<1];
int n,m,cnt,tot,ls[N],rfn[N],ed[N],fa[N];
int dep[N],son[N],siz[N],top[N];
struct SegTree{int w[N<<2],lazy[N<<2];void Downdata(int x){if(!lazy[x])return;w[x*2]+=lazy[x];lazy[x*2]+=lazy[x];w[x*2+1]+=lazy[x];lazy[x*2+1]+=lazy[x];lazy[x]=0;return;}void Change(int x,int L,int R,int l,int r,int val){if(L==l&&R==r){w[x]+=val;lazy[x]+=val;return;}int mid=(L+R)>>1;Downdata(x);if(r<=mid)Change(x*2,L,mid,l,r,val);else if(l>mid)Change(x*2+1,mid+1,R,l,r,val);else Change(x*2,L,mid,l,mid,val),Change(x*2+1,mid+1,R,mid+1,r,val);w[x]=max(w[x*2],w[x*2+1]);return;}int Ask(int x,int L,int R,int l,int r){if(L==l&&R==r)return w[x];int mid=(L+R)>>1;Downdata(x);if(r<=mid)return Ask(x*2,L,mid,l,r);if(l>mid)return Ask(x*2+1,mid+1,R,l,r);return max(Ask(x*2,L,mid,l,mid),Ask(x*2+1,mid+1,R,mid+1,r));}
}Tr;
struct LinkCutTree{int t[N][2],fa[N];bool Nroot(int x){return fa[x]&&(t[fa[x]][0]==x||t[fa[x]][1]==x);}bool Direct(int x){return t[fa[x]][1]==x;}void Rotate(int x){int y=fa[x],z=fa[y];int xs=Direct(x),ys=Direct(y);int w=t[x][xs^1];if(Nroot(y))t[z][ys]=x;t[x][xs^1]=y;t[y][xs]=w;if(w)fa[w]=y;fa[y]=x;fa[x]=z;return;}void Splay(int x){while(Nroot(x)){int y=fa[x];if(!Nroot(y))Rotate(x);else if(Direct(x)==Direct(y))Rotate(y),Rotate(x);else Rotate(x),Rotate(x);}return;}int FindRoot(int x){while(t[x][0])x=t[x][0];return x;}void Access(int x){for(int y=0;x;y=x,x=fa[x]){Splay(x);int z=t[x][1];if(z)z=FindRoot(z),Tr.Change(1,1,n,rfn[z],ed[z],1);if(y)z=FindRoot(y),Tr.Change(1,1,n,rfn[z],ed[z],-1);t[x][1]=y;}return;}
}T;
void addl(int x,int y){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;return;
}
void dfs1(int x){siz[x]=1;rfn[x]=++cnt;dep[x]=dep[fa[x]]+1;T.fa[x]=fa[x];Tr.Change(1,1,n,cnt,cnt,dep[x]);for(int i=ls[x];i;i=a[i].next){int y=a[i].to;if(y==fa[x])continue;fa[y]=x;dfs1(y);siz[x]+=siz[y];if(siz[y]>siz[son[x]])son[x]=y;}ed[x]=cnt;
}
void dfs2(int x){if(son[x]){top[son[x]]=top[x];dfs2(son[x]);}for(int i=ls[x];i;i=a[i].next){int y=a[i].to;if(y==fa[x]||y==son[x])continue;top[y]=y;dfs2(y);}return;
}
int LCA(int x,int y){while(top[x]!=top[y]){if(dep[top[x]]<dep[top[y]])swap(x,y);x=fa[top[x]];}return dep[x]<dep[y]?x:y;
}
int main()
{// freopen("paint1.in","r",stdin);scanf("%d%d",&n,&m);for(int i=1;i<n;i++){int x,y;scanf("%d%d",&x,&y);addl(x,y);addl(y,x);}dfs1(1);dfs2(1);while(m--){int op,x,y;scanf("%d%d",&op,&x);if(op==1)T.Access(x);else if(op==2){scanf("%d",&y);int lca=LCA(x,y);int p1=Tr.Ask(1,1,n,rfn[x],rfn[x]);int p2=Tr.Ask(1,1,n,rfn[y],rfn[y]);int p3=Tr.Ask(1,1,n,rfn[lca],rfn[lca]);printf("%d\n",p1+p2-p3*2+1);}else printf("%d\n",Tr.Ask(1,1,n,rfn[x],ed[x]));}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/319602.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

牛客练习赛 68

A.牛牛的mex n,q≤105,0≤ai<n且ai互不相同n,q≤10 ^5,0≤a i<n 且 a_i互不相同n,q≤105,0≤ai<n且ai​互不相同 后面两个条件非常重要&#xff0c;通过后面两个条件将问题转化为一个区间内最小未出现的自然数就等于不在这个区间内最小出现的自然数对于区间[l,r][l,r]…

牛客网【每日一题】7月31日题目精讲—兔子的区间密码

来源&#xff1a;牛客网 &#xff1a; 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 32768K&#xff0c;其他语言65536K 64bit IO Format: %lld文章目录题目描述题解&#xff1a;代码&#xff1a;题目描述 有一只可爱的兔子被困在了密室了&a…

Redis 实战之主从复制、高可用、分布式

简介本节内容基于 CentOS 7.4.1708&#xff0c;Redis 3.2.12 环境实验。Redis 是一个开源的高性能键值对数据库。安装&#xff1a;yum install -y redis特性&#xff1a;高性能 Key-Value 服务器多种数据结构丰富功能缓存(get|set)计数器(incre)消息队列(publish|subcrib)高可用…

纪中A组模拟赛总结(2021.7.16)

成绩 rankrankranknamenamenamescorescorescoreT1T1T1T2T2T2T3T3T3T4T4T4525252lyflyflyf000000000000000 前言&#xff1a;暴零了&#xff0c;太棒了呀(&#xff3e;&#xff0d;&#xff3e;) 总结 开考后再四题里来回看&#xff0c;都没有找到解题方法&#xff08;主要是有…

CF710F-String Set Queries【AC自动机,二进制分组】

正题 题目链接:https://www.luogu.com.cn/problem/CF710F 题目大意 TTT次操作 往集合中加入一个字符串往集合中删除一个字符串给出一个模式串求出现的集合里面的字符串个数 解题思路 删除的话改成加入一个权值为−1-1−1的字符串就是全都是加入操作了。 然后就可以像[SDOI2…

牛客网 【每日一题】8月5日题目精讲—蓝魔法师

来源&#xff1a;牛客网&#xff1a; 文章目录题目描述题解&#xff1a;代码&#xff1a;题目描述 “你&#xff0c;你认错人了。我真的&#xff0c;真的不是食人魔。”–蓝魔法师 给出一棵树&#xff0c;求有多少种删边方案&#xff0c;使得删后的图每个连通块大小小于等于k…

AtCoder Beginner Contest 178 总结

A - Not 签到题 #define IO ios::sync_with_stdio(false);cin.tie();cout.tie(0) #pragma GCC optimize(2) #include<iostream> #include<algorithm> using namespace std; const int N100010; int main() {IO;int T1;//cin>>T;while(T--){int x;cin>&g…

.NET Core 2.1.5和.NET Core SDK 2.1.403发布

2018 年 10 月 1 日&#xff0c;.NET Core 2.0 已宣布结束生命周期。这意味着 2.0 系列将不会提供更新。2.1 系列上月作为长期支持版发布&#xff0c;.NET Core 2.1.5 和 .NET Core SDK 2.1.403 版本包含许多问题修复。主要更新内容如下&#xff1a;CoreCLR2018-09-11 - [cb730…

P4859-已经没有什么好害怕的了【容斥,dp】

正题 题目链接:https://www.luogu.com.cn/problem/P4859 题目大意 两个长度为nnn的序列a,ba,ba,b两两匹配&#xff0c;求ai>bia_i>b_iai​>bi​的组数比ai<bia_i<b_iai​<bi​的组数多kkk的方案数。 保证输入数字两两不同 解题思路 其实就是求恰好有nk2\f…

(牛客网)树型dp

树型dp 视频链接 &#xff08;如果想购买网课&#xff0c;可以用我的邀请码&#xff09; 用我的链接购买&#xff0c;我再反你10&#xff0c;一共花54多值 购买链接 不放心可以先加我好友2830872914 总试题链接 文章目录树型dp例题NC15033 小G有一个大树NC511788 没有上司的舞…

Node 源项目定制化、打包并使用全过程讲解

&#x1f468;&#x1f3fb;‍&#x1f4bb; 热爱摄影的程序员 &#x1f468;&#x1f3fb;‍&#x1f3a8; 喜欢编码的设计师 &#x1f9d5;&#x1f3fb; 擅长设计的剪辑师 &#x1f9d1;&#x1f3fb;‍&#x1f3eb; 一位高冷无情的编码爱好者 大家好&#xff0c;我是全栈工…

纪中A组模拟赛总结(2021.7.17)

成绩 rankrankranknamenamenamescorescorescoreT1T1T1T2T2T2T3T3T3222222lyflyflyf999999000000总结 看到T1每一段有代价&#xff0c;每条边都至少要便利1次&#xff0c;感觉是树上线性规划&#xff08;昨天刚讲线性规划&#xff09;&#xff0c;觉得正着没那么好做&#xff0…

.Net Core中的Api版本控制

原文链接&#xff1a;API Versioning in .Net Core作者&#xff1a;Neel Bhatt简介Api的版本控制是Api开发中经常遇到的问题, 在大部分中大型项目都需要使用到Api的版本控制在本篇博客中&#xff0c;我们将说明一下如何在.Net Core Api项目中使用Api版本控制。本篇博客中测试项…

Educational Codeforces Round 95 (Rated for Div. 2)

昨天本来想打一下&#xff0c;但是今天早上课很早&#xff0c;就没有打&#xff0c;只是看了看前三个题写了个代码&#xff0c;今天中午课结束交了一下都AC了。y1s1 A题第一次写就出来了&#xff0c;但是答案一直不对&#xff0c;最后结果是样例错了-。- A - Buying Torches …

【每日一题】8月4日题目精讲—购物

来源&#xff1a;牛客网&#xff1a; 文章目录购物题目描述题解&#xff1a;代码&#xff1a;购物 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 32768K&#xff0c;其他语言65536K 64bit IO Format: %lld题目描述 在遥远的东方&#xff0c;有…

P4606-[SDOI2018]战略游戏【圆方树,虚树】

正题 题目链接:https://www.luogu.com.cn/problem/P4606 题目大意 给出nnn个点mmm条边的一张图&#xff0c;qqq次询问给出一个点集&#xff0c;询问有多少个点割掉后可以是点集中至少一个点对不连通。 解题思路 就是问圆方树上的虚树中的圆点数量&#xff0c;照着统计就好了…

【二分】雪(luogu 7405)

正题 luogu 7405 题目大意 坐标轴上有n个雪球&#xff0c;初始重量为0&#xff0c;每一条线段上有重量为1的雪&#xff0c;当雪球经过时&#xff0c;会加上这些雪&#xff0c;而地上就没有雪了 共有m个时刻&#xff0c;每个时刻会使所有雪球向左/右移动wiw_iwi​格&#xff…

Visual Studio 2017 与 Visual Studio for Mac 支持更新

微软在博客中简单介绍了关于 VS 2017 和 VS for Mac 项目的支持计划&#xff1a;https://blogs.msdn.microsoft.com/visualstudio/2018/10/05/visual-studio-2017-and-visual-studio-for-mac-support-updates/。微软表示&#xff0c;在目前努力开发 Visual Studio 2019 的同时&…

牛客练习赛 67——ST表

A.牛牛爱字符串 注意原字符串有空格&#xff0c;不要用cin #define IO ios::sync_with_stdio(false);cin.tie();cout.tie(0) #pragma GCC optimize(2) #include<string> #include<iostream> #include<algorithm> using namespace std; int main() {IO;int …

【每日一题】8月3日题目精讲—小A的最短路

来源&#xff1a;牛客网&#xff1a; 文章目录题目描述题解&#xff1a;代码&#xff1a;时间限制&#xff1a;C/C 3秒&#xff0c;其他语言6秒 空间限制&#xff1a;C/C 262144K&#xff0c;其他语言524288K 64bit IO Format: %lld题目描述 小A这次来到一个景区去旅游&#xf…