P4211 [LNOI2014]LCA & CF757G Can Bash Save the Day? 解题报告

最新推荐文章于 2024-06-21 23:28:17 发布
最新推荐文章于 2024-06-21 23:28:17 发布
阅读量116 收藏
点赞数
分类专栏： OI的那些事 文章标签： 树链剖分 数据结构 树 可持久化
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41996523/article/details/112852806
版权

P4211 [LNOI2014]LCA & CF757G Can Bash Save the Day? 解题报告

这两道题很有代表性，代表了一种有关lca问题的解法。

P4211 [LNOI2014]LCA

链接

题意很清晰，就不转化了。

要求 ∑ i = l r d e p [ l c a ( i , z ) ] \sum_{i=l}^r\mathrm{dep}[\mathrm{lca}(i,z)] ∑i=lr​dep[lca(i,z)]，考虑差分。只要能求 ∑ i = 1 x d e p [ l c a ( i , z ) ] \sum_{i=1}^x\mathrm{dep}[\mathrm{lca}(i,z)] ∑i=1x​dep[lca(i,z)] 就行了。

这里用到一个非常巧妙的转化， lca的深度就是两点到根的重叠部分的长度。于是，只要将 1 ~ x x x 到根的路径全部加1，再查询 z z z 到根的路径上的值即可。这可以用树链剖分完成。

只要离散+树剖线段树维护即可。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
char In[1 << 20], *ss = In, *tt = In;
#define getchar() (ss == tt && (tt = (ss = In) + fread(In, 1, 1 << 20, stdin), ss == tt) ? EOF : *ss++)
ll read() {
	ll 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 + int(ch - '0');
	return x * f;
}
const int MAXN = 50005, MAXM = 50005, P = 201314;
int n, q, fa[MAXN], ver[MAXN], nxt[MAXN], cnt, head[MAXN], tot, dfn[MAXN], tim, sz[MAXN], son[MAXN], top[MAXN], dep[MAXN];
ll ans[MAXM];
struct Qry {int u, z, id, k;}qry[MAXM << 1];
bool operator < (const Qry& a, const Qry& b) {return a.u < b.u;}

void addedge(int u, int v) {
	ver[++cnt] = v; nxt[cnt] = head[u]; head[u] = cnt;
}
//Segment_tree
#define ls o << 1
#define rs o << 1 | 1
struct SgtNode {ll v, a;}e[MAXN << 2];
void pushup(int o) {e[o].v = e[ls].v + e[rs].v;}
void addpoint(int o, int l, int r, ll k) {e[o].v += (r-l+1)*k; e[o].a += k;}
void pushdown(int o, int l, int r) {
	if(e[o].a) {
		int m = (l + r) >> 1;
		addpoint(ls, l, m, e[o].a); addpoint(rs, m+1, r, e[o].a);
		e[o].a = 0;
	}
}
void addrange(int o, int l, int r, int x, int y, ll k) {
	if(x <= l && r <= y) {addpoint(o, l, r, k); return ;}
	pushdown(o, l, r);
	int m = (l + r) >> 1;
	if(x <= m) addrange(ls, l, m, x, y, k);
	if(y > m) addrange(rs, m+1, r, x, y, k);
	pushup(o);
}
ll qryrange(int o, int l, int r, int x, int y) {
	if(x <= l && r <= y) return e[o].v;
	pushdown(o, l, r);
	int m = (l + r) >> 1;
	if(y <= m) return qryrange(ls, l, m, x, y);
	else if(x > m) return qryrange(rs, m+1, r, x, y);
	else return qryrange(ls, l, m, x, y) + qryrange(rs, m+1, r, x, y);
}
#undef ls 
#undef rs
//Tree cut
void dfs1(int u) {
	dep[u] = dep[fa[u]] + 1; sz[u] = 1; son[u] = 0;
	for(int i = head[u]; i; i = nxt[i]) {
		int v = ver[i]; dfs1(v); sz[u] += sz[v];
		if(sz[v] > sz[son[u]]) son[u] = v;
	}
}
void dfs2(int u, int tprt) {
	top[u] = tprt; dfn[u] = ++tim; if(son[u]) dfs2(son[u], tprt);
	for(int i = head[u]; i; i = nxt[i]) if(ver[i] != son[u]) {
		int v = ver[i]; dfs2(v, v);
	}
}
void addroute(int u, int v, ll k) {
	while(top[u] != top[v]) {
		if(dep[top[u]] < dep[top[v]]) swap(u, v);
		addrange(1, 1, n, dfn[top[u]], dfn[u], k);
		u = fa[top[u]];
	}
	if(dep[u] < dep[v]) swap(u, v);
	addrange(1, 1, n, dfn[v], dfn[u], k);
}
ll qryroute(int u, int v) {
	ll ans = 0;
	while(top[u] != top[v]) {
		if(dep[top[u]] < dep[top[v]]) swap(u, v);
		ans += qryrange(1, 1, n, dfn[top[u]], dfn[u]);
		u = fa[top[u]];
	}
	if(dep[u] < dep[v]) swap(u, v);
	ans += qryrange(1, 1, n, dfn[v], dfn[u]);
	return ans;
}


