Query on a tree V VI VII [LCT]

最新推荐文章于 2020-01-17 20:31:58 发布
原创 最新推荐文章于 2020-01-17 20:31:58 发布 · 283 阅读 · 0 ·
本内容遵循CC 4.0 BY-SA版权协议
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
收录于
FSY的好题汇总

当前文章被收录于：

FSY的好题汇总
97 篇文章 2 人学习

FSY暑假写的比较认真的博客就放在这儿了吧

订阅专栏 查看详情
LCT

当前文章被收录于：

LCT
19 篇文章 0 人学习
订阅专栏 查看详情
本文深入探讨了QTREE5、QTREE6和QTREE7算法，详细解释了使用Splay树进行点颜色反转、查询最近白点距离、同色连通块大小及最大值等操作的实现原理与代码细节。

QTREE5

操作 : 反转某个点的颜色, 查询某个点的最近白点

我们维护Splay中, 最浅节点到最近白点的距离(lsum), 以及最深节点到它的最近白点的距离(rsum)

那么查询x的答案很明显 Accss(x), Splay(x), x -> rsum 就是答案

考虑如何 Pushup

lsum_x=min(lsum_{ls}, siz_{ls}+min(Minson_x, lsum_{rs}+1))

Minson是虚子树答案, 用一个set维护

#include<bits/stdc++.h>
#define N 100050
using namespace std;
int read(){
	int cnt = 0; char ch = 0;
	while(!isdigit(ch)) ch = getchar();
	while(isdigit(ch)) cnt = cnt*10 + (ch-'0'), ch = getchar();
	return cnt;
}
const int inf = 0x3fffffff;
vector<int> v[N];
int n, m, col[N];
int ch[N][2], fa[N], siz[N], lsum[N], rsum[N];
multiset <int> S[N]; 
void dfs(int u, int f){
	for(int i=0; i<v[u].size(); i++){
		int t = v[u][i]; if(t == f) continue;
		fa[t] = u; S[u].insert(inf + 1); dfs(t, u);
	} 
}
#define ls ch[x][0]
#define rs ch[x][1]
bool isRoot(int x){ return ch[fa[x]][0] != x && ch[fa[x]][1] != x;}
int Top(int x){ if(S[x].empty()) return inf; return *S[x].begin();}
void Pushup(int x){
	if(!x) return; lsum[x] = rsum[x] = inf;
	siz[x] = siz[ls] + siz[rs] + 1;
	lsum[x] = min(lsum[ls], siz[ls] + min(col[x] ? 0 : inf, min(Top(x), lsum[rs] + 1)));
	rsum[x] = min(rsum[rs], siz[rs] + min(col[x] ? 0 : inf, min(Top(x), rsum[ls] + 1)));
}
void rotate(int x){
	int y = fa[x], z = fa[y], k = ch[y][1] == x;
	if(!isRoot(y)) ch[z][ch[z][1] == y] = x; fa[x] = z;
	ch[y][k] = ch[x][k^1]; fa[ch[x][k^1]] = y;
	ch[x][k^1] = y; fa[y] = x; Pushup(y); Pushup(x);
}
void Splay(int x){
	while(!isRoot(x)){
		int y = fa[x], z = fa[y];
		if(!isRoot(y)) rotate((ch[y][0] == x) ^ (ch[z][0] == y) ? x : y);
		rotate(x);
	}
}
void Access(int x){
	for(int y = 0; x; y = x, x = fa[x]){
		Splay(x); 
		if(rs) S[x].insert(lsum[rs] + 1);
		rs = y;
		if(rs) S[x].erase(S[x].find(lsum[rs] + 1));
		Pushup(x);
	}
}
int main(){
	n = read(); lsum[0] = rsum[0] = inf;
	for(int i=1; i<n; i++){
		int x = read(), y = read();	
		v[x].push_back(y); v[y].push_back(x); 
	} dfs(1, 0);
	m = read();
	while(m--){
		int op = read(), x = read();
		if(op == 0){ Access(x); Splay(x); col[x] ^= 1; Pushup(x);}
		if(op == 1){ Access(x); Splay(x); printf("%d\n", (rsum[x] == inf) ? -1 : rsum[x]);}
	} return 0;
}

