树上点差分的经典应用 LuoguP3258松鼠的新家
树上点差分的核心就是如何避免重复,即正确的运用差分数组
例如a,b点路径上点权值加1,则把a,b路径找到,并找到其LCA,此时可以把a到根,b到根这两条路径看出两条链,把每条链看出我们熟悉的
顺序差分结构.以其中一条链为例子,把a当成数组的起点,根当成数组的末尾,进行差分,显然有C[a]++,C[f[lca][0]]--(这里f[lca][0]为lca的父节点)
附上图片理解(摘自Oi-wiki)

故对a到b的路径上的差分修改的完整过程为C[a]++,C[b]++,C[lca]--(因为加了两次,要减去1次),C[f[lca][0]]--;
最后附上这条代码(里面还有一些具体应用的细节)
点击查看代码
//本题有点特殊,因为同时作为路径起点与终点的点差分值多算了一次,但仅限于点差分!!!
#include <bits/stdc++.h>
using namespace std;
typedef struct edge//链式前向星
{
int to,next;
}EDGE;
EDGE e[600005];//边
int head[300001];//以i为起点的第一条边
int cnt;//边的数量
int number[300001];//答案,同时兼具差分功能
int a[300001];//访问顺序
int f[300001][30];//倍增父节点
int depth[300001];//深度,服务于倍增
inline int read()//快读
{
int s = 0;
char ch = getchar();
while(ch > '9'||ch < '0') ch = getchar();
while(ch >= '0'&&ch <= '9')
{
s = s*10 + ch - '0';
ch = getchar();
}
return s;
}
void add(int x,int y)//加边
{
e[++cnt].to = y;
e[cnt].next = head[x];
head[x] = cnt;
}
void DFS(int x,int fa)//建树及预处理LCA
{
depth[x] = depth[fa] + 1;//深度加1
f[x][0] = fa;//第一个父节点
for (int i = 1;(1<<i)<=depth[x];i++)//倍增父节点
{
f[x][i] = f[f[x][i-1]][i-1];
}
for (int i = head[x];i;i = e[i].next)//遍历子树
{
if(e[i].to == fa) continue;
DFS(e[i].to,x);
}
}
void LCA(int a,int b)//求LCA兼具差分
{
if(depth[a]<depth[b]) swap(a,b);//不妨a的深度大于b
number[a]++;//差分思想
number[b]++;//差分思想
while(depth[a]>depth[b])//向上跳至深度相同
{
a = f[a][(int)(log(depth[a]-depth[b])/log(2))];//以后尽量不用用log写法,容易出错且不美观
}
if(a == b)
{
number[a]--;//祖先减1
number[f[a][0]]--;//祖先的父亲减1,这样保证了LCA只减了1
return;一定要返回
}
for (int i = (int)(log(depth[a])/log(2));i>=0;i--)//以后尽量不用用log写法,容易出错且不美观
{
if(f[a][i] != f[b][i])
{
a = f[a][i];
b = f[b][i];
}
}
number[f[a][0]]--;//同上原理
number[f[f[a][0]][0]]--;
}
void Solve(int x,int fa)//计算子树和
{
for (int i = head[x];i;i = e[i].next)
{
if(e[i].to == fa) continue;
Solve(e[i].to,x);
number[x] += number[e[i].to];//及求区间和,差分思想
}
}
int main()
{
int n = read();//读入节点数
for (int i = 1;i<=n;i++)//读入访问顺序
{
a[i] = read();
}
for (int i = 1;i<=n-1;i++)//建树
{
int x = read();
int y = read();
add(x,y);
add(y,x);
}
depth[0] = -1;
DFS(1,0);//预处理
for (int i = 1;i<=n-1;i++)//求LCA及差分
{
LCA(a[i],a[i+1]);
}
Solve(1,0);//计算子树和,计算区间和(注意这里是先计算区间和,再对多加点进行减法操作)
for (int i = 2;i<=n;i++) number[a[i]]--;//同时为起点终点的多加了一次,要减1,同时最后一点为餐厅,不用糖果,减1
for (int i = 1;i<=n;i++)//输出
{
printf("%d\n",number[i]);
}
return 0;
}
码字不易,多多支持!!
树上点差分的经典应用 LuoguP3258松鼠的新家的更多相关文章
- luoguP3258 [JLOI2014]松鼠的新家 题解(树上差分)
P3258 [JLOI2014]松鼠的新家 题目 树上差分:树上差分总结 #include<iostream> #include<cstdlib> #include<c ...
- [BZOJ3631]:[JLOI2014]松鼠的新家(LCA+树上差分)
题目传送门 题目描述: 松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的.天哪,他居然真的住在“树”上.松鼠想邀 ...
- [填坑]树上差分 例题:[JLOI2014]松鼠的新家(LCA)
今天算是把LCA这个坑填上了一点点,又复习(其实是预习)了一下树上差分.其实普通的差分我还是会的,树上的嘛,也是懂原理的就是没怎么打过. 我们先来把树上差分能做到的看一下: 1.找所有路径公共覆盖的边 ...
- 【洛谷】【lca+树上差分】P3258 [JLOI2014]松鼠的新家
[题目描述:] 松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n(2 ≤ n ≤ 300000)个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的.天哪,他居然真 ...
- [JLOI2014] 松鼠的新家 (lca/树上差分)
[JLOI2014]松鼠的新家 题目描述 松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的.天哪,他居然真的住在 ...
- 洛谷P3258 [JLOI2014]松鼠的新家(树上差分+树剖)
题目描述 松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的.天哪,他居然真的住在”树“上. 松鼠想邀请小熊维尼前 ...
- BZOJ3631 [JLOI2014]松鼠的新家 【树上差分】
题目 松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的.天哪,他居然真的住在"树"上.松鼠想 ...
- P3258[JLOI2014]松鼠的新家(LCA 树上差分)
P3258 [JLOI2014]松鼠的新家 题目描述 松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的.天哪,他 ...
- BZOJ 3631: [JLOI2014]松鼠的新家 树上差分 + LCA
Description 松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的.天哪,他居然真的住在“树”上.松鼠想邀 ...
- 【bzoj3631】[JLOI2014]松鼠的新家 LCA+差分数组
题目描述 松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的.天哪,他居然真的住在“树”上.松鼠想邀请小熊维尼前来 ...
随机推荐
- MindSpore自定义算子中的张量维度问题
技术背景 在前面的几篇博客中,我们介绍了MindSpore框架下使用CUDA来定义本地算子的基本方法,以及配合反向传播函数的使用,这里主要探讨一下MindSpore框架对于CUDA本地算子的输入输出的 ...
- [VueJsDev] 日志 - BBTime-LOG
[VueJsDev] 目录列表 https://www.cnblogs.com/pengchenggang/p/17037320.html BBTime-LOG ::: details 目录 目录 B ...
- 被 两个 vetur 玩惨了的几天,可不是vue-helper的锅
最近由idea 转 vscode,因为 由原来一个项目变成7个项目了,2G内存不够看. 然后发现点击组件跳转不好使,装个vue-helper也是有一部分不好使. 下载vue-helper源码一通查,发 ...
- python3中print()函数打印多个变量值
第一种方法: print("变量1", file_name, "变量2", new_name) print("变量1", file_name ...
- GDB调试之多线程
1.set scheduler-locking off/on/step 调试时除了当前线程在运行,要想规定其他线程的运行情况用这个命令,有三个选择: set scheduler-locking off ...
- Android native项目编译v7a时-Werror参数导致报错
在GitHub上下载了一个demo,自己编译运行时,v8a的so是可以正常编译运行的,但是就是v7a的不能也一直报错,最后发现在Android.mk文件里去掉-Werror字段就可以了.于是上网了解了 ...
- 广西首次!3DCAT实时云渲染助力南宁数字气象科普馆上线
11月17日,南宁市气象局举行数字气象科普馆上线发布仪式.目前,以气象为主题的三维虚拟数字气象科普馆在广西尚属首创.17日起,在微信搜索小程序南宁气象科普馆&就能在3DCAT实时云渲染技术助力 ...
- npm install安装依赖包时报错npm ERR! command C:\Windows\system32\cmd.exe /d /s /c node install.js,npm ERR! ChromeDriver installation failed Error with http(s) request: Error: read ECONNRESET
PS E:\20231213\uirecorder> PS E:\20231213\uirecorder> PS E:\20231213\uirecorder> PS E:\2023 ...
- Scala 简单分词求和
1 package chapter07 2 3 object Test17_CommonWordCount { 4 def main(args: Array[String]): Unit = { 5 ...
- 性能对比 Go、Python、Perl、Ruby、Rust、C/C++、PHP、Node.js、Java.. 等多编
1. 有人说 Python 性能没那么 Low? 这个我用 pypy 2.7 确认了下,确实没那么差, 如果用 NumPy 或 其它版本 Python 的话,性能更快.但 pypy 还不完善,pypy ...