int main() {
	n = read(); q = read();
	for(int i = 2; i <= n; i++) {
		fa[i] = read()+1; addedge(fa[i], i);//father to son
	}
	for(int i = 1; i <= q; i++) {
		int l = read()+1, r = read()+1, z = read()+1;
		qry[++tot] = (Qry){r, z, i, 1}; 
		if(l != 1) qry[++tot] = (Qry){l-1, z, i, -1};
	}
	sort(qry + 1, qry + 1 + tot);
	dfs1(1); dfs2(1, 1);
	int cur = 1;
	for(int i = 1; i <= n && cur <= tot; i++) {
		addroute(i, 1, 1);
		for(; cur <= tot && qry[cur].u <= i; cur++) {
			ans[qry[cur].id] += qry[cur].k * qryroute(qry[cur].z, 1);
		}
	}
	for(int i = 1; i <= q; i++)
		printf("%lld\n", ans[i] % P);
	return 0;
}

CF757G Can Bash Save the Day

链接

一棵 n n n 个点的树和一个排列 p i p_i pi​，边有边权，支持两种操作：

  • l r x，询问 ∑ i = l r d i s ( p i , x ) \sum\limits_{i=l}^{r} dis(p_i,x) i=l∑r​dis(pi​,x)。
  • x，交换 p x , p x + 1 p_x,p_{x+1} px​,px+1​
  • n , q ≤ 2 × 1 0 5 n,q\le 2\times 10^5 n,q≤2×105，强制在线。

（题目描述摘自洛谷）

其实这道题和上道题没有太大区别，只是上道题可以离散化，这道题强制在线。

设 D i s ( x ) Dis(x) Dis(x) 表示 x x x 到根的距离。这里的根可以随意定。

依然考虑差分，则

∑ i = 1 t d i s ( p i , x ) = ∑ i = 1 t [ D i s ( p i ) + D i s ( x ) − 2 ⋅ D i s ( l c a ( p i , x ) ) ] \sum\limits_{i=1}^tdis(p_i,x)=\sum\limits_{i=1}^t[Dis(p_i)+Dis(x)-2\cdot Dis(\mathrm{lca}(p_i,x))] i=1∑t​dis(pi​,x)=i=1∑t​[Dis(pi​)+Dis(x)−2⋅Dis(lca(pi​,x))]

= ∑ i = 1 t D i s ( p i ) + t ⋅ D i s ( x ) − 2 ∑ i = 1 t D i s ( l c a ( p i , x ) ) =\sum\limits_{i=1}^tDis(p_i)+t\cdot Dis(x)-2\sum\limits_{i=1}^tDis(lca(p_i, x)) =i=1∑t​Dis(pi​)+t⋅Dis(x)−2i=1∑t​Dis(lca(pi​,x))

第一项可以前缀和维护，第二项可以直接 O ( 1 ) O(1) O(1) 算出，而第三项可以用类似于上题的方法得到。但是要用可持久化线段树维护。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
char In[1 << 20], *ss = In, *tt = In;
#define getchar() (ss == tt && (tt = (ss = In) + fread(In, 1, 1 << 20, stdin), ss == tt) ? EOF : *ss++)
ll read() {
	ll 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 + int(ch - '0');
	return x * f;
}

const int MAXN = 2e5 + 5;
const int MAXNODE = 3e7;
int n, q, a[MAXN], head[MAXN], ver[MAXN << 1], nxt[MAXN << 1], cnt, dep[MAXN], fa[MAXN], sz[MAXN], son[MAXN], top[MAXN], dfn[MAXN], tim;
ll edg[MAXN << 1], dist[MAXN], suma[MAXN], sw[MAXN], ew[MAXN];
void addedge(int u, int v, ll w) {
	ver[++cnt] = v; nxt[cnt] = head[u]; edg[cnt] = w; head[u] = cnt;
}

