【NOIP】提高组2016 天天爱跑步
【题意】n个点的树,有m个人同时开始走链,每一步花一秒,n个点都有观察员在ai秒观察,求每个观察员观察到的人数。
【算法】树上差分(主席树||线段树合并)
【题解】一个人的走链可以拆成u-lca和lca-v两部分,可以发现在u-lca链上的点能观察到这个人的w[x],满足所有deep[x]+w[x]相等。同理,在lca-v链上需满足deep[x]-w[x]相等。(画个图,将深度和秒数都标出来就可以得到结论了)
所以可以得到一个点观察到链u-v的条件是【在u-lca上且deep[x]+w[x]=deep[u]】或【在lca-v上且deep[x]-w[x]=deep[v]-time】,其中time=deep[u]+deep[v]-2*deep[lca]。
为了判断是否满足条件,我们可以维护数组A(针对deep[x]+w[x])和数组B(针对deep[x]-w[x])作为判断数组。对于每条链A[deep[u]]+1和B[deep[v]-time]+1,然后点x就可以直接判断满足多少链的要求。
接下来需要解决是否在u-lca或lca-v上的问题,可以运用树上差分。对于一条链,在u点标记A[deep[u]]+1的操作,在v点标记B[deep[v]-time]+1的操作,在lca标记A[deep[u]]-1的操作,在fa[lca]标记B[deep[v]-time]-1的操作(否则lca点同时满足两个要求会计算两次),标记操作需要用类似邻接表的链表接在每个节点处。
然后总的进行一次dfs,每个节点按顺序进行:统计ans[x],将标记操作,dfs子节点,ans[x]=统计ans[x]-ans[x](原)。
★注意:
1.树上差分只能是u+1,v+1,lca-2(或lca-1&&fa(lca)-1),绝对不能是lca处加,因为lca处加会影响到所有子树而不止一条链。
2.树上差分的操作顺序必须是:【统计ans】【操作标记】【DFS子节点】【统计ans-原ans】
对点x,最后统计ans是x的子树和x上面已统计的部分。减原ans相当于减掉x上面已统计的部分。
为什么必须这样做?要是在返回的时候操作标记,就会干扰父节点另一棵子树,所以必须用 [ 当前已统计的所有-x上方已统计的所有 ] 这种后统计-先统计的强操作。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
using namespace std;
int read(){
char c;int s=,t=;
while(!isdigit(c=getchar()))if(c=='-')t=-;
do{s=s*+c-'';}while(isdigit(c=getchar()));
return s*t;
}
const int maxn=;
int first[maxn],firstA[maxn],firstB[maxn],a[maxn],markA[maxn*],markB[maxn*],tot,totA,totB;
int deep[maxn],f[maxn][],ans[maxn],n,m;
struct edge{int v,w,from;}e[maxn*],A[maxn*],B[maxn*];
void insert(int u,int v){tot++;e[tot].v=v;e[tot].from=first[u];first[u]=tot;}
void insertA(int u,int v,int w){totA++;A[totA].v=v;A[totA].w=w;A[totA].from=firstA[u];firstA[u]=totA;}
void insertB(int u,int v,int w){totB++;B[totB].v=v;B[totB].w=w;B[totB].from=firstB[u];firstB[u]=totB;}
void DFS(int x,int fa){
for(int j=;(<<j)<=deep[x];j++)f[x][j]=f[f[x][j-]][j-];
for(int i=first[x];i;i=e[i].from)if(e[i].v!=fa){
deep[e[i].v]=deep[x]+;
f[e[i].v][]=x;
DFS(e[i].v,x);
}
}
int lca(int x,int y){
if(deep[x]<deep[y])swap(x,y);
int d=deep[x]-deep[y];
for(int j=;j<=;j++)if((<<j)&d)x=f[x][j];
if(x==y)return x;
for(int i=;i>=;i--)if((<<i)<=deep[x]&&f[x][i]!=f[y][i]){
x=f[x][i];y=f[y][i];
}
return f[x][];
}
void dfs(int x,int fa){
ans[x]=markA[a[x]+deep[x]]+markB[a[x]-deep[x]+n];//
for(int i=firstA[x];i;i=A[i].from){
if(A[i].w)markA[A[i].v]--;else markA[A[i].v]++;
}
for(int i=firstB[x];i;i=B[i].from){
if(B[i].w)markB[B[i].v]--;else markB[B[i].v]++;
}
for(int i=first[x];i;i=e[i].from)if(e[i].v!=fa)dfs(e[i].v,x);
ans[x]=markA[a[x]+deep[x]]+markB[a[x]-deep[x]+n]-ans[x];
}
int main(){
n=read();m=read();
for(int i=;i<n;i++){
int u=read(),v=read();
insert(u,v);insert(v,u);
}
for(int i=;i<=n;i++)a[i]=read();
DFS(,);f[][]=;
for(int i=;i<=m;i++){
int u=read(),v=read();
int w=lca(u,v);
insertA(u,deep[u],);//0plus 1minus
insertA(f[w][],deep[u],);
insertB(v,deep[u]-*deep[w]+n,);
insertB(w,deep[u]-*deep[w]+n,);
}
for(int i=firstA[];i;i=A[i].from)if(A[i].w)markA[A[i].v]--;else markA[A[i].v]++;
for(int i=firstB[];i;i=B[i].from)if(A[i].w)markB[B[i].v]--;else markB[B[i].v]++;
dfs(,);
for(int i=;i<n;i++)printf("%d ",ans[i]);printf("%d",ans[n]);
return ;
}
顺便一提主席树和线段树合并的做法。
主席树:对每个点询问子树内等于某个值的起点数和等于某个值的终点树,转化为dfs序上的区间权值询问,可以用可持久化权值线段树解决。
线段树合并:线段树下标记录关键值(和差分的两个值一样),然后从叶子往上进行线段树合并和查询,非常套路了。
【NOIP】提高组2016 天天爱跑步的更多相关文章
- NOIP提高组2016 D1T2 【天天爱跑步】
码了一个下午加一个晚上吧...... 题目描述: 小c同学认为跑步非常有趣,于是决定制作一款叫做<天天爱跑步>的游戏.<天天爱跑步>是一个养成类游戏,需要玩家每天按时上线,完成 ...
- BZOJ4719[NOIP2016提高组Day1T2] 天天爱跑步
#261. [NOIP2016]天天爱跑步 描述 提交 自定义测试 小C同学认为跑步非常有趣,于是决定制作一款叫做<天天爱跑步>的游戏.<天天爱跑步>是一个养成类游戏,需要玩家 ...
- NOIP2016提高组D1T2 天天爱跑步
n<=300000个点的树,每个点有个人于第Ti秒观测,有m<=300000个人于时间0开始从Sj跑到Tj,速度1个点每秒,输出每个点上的人观察到的跑步的人的数量. 前25分:直接模拟每条 ...
- NOIP2016提高组Day1T2 天天爱跑步 树链剖分 LCA 倍增 差分
原文链接https://www.cnblogs.com/zhouzhendong/p/9275606.html 题目传送门 - 洛谷P1600 题目传送门 - LOJ#2359 题目传送门 - Vij ...
- 题解——洛谷P2827 NOIP提高组 2016 蚯蚓
队列模拟 详细题解待填坑 #include <cstdio> #include <algorithm> #include <queue> #include < ...
- NOIP提高组2016总结
前言 大翻车! 300--: day1 8:30~9:00, 照常看题,思考. 9:00~9:15, 搞定第一题,很水. 9:15~9:45, 思考第二题,我考虑用分深度来处理,想出个个玄学暴力,但刚 ...
- NOIP提高组2016 D2T3 【愤怒的小鸟】
貌似还没有写过状压DP的题目,嗯,刚好今天考了,就拿出来写一写吧. 题目大意: 额,比较懒,这次就不写了... 思路分析: 先教大家一种判断题目是不是状压DP的方法吧. 很简单,那就是--看数据范围! ...
- [NOIP]2016天天爱跑步
[NOIP]2016天天爱跑步 标签: LCA 树上差分 NOIP Description 小C同学认为跑步非常有趣,于是决定制作一款叫做<天天爱跑步>的游戏.<天天爱跑步>是 ...
- NOIP提高组初赛难题总结
NOIP提高组初赛难题总结 注:笔者开始写本文章时noip初赛新题型还未公布,故会含有一些比较老的内容,敬请谅解. 约定: 若无特殊说明,本文中未知数均为整数 [表达式] 表示:在表达式成立时它的值为 ...
随机推荐
- Zigbee安全基础篇Part.1
原文地址: https://www.4hou.com/wireless/14211.html 导语:ZigBee是一种开源无线技术,用于低功耗嵌入式设备(无线电系统).本文探讨了ZigBee协议的可用 ...
- Bookmark Sentry – 检查重复、删除死链书签 Chrome扩展
Bookmark Sentry 的用处,就是 处理重复的收藏夹的死链 . 重复链收藏.具体,请百度. Bookmark Sentry 下载 : https://files.cnblogs.com/f ...
- Kafka性能之道
Kafka高性能之道 高效使用磁盘 零拷贝 批处理和压缩 Partition ISR 高效使用磁盘 >顺序写cipan >Append Only(数据不更新,无记录级的数据删除,只会整个s ...
- 【Python】python基础_代码编写注意事项
1. 说明使用的编译方式 1 #!/usr/bin/python 2. 说明字符编码方式 1 #coding=utf-8 3. print 默认输出是换行的,如果要实现不换行需要在变量末尾加上逗号 # ...
- Nginx + Keepalived使用文档
第一步: 下载keepalived地址:http://www.keepalived.org/download.html 解压安装: tar -zxvf keepalived-1.2.18.tar.gz ...
- 关联容器 // append方法
关联容器和顺序容器的差别在于:关联容器通过键(key)存储和读取元素,而顺序容器则通过元素在容器中的位置顺序存储和访问元素. 1.关联容器支持通过键来高效地查找和读取元素.两个基本的关联容器类型是ma ...
- 【bzoj3732】Network 最小生成树+倍增LCA
题目描述 给你N个点的无向图 (1 <= N <= 15,000),记为:1…N. 图中有M条边 (1 <= M <= 30,000) ,第j条边的长度为: d_j ( 1 & ...
- poj 1719 Shooting Contest (二分匹配)
Shooting Contest Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 3812 Accepted: 1389 ...
- JavaScript中:表达式和语句的区别
JavaScript中:表达式和语句的区别 Javascript语言精粹:表达式是由运算符构成,并运算产生结果的语法结构.程序是由语句构成,语句则是由“:(分号)”分隔的句子或命令.如果在表达式后面加 ...
- POJ2142:The Balance——题解
http://poj.org/problem?id=2142 题目大意:有一天平和两种数量无限的砝码(重为a和b),天平左右都可以放砝码,称质量为c的物品,要求:放置的砝码数量尽量少:当砝码数量相同时 ...