QTREE6 

操作 : 反转颜色, 查询同色连通块的大小

首先开两颗LCT, 反转颜色相当于在一边暴力 Cut, 一边暴力Link

我们考虑只Cut x与它的fa, 也之Link x与它的fa, 这样一个连通块只有可能 rt与其颜色不同, 特判一下就好

然后就是维护虚子树的Siz

#include<bits/stdc++.h>
#define N 1000050
using namespace std;
int read(){
	int cnt = 0; char ch = 0;
	while(!isdigit(ch)) ch = getchar();
	while(isdigit(ch)) cnt = cnt*10 + (ch-'0'), ch = getchar();
	return cnt;
}
int first[N], nxt[N], to[N], tot;
void add(int x, int y){ nxt[++tot] = first[x], first[x] = tot, to[tot] = y;}
int n, m, col[N];
struct LCT{
	struct Node{ int ch[2], fa, siz; } t[N];
	int si[N];
	#define ls t[x].ch[0] 
	#define rs t[x].ch[1]
	bool isRoot(int x){ int fa = t[x].fa; return t[fa].ch[0] != x && t[fa].ch[1] != x; }
	void Pushup(int x){ t[x].siz = t[ls].siz + t[rs].siz + si[x] + 1;}
	void rotate(int x){
		int y = t[x].fa, z = t[y].fa, k = t[y].ch[1] == x;
		if(!isRoot(y)) t[z].ch[t[z].ch[1] == y] = x; t[x].fa = z;
		t[y].ch[k] = t[x].ch[k^1]; t[t[x].ch[k^1]].fa = y;
		t[x].ch[k^1] = y; t[y].fa = x; Pushup(y); Pushup(x);
	}
	void Splay(int x){
		while(!isRoot(x)){
			int y = t[x].fa, z = t[y].fa;
			if(!isRoot(y)) rotate((t[y].ch[0] == x) ^ (t[z].ch[0] == y) ? x : y);
			rotate(x);
		} Pushup(x);
	}
	void Access(int x){
		for(int y=0; x; y=x, x=t[x].fa){ 
			Splay(x);
			si[x] += t[rs].siz;
			rs = y;
			si[x] -= t[rs].siz;
			Pushup(x);
		} 
	}
	void Link(int x, int y){ 
		if(!y) return;
		Access(y); Splay(y); Splay(x); t[x].fa = y; si[y] += t[x].siz; Pushup(y); 
	}
	void Cut(int x, int y){
		if(!y) return;
		Access(x); Splay(x); t[ls].fa = 0; ls = 0; Pushup(x);
	}
	int Findroot(int x){
		Access(x); Splay(x);
		while(ls) x = ls; Splay(x);
		return x; 
	}
	int Quary(int x){
		int c = col[x]; x = Findroot(x);
		return col[x] == c ? t[x].siz : t[rs].siz;
	}
}lct[2]; 
int fa[N];
void dfs(int u, int f){
	for(int i=first[u];i;i=nxt[i]){
		int t = to[i]; if(t == f) continue;
		fa[t] = u; lct[0].Link(t, u); dfs(t, u);
	}
}
int main(){
	n = read();	
	for(int i=1; i<n; i++){
		int x = read(), y = read();
		add(x, y); add(y, x);
	} dfs(1, 0); m = read();
	while(m--){
		int op = read(), x = read();
		if(op == 0) printf("%d\n", lct[col[x]].Quary(x));
		if(op == 1){ lct[col[x]].Cut(x, fa[x]); col[x] ^= 1; lct[col[x]].Link(x, fa[x]);}
	} return 0;
}

QTREE7 

与上一道题类似, 维护同色连通块的最大值, multiset就可以了