void dfs1(int u, int f) {
	fa[u] = f; dep[u] = dep[f] + 1; son[u] = 0; sz[u] = 1;
	for(int i = head[u]; i; i = nxt[i]) if(ver[i] != f) {
		int v = ver[i]; 
		dist[v] = dist[u] + edg[i]; ew[v] = edg[i]; dfs1(v, u); sz[u] += sz[v];
		if(sz[v] > sz[son[u]]) son[u] = v;
	}
}
void dfs2(int u, int tprt) {
	top[u] = tprt; dfn[u] = ++tim; sw[tim] = sw[tim-1] + ew[u]; if(son[u]) dfs2(son[u], tprt);
	for(int i = head[u]; i; i = nxt[i]) if(ver[i] != fa[u] && ver[i] != son[u]) {
		int v = ver[i]; dfs2(v, v);
	}
}

int tot, limtot, rt[MAXN];
struct SgtNode {ll v, a; int ls, rs;}e[MAXNODE];
void addrange(int &o, int pre, int l, int r, int x, int y) {
	if(pre <= limtot) o = ++tot, e[o] = e[pre];
	else o = pre;
	if(x <= l && r <= y) {e[o].a ++; e[o].v += sw[r] - sw[l-1]; return;}
	int m = (l + r) >> 1;
	if(x <= m) addrange(e[o].ls, e[pre].ls, l, m, x, y);
	if(y > m) addrange(e[o].rs, e[pre].rs, m+1, r, x, y);
	e[o].v = e[e[o].ls].v + e[e[o].rs].v + e[o].a * (sw[r] - sw[l-1]);
}
ll qryrange(int o, int l, int r, int x, int y) {
	if(x <= l && r <= y) return e[o].v;
	int m = (l + r) >> 1; ll ans = e[o].a * (sw[y] - sw[x-1]);
	if(y <= m) return ans + qryrange(e[o].ls, l, m, x, y);
	else if(x > m) return ans + qryrange(e[o].rs, m+1, r, x, y);
	else return ans + qryrange(e[o].ls, l, m, x, m) + qryrange(e[o].rs, m+1, r, m+1, y);
}
void addroute(int vs, int u, int v, ll k) {
	limtot = tot;
	while(top[u] != top[v]) {
		if(dep[top[u]] < dep[top[v]]) swap(u, v);
		addrange(rt[vs], rt[vs], 1, n, dfn[top[u]], dfn[u]);
		u = fa[top[u]];
	}
	if(u != v) {
		if(dep[u] > dep[v]) swap(u, v);
		addrange(rt[vs], rt[vs], 1, n, dfn[u]+1, dfn[v]);
	}
}
ll qryroute(int vs, int u, int v) {
	ll ans = 0;
	while(top[u] != top[v]) {
		if(dep[top[u]] < dep[top[v]]) swap(u, v);
		ans += qryrange(rt[vs], 1, n, dfn[top[u]], dfn[u]);
		u = fa[top[u]];
	}
	if(u != v) {
		if(dep[u] > dep[v]) swap(u, v);
		ans += qryrange(rt[vs], 1, n, dfn[u]+1, dfn[v]);
	}
	return ans;
}

ll lans;
ll getval(ll a) {return (lans % (1 << 30ll)) ^ a;}
ll calc(int vs, int v) {
	return suma[vs] + vs * dist[v] - 2 * qryroute(vs, v, 1);
}
int main() {
	n = read(); q = read();
	for(int i = 1; i <= n; i++) a[i] = read();
	for(int i = 1; i < n; i++) {
		int u = read(), v = read(); ll w = read();
		addedge(u, v, w); addedge(v, u, w);
	}
	dfs1(1, 0); dfs2(1, 1);
	for(int i = 1; i <= n; i++) suma[i] = suma[i-1] + dist[a[i]];
	for(int i = 1; i <= n; i++) {
		rt[i] = rt[i-1]; 
		addroute(i, a[i], 1, dist[a[i]]);
	}
	for(int i = 1; i <= q; i++) {
		int t = read();
		if(t == 1) {
			int l = getval(read()), r = getval(read()), v = getval(read());
			printf("%lld\n", lans = calc(r, v) - calc(l-1, v));
		} else {
			int x = getval(read());
			swap(a[x], a[x+1]); suma[x] = suma[x-1] + dist[a[x]];
			if(tot < 2.7e7) {
				rt[x] = rt[x-1]; addroute(x, a[x], 1, dist[a[x]]);
			} else {
				tot = 0;
				for(int i = 1; i <= n; i++) {
					rt[i] = rt[i-1];
					addroute(i, a[i], 1, dist[a[i]]);
				}
			}
		}
	}
	return 0;
}
日居月诸Rijuyuezhu
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
cf757G. Can Bash Save the Day?
miaom的博客
03-02 907
简述题意：给一棵n个节点的数和一个n的数组（n的排列），支持两个操作：1 询问a[l]~a[r]的所有点到某个点的距离，2 交换a[x]与a[x+1]。 一开始想到了这道题bzoj4012: [HNOI2015]开店，大概看了之后发现是开店的加强版。我想到像开店一样使用主席树处理，显然操作2很容易。然而 他的时间是一个log。大概看了题解，感觉套路~   毕竟开店的题解大多是log^2。。
CF757G Can Bash Save the Day?
weixin_30477797的博客
04-03 84
CF757G Can Bash Save the Day? #include<bits/stdc++.h> #define RG register #define IL inline #define _ 200100 #define ycb (1<<30) #define ll long long using namespace std; IL int gi(){ ...
[Codeforces757G]Can Bash Save the Day?——动态点分治(可持久化点分树)
weixin_33849942的博客
12-26 149
题目链接： Codeforces757G 题目大意：给出一棵n个点的树及一个1~n的排列pi，边有边权，有q次操作： 1 l r x 求 $\sum\limits_{i=l}^{r}dis(p_{i},x)$ 2 x $swap(p_{x},p_{x+1})$ $n,q<=2*10^5$，强制在线 如果多次询问一个点到所有点的距离和，我们可以点分树解决，在点分树上每个点x维护点...
【CF757G】Can Bash Save the Day? 可持久化点分树
aodanchui1057的博客
02-25 187
【CF757G】Can Bash Save the Day? 题意：给你一棵n个点的树和一个排列${p_i}$，边有边权。有q个操作： 1 l r x：询问$\sum\limits_{i=l}^r dist(p_i,x)$2 x：$swap(a_x,a_{x+1})$ $n,q\le 2\times 10^5$，强制在线。题解：学到了新姿势：可持久化点分树。我们先将询问变成前缀相减...
【可持久化点分树？？】CF757G Can Bash Save the Day?
MEET YOU FIRST TIME
01-28 558
【题目】 CF 一棵nnn个点的树和一个排列pip_ipi​，边有边权，支持两种操作： l&amp;nbsp;r&amp;nbsp;xl\ r\ xl&amp;nbsp;r&amp;nbsp;x，询问∑i=lrdis(pi,x)\sum\limits _{i=l}^{r} dis(p_i,x)i=l∑r​dis(pi​,x) xxx，交换px,px+1p_x,p_{x+1}px​,px+1​ n,q≤105n,q\leq 1...
算法文档无代码RMQ&LCA问题
04-14
算法文档无代码RMQ&LCA问题提取方式是百度网盘分享地址
LCA (最近公共祖先) Tarjan & 倍增
01-03
LCA Tarjan： 实现原理 理解：离线算法，建好树后再查询，一次DFS 吧所有查询解决完。 时间复杂度：O(n+q)； n个点 q次询问 补一下：链式向前星，并查集 ，Tarjan 代码 #include #include #include #include #...
RMQ & LCA 问题及其相互关系
10-26
一份很不错RMQ与LCA问题的讲解，及其相互关系。推荐。
LCA.rar_LCA
09-19
最近公共祖先(LCA)，转化为 RMQ 用线段树解决
lca.rar_LCA
09-24
设计一个算法，对于给定的树中2 结点返回它们的最近公共祖先。
CF757G Can Bash Save the Day? 题解
莫莫 dicol 的博客
08-31 107
CF757G Can Bash Save the Day? 时间复杂度分析好题（大雾）。 将询问拆成几部分，也就是 dis(u,v)=dep(u)+dep(v)−2∗dep(Lca(u,v))dis(u, v) = dep(u) + dep(v) - 2 * dep(Lca(u, v))dis(u,v)=dep(u)+dep(v)−2∗dep(Lca(u,v))。 显然对于 LcaLcaLca 我们直接进行树链剖分即可。 我们对于每一个 vvv 维护 ≤l\le l≤l 的 dep(Lca)dep(Lc
Codeforces757G Can Bash Save the Day?【可持久化点分树】
Master.Yi的博客
05-11 199
题目描述： 给出nnn个点的树和排列pip_ipi​，有边权，Q次操作： 1 l r x1~l~r~x1 l r x：求∑i=lrdis(pi,x)\sum_{i=l}^rdis(p_i,x)∑i=lr​dis(pi​,x) 2 x2~x2 x：swap(px,px+1)swap(p_x,p_{x+1})swap(px​,px+1​) n,Q≤200000n,Q\le200000n,Q≤200000 题目分析： 首先肯定需要
CF757G Can Bash Save the Day?（可持久化点分树）
Freopen的博客
05-03 290
题目 这个题在CF上跑了15min才跑出来、、 真就可持久化呗。 考虑求一个点到一个点集的距离和。 可以用动态点分治也就是点分树来维护点集加删点和询问操作。 对于这个题，差分一下变成： 多次询问前1...r1...r1...r个组成的点集到一个点的距离。 那么就可以可持久化点分树维护1...r1...r1...r的每个版本的点分树。 对于修改操作，发现其实就是对于xxx这个位置的单点修改，后面的版...
Http请求与响应 contentType
编程日常记录与总结 https://franciswmf.iteye.com
08-23 186
Http请求 1、方法-URI-协议/版本; 2、请求头信息; 3、请求正文; Http响应 1、协议-状态码-描述; 2、响应头信息； 3、响应正文; HttpServletResponse --回的ContentType类型 [url]http://ewf-momo.iteye.com/blog/1902634[/url] google pb: Content-T...
Oracle中先进后出数据结构的例子
codemami的博客
06-19 465
在Oracle数据库中，虽然数据库本身并不直接提供特定于“先进后出”（Last In First Out, LIFO）数据结构的数据类型或功能，但我们可以使用表（table）和过程（procedure）来模拟这种数据结构的行为。上述示例中的SQL语句是基于MySQL的语法，但在Oracle中可能需要进行一些调整，例如使用SEQUENCE来替代AUTO_INCREMENT，以及使用不同的锁定机制来确保并发安全。-- 此时v_value中存储的就是出栈的元素值，可以在这里进行进一步的处理或返回给调用者。
从零开始学数据结构系列之第三章《后续线索二叉树线索化理论》
铁匠匠匠的博客
06-19 838
仔细想想，我们中序遍历这个是 左 右 中，而我们这个先序遍历是直接寻找最左边的节点，所以我们以防万一下图这种情况。我们继续还是用这个图，我们根据前面学习的中序线索二叉树和先序线索二叉树，来看看着后续线索二叉树怎么写。如果是这种情况，我们想上述的代码是不能有效的完成工作的，所以我们需要一些简单的修改。答案是没有的，可以直接正常运行的，那我们接下来看看下面的是否可以直接用。​ 我们再看到中序线索二叉树是怎么写的，看看稍微修改下直接用行不行。就例如创树和结构初始化都是一样的。​ 这是我们中序线索二叉树写的。
【408考点之数据结构】算法和算法评价（时间空间复杂度）
gygkhd的博客
06-18 584
例如，时间复杂度为O(n)表示算法的运行时间随着输入规模n线性增长。常见的时间复杂度还有O(1)、O(log n)、O(n^2)等。例如，空间复杂度为O(1)表示算法只需要常量的额外空间，而空间复杂度为O(n)表示算法需要线性增长的存储空间。在实际应用中，选择算法时需要根据具体问题和资源限制，综合考虑时间和空间复杂度。通过理解算法的基本概念和效率度量，我们可以更好地设计和选择适合具体问题的算法，提高程序的运行效率。函数的空间复杂度为O(1)，因为只使用了一个额外的临时变量。，与输入数组的大小无关。
手撕排序2--选择排序(直接选择+堆排序
最新发布
2302_81153999的博客
06-21 264
把最大的换到最后，不把其看做是堆里的。前n-1个数向下调整，选出次大的数，再跟倒数第二个位置交换........每一次从待排序的数据元素中选出最小（或最大）的一个元素，存放在序列的起始位置，直到全部待排序的数据元素排完。堆排序是利用堆 这种数据结构而设计的一种排序算法，它是选择排序的一种。将该节点与左右子树的节点进行比较，重新构建成一个大堆或者小堆。每个结点的值都小于或等于其左右孩子结点的值，称为。每个结点的值都大于或等于其左右孩子结点的值。物理结构是：一个数组。1. 对待排序的数组，先构建成 堆。
LCA features enabled by '-lca' argument on the command line
03-31
However, based on my understanding, the '-lca' argument on the command line typically enables the use of lowest common ancestor (LCA) features in a given program or application. LCA is a term used ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

日居月诸Rijuyuezhu

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

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

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

打赏作者

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

抵扣说明：

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

余额充值