#include<bits/stdc++.h>
#define N 200050
using namespace std;
int read(){
	int cnt = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1;}
	while(isdigit(ch)) cnt = cnt*10 + (ch-'0'), ch = getchar();
	return cnt * f;
}
int first[N], nxt[N], to[N], tot;
void add(int x, int y){nxt[++tot] = first[x], first[x] = tot, to[tot] = y;}
int n, m, col[N], f[N], w[N];
struct Node{
	int ch[N][2], fa[N], Max[N];
	multiset<int> S[N];
	bool isRoot(int x){ return ch[fa[x]][0] != x && ch[fa[x]][1] != x;}
	void Pushup(int x){
		Max[x] = max(w[x], max(Max[ch[x][0]], Max[ch[x][1]]));
		if(S[x].size()) Max[x] = max(Max[x], *S[x].rbegin());
	}
	void rotate(int x){
		int y = fa[x], z = fa[y], k = ch[y][1] == x;
		if(!isRoot(y)) ch[z][ch[z][1] == y] = x; fa[x] = z;
		ch[y][k] = ch[x][k^1]; fa[ch[x][k^1]] = y;
		ch[x][k^1] = y; fa[y] = x; Pushup(y); Pushup(x);
	}
	void Splay(int x){
		while(!isRoot(x)){
			int y = fa[x], z = fa[y];
			if(!isRoot(y)) rotate((ch[y][1]==x) ^ (ch[z][1]==y) ? x : y);
			rotate(x); 
		} Pushup(x);
	}
	void Access(int x){
		for(int y = 0; x; y = x, x = fa[x]){
			Splay(x); 
			if(ch[x][1]) S[x].insert(Max[ch[x][1]]);
			ch[x][1] = y; 
			if(ch[x][1]) S[x].erase(Max[ch[x][1]]);
			Pushup(x);
		}
	}
	int Findroot(int x){
		Access(x); Splay(x); while(ch[x][0]) x = ch[x][0]; Splay(x); return x;
	}
	void Link(int x, int y){
		if(!y) return;
		Access(y); Splay(y); Splay(x);
		fa[x] = y; S[y].insert(Max[x]); Pushup(y);
	}
	void Cut(int x, int y){
		if(!y) return;
		Access(x); Splay(x); fa[ch[x][0]] = 0; ch[x][0] = 0; Pushup(x);
	}
	int Quary(int x){
		int c = col[x]; x = Findroot(x);
		return col[x] == c ? Max[x] : Max[ch[x][1]];
	}
	void Modify(int x){
		Access(x); Splay(x); w[x] = read(); Pushup(x);
	}
}T[2];
void dfs(int u, int ff){
	for(int i=first[u];i;i=nxt[i]){
		int t = to[i]; if(t == ff) continue;
		T[col[t]].Link(t, u); f[t] = u; dfs(t, u); 
	}
}
int main(){
	n = read();
	for(int i=1; i<n; i++){
		int x = read(), y = read();
		add(x, y); add(y, x);
	} 
	for(int i=1; i<=n; i++) col[i] = read();
	for(int i=1; i<=n; i++) w[i] = read();
	T[0].Max[0] = T[1].Max[0] = -2e9;
	dfs(1, 0); m = read();
	while(m--){
		int op = read(), u = read();
		if(op == 0) printf("%d\n", T[col[u]].Quary(u));
		if(op == 1){ T[col[u]].Cut(u, f[u]); col[u] ^= 1; T[col[u]].Link(u, f[u]);}
		if(op == 2){ T[col[u]].Modify(u);}
	} return 0;
}

 

SP16580 QTREE7 - Query on a tree VII（LCT 维护树上同色点连通块 + LCT维护子树最大点权 + multiset）
qq_41997978的博客
02-26 404
建立两棵 LCT，分别维护颜色为 0 的图和颜色为 1 的图，维护的方法同 SP16549 QTREE6 - Query on a tree VI。 对 LCT 上的每个节点，要维护子树点权最大值，虚子树的信息用 multiset 维护，对于操作3要修改两棵 LCT 的点权，而不是一棵。 如果 pushup 采用两段直接合并，注意0号点的 valvalval 和 mxvalmxvalmxval...
BZOJ 3639: Query on a tree VII LCT_set维护子树信息
EM LGH
12-22 174
Code: #include &lt;cstring&gt; #include &lt;cstdio&gt; #include &lt;algorithm&gt; #include &lt;set&gt; #include &lt;string&gt; #define setIO(s) freopen(s".in","r",stdin), freopen(s".ou
[SP16580 QTREE7 - Query on a tree VII]
LPA的博客
05-24 418
洛谷传送门 题目描述 You are given a tree (an acyclic undirected connected graph) with n nodes. The tree nodes are numbered from 111 to nnn. Each node has a color, white or black, and a weight. We will ask yo...
LCT 动态维护 LCA
FSYo 的博客
02-12 1282
步骤如下: 1. 将当前查询需要的根makeroot, 表示为当前根下的lca 2. 将x Access到根, 此时x到根的路径是一棵平衡树 3. 将y Access, 中途遇到了与根一棵平衡树的点就记录下来(即与x到根的路径相遇), 返回即可  模板传送门 #include&lt;bits/stdc++.h&gt; #define N 500050 #define ls t[x].c...
近期刷题总结 [19 03 17]
FSYo 的博客
03-17 1093
目录 P3703 [SDOI2017]树点涂色[LCT + 线段树] CF739B Alyona and a tree[主席树] P4098 [HEOI2013]ALO[双向链表 + 可持久化0/1trie] P5022 旅行[基环树] P1453 城市环路[基环树DP] CF837G Functions On The Segments[主席树] CF837D Round...
LCT 进阶指南
FSYo 的博客
04-21 722
这里总结了LCT的进阶应用, 以及几道比较难的例题 强烈推荐 :https://www.cnblogs.com/flashhu/p/9498517.html 目录 LCT 维护双强连通分量 LCT 维护子树信息 例题 : P2542 [AHOI2005]航线规划 U19464 山村游历(Wander) P4299 首都 BZOJ4998星球联盟 BZOJ2959长跑 P3...
近期刷题总结 [19 03 24]
FSYo 的博客
03-24 705
目录 CF14D Two Paths[树的直径] P3174 [HAOI2009]毛毛虫[树的直径] P2312 解方程[秦九韶算法] CF888D Almost Identity Permutations[ DP ] CF888G Xor-MST[0/1 trie] CF863D Yet Another Array Queries Problem[ splay 模板 ] ...
城市旅行[LCT]
FSYo 的博客
10-10 593
题目描述 W国地人物博，有n座城市组成，共n-1条双向道路连接其中的两座城市，且任意两座城市都可相互到达。 风景秀美的w国吸引了无数慕名而来的游客，根据游客对每座城市的打分，我们定义第i座城市的美丽度为a_i。一次从城市x到城市y的旅行，所获得的的偷悦指数为从城市x到城市y所有城市的美丽度之和（包括X和y）。我们诅义这个值为H(x,y)。 现在小A在城市X，Sharon在城市Y，他们想知道如...
NOIP 2018 保卫王国 （动态 DP）（LCT）
FSYo 的博客
11-13 358
传送门 跟动态 DP 的模板差不多 考虑朴素的树形 dpdpdp fu,0=∑vfv,1f_{u,0}=\sum_{v}f_{v,1}fu,0​=v∑​fv,1​ fu,1=val[u]+∑vmin(fv,0,fv,1)f_{u,1}=val[u]+\sum_{v}min(f_{v,0},f_{v,1})fu,1​=val[u]+v∑​min(fv,0​,fv,1​) 令 gu,0/1g_{u,0...
[SCOI2018] Tree [LCT]
FSYo 的博客
07-08 348
传送门 类似 QTREE 的套路, 我们维护到Splay最浅点的最长路径Lmax, 最深点的最长路径Rmax 考虑如何Pushup? son 是虚子树, 用一个 multiset 维护虚子树 #include<bits/stdc++.h> #define N 200050 using namespace std; typedef long long ll; const ll...
近期刷题总结[2019 03 09]
FSYo 的博客
03-09 343
目录 寻找段落[分数规划+单调队列] P4319 变化的道路[LCT + 线段树分治] P2485 [SDOI2011]计算器[exgcd, bsgs] [模板] P4458 [BJOI2018]链上二次求和[线段树][二次前缀和] 严格次小生成树[最小生成树+倍增] P2120 [ZJOI2007]仓库建设[斜率优化板子] P2155 [SDOI2008]沙拉公主的困惑...
弹飞绵羊[LCT]
FSYo 的博客
10-08 324
题目描述 某天，Lostmonkey发明了一种超级弹力装置，为了在他的绵羊朋友面前显摆，他邀请小绵羊一起玩个游戏。游戏一开始，Lostmonkey在地上沿着一条直线摆上n个装置，每个装置设定初始弹力系数ki，当绵羊达到第i个装置时，它会往后弹ki步，达到第i+ki个装置，若不存在第i+ki个装置，则绵羊被弹飞。绵羊想知道当它从第i个装置起步时，被弹几次后会被弹飞。为了使得游戏更有趣，Lostmo...
差值最小生成树[LCT]
FSYo 的博客
10-29 310
传送门 将边按边权从小到大排序; 如果插入的两个点在一棵树上,就找他们路径上权值最小的那条边并删除,然后在新加一条边 否则直接加就好了 路径需要以点的形式插入到LCT,第i跳路径对于i+n号节点,Link(x,y) 就Link(x,i+n),Link(i+n,y) 就好了 维护一下最小值 #include&lt;bits/stdc++.h&gt; #define N 1000005 ...
暑假好题汇总
FSYo 的博客
07-07 303
目录 [JSOI2007]文本生成器 [AC自动机 + DP] BZOJ3687简单题[Bitset] 恨 7 不成妻[数位 DP] [SHOI2011]双倍回文[Manacher] [NOI2007]货币兑换[CDQ+斜率优化DP] BZOJ2655 calc[ 生成函数 + DP + 拉格朗日差值 ] P3401 洛谷树[树链剖分] [BZOJ4205][WOJ3...
极地旅行社[LCT]
FSYo 的博客
10-23 293
题目描述 不久之前，Mirko建立了一个旅行社，名叫“极地之梦”。这家旅行社在北极附近购买了N座冰岛，并且提供观光服务。 当地最受欢迎的当然是帝企鹅了，这些小家伙经常成群结队的游走在各个冰岛之间。Mirko的旅行社遭受一次重大打击，以至于观光游轮已经不划算了。旅行社将在冰岛之间建造大桥，并用观光巴士来运载游客。 Mirko希望开发一个电脑程序来管理这些大桥的建造过程，以免有不可预料的错误发生...
【51nod1600】Simple KMP（SAM）（LCT）（差分）
FSYo 的博客
01-17 289
题意： 对于一个字符串 ∣S∣|S|∣S∣，我们定义 fail[i]fail[i]fail[i]，表示最大的 x 使得 S[1..x]=S[i−x+1..i]S[1..x]=S[i-x+1..i]S[1..x]=S[i−x+1..i]，满足 (x<i)(x<i)(x<i) 显然对于一个字符串，如果我们将每个 0<=i<=∣S∣0<=i<=|S|0<=...
P2486 [SDOI2011]染色 [LCT]
FSYo 的博客
02-13 257
传送门 修改时提取路径打tag就可以 询问时提取路径直接返回splay根节点的值 值得一提的是Pushup, 我们发现如果x 与 x的前驱颜色相同, 那么就算一条路径, 后继同理 为了高效地找到前驱后继的颜色, 我们维护一个lcol , rcol 表示前驱后继的颜色 t[x].lcol = t[ls].lcol t[x].rcol = t[rs].rcol 如果 t[x].col =...
实用代码脚本易语言源码时间
最新发布
07-26
实用代码脚本易语言源码时间
实用代码脚本易语言源码设置控件锚点
07-26
实用代码脚本易语言源码设置控件锚点
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值