暴力移步 http://www.cnblogs.com/TheRoadToTheGold/p/6673430.html

首先解决本题应用的知识点:

dfs序——将求子树的信息(树形)转化为求一段连续区间信息(线形)

线段树——求区间信息

树上差分——统计答案

lca——拆分路径

树链剖分——求lca

另deep[]表示节点的深度,watch[]表示观察者的出现时间,s表示玩家起点,t表示终点

固定节点的观察者出现的时间固定,说明对这个观察者有贡献的点是有限且固定的

只有满足  观察者出现时间=玩家起点与观察者距离的  玩家才对观察者有贡献

每条路径拆成   起点到lca(向上跑)  和   终点到lca的子节点(向下跑)  的两条路径

对于向上跑的,如果玩家能被观察员i观察到,那么deep[s]-deep[i]=watch[i]   式①

对于向下跑的,就是 deep[s]+deep[i]-2*deep[lca(s,i)]=watch[i]  式②

等号左边是玩家起点与观察者的距离,等号右边是观察者出现时间

向上跑的很显然,向下跑的如何理解?

假设我们知道点a,b到lca(a,b)的距离分别为da,db,那么a,b之间的距离=da+db

但这里的deep不是到lca的距离,是深度,即到根节点的距离+1

deep[s]+deep[i]包含2段信息,1、s到i的距离,   2、lca(s,i)到根节点的距离+1

第2段包含了2次,所以减去

先看向上跑的

玩家路径:玩家起点 到 起点与终点的lca

将式①移项,deep[s]=deep[i]+watch[i]

发现等号右边是定值

也就是说对与观察者i,他所能观察到的向上跑的玩家,是所有的起点=deep[i]+watch[i]的玩家

换句话说,以i为根的子树中,所有深度为deep[i]+watch[i]的玩家都能被i观察到

我们如果搞一个dfs序,i的在a时入栈,在b时出栈,

那么以i为根的子树就可以转化为区间[a,b]

深度咋整?

我们对每个深度建立一颗线段树(动态加点)

那么问题就转化为了  在深度为deep[i]+watch[i]的线段树中,查询区间[a,b]的玩家个数

现在就差玩家个数了

很容易想到在起点处+1

但是还要在起点与终点的lca的父节点处-1

差分惯用思想

用sum[]统计这些1和-1的和

那么问题就转化为了  在深度为deep[i]+watch[i]的线段树中,查询区间[a,b]的sum和

提问:为什么是起点处+1,lca的父节点处-1,可以反过来吗?

不可以。

因为起点的深度深,lca的父节点深度浅,在深度深的节点处+1,以深度深度浅的点为根的子树可以包含这个点

想想序列上的差分,是左端点+1,右端点后面的点-1

因为序列差分与前缀和相联系,前面的点的信息对后面的点会产生影响,所以只需加一个1

这里查询的是子树信息,是这个点深度及以下的信息

对照理解即可

向下跑的同理,只简单说怎么做

玩家路径:lca的子节点到玩家终点

把式②移项 deep[s]-2*deep[lca(s,i)]=watch[i]-deep[i]

在watch[i]-deep[i]深度为deep[s]-2*deep[lca(s,i)]的线段树中,终点处+1,lca处-1

查询时查深度为watch[i]-deep[i]的线段树即可

2个小问题:

1、做完向上跑的后,不要忘了清空线段树

2、向下跑的deep[s]-2*deep[lca(s,i)]可能会产生负数,所以全体后移一定长度,root[]数组开大

我后移了2*n,那么root[]数组要开3倍

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 300401
using namespace std;
int n,m,fa[N],son[N],deep[N],bl[N],sz,id[N],ans[N];
int in[N],out[N],watch[N];
int front[N],nextt[N*],to[N*];
int root[N*],lc[N*],rc[N*],sum[N*],tot,cnt;
struct node
{
int s,t,lca;
}runner[N];
void add(int u,int v)
{
to[++cnt]=v; nextt[cnt]=front[u]; front[u]=cnt;
}
void dfs1(int now)
{
son[now]++;
for(int i=front[now];i;i=nextt[i])
{
if(to[i]==fa[now]) continue;
deep[to[i]]=deep[now]+;
fa[to[i]]=now;
dfs1(to[i]);
son[now]+=son[to[i]];
}
}
void dfs2(int now,int chain)
{
id[now]=++sz;
in[now]=sz;
bl[now]=chain;
int y=;
for(int i=front[now];i;i=nextt[i])
{
if(to[i]==fa[now]) continue;
if(son[to[i]]>son[y]) y=to[i];
}
if(!y)
{
out[now]=sz;
return;
}
dfs2(y,chain);
for(int i=front[now];i;i=nextt[i])
{
if(to[i]==fa[now]||to[i]==y) continue;
dfs2(to[i],to[i]);
}
out[now]=sz;
}
int getlca(int u,int v)
{
while(bl[u]!=bl[v])
{
if(deep[bl[u]]<deep[bl[v]]) swap(u,v);
u=fa[bl[u]];
}
return deep[u]<deep[v] ? u:v;
}
void change(int & now,int l,int r,int pos,int w)
{
if(!pos) return;
if(!now) now=++tot;
sum[now]+=w;
if(l==r) return;
int mid=l+r>>;
if(pos<=mid) change(lc[now],l,mid,pos,w);
else change(rc[now],mid+,r,pos,w);
}
int query(int now,int l,int r,int opl,int opr)
{
if(!now) return ;
if(l==opl&&r==opr) return sum[now];
int mid=l+r>>;
if(opr<=mid) return query(lc[now],l,mid,opl,opr);
else if(opl>mid) return query(rc[now],mid+,r,opl,opr);
else return query(lc[now],l,mid,opl,mid)+query(rc[now],mid+,r,mid+,opr);
}
void clear()
{
tot=;
memset(lc,,sizeof(lc));
memset(rc,,sizeof(rc));
memset(sum,,sizeof(sum));
memset(root,,sizeof(root));
}
int main()
{
/*freopen("runninga.in","r",stdin);
freopen("runninga.out","w",stdout);*/
scanf("%d%d",&n,&m);
int u,v;
for(int i=;i<n;i++)
{
scanf("%d%d",&u,&v);
add(u,v);add(v,u);
}
for(int i=;i<=n;i++) scanf("%d",&watch[i]);
for(int i=;i<=m;i++) scanf("%d%d",&runner[i].s,&runner[i].t);
dfs1();
dfs2(,);
for(int i=;i<=m;i++) runner[i].lca=getlca(runner[i].s,runner[i].t);
int now;
for(int i=;i<=m;i++)
{
now=deep[runner[i].s];
change(root[now],,n,id[runner[i].s],);
change(root[now],,n,id[fa[runner[i].lca]],-);
}
for(int i=;i<=n;i++) ans[i]=query(root[deep[i]+watch[i]],,n,in[i],out[i]);
clear();
for(int i=;i<=m;i++)
{
now=deep[runner[i].s]-deep[runner[i].lca]*+n*;
change(root[now],,n,id[runner[i].t],);
change(root[now],,n,id[runner[i].lca],-);
}
for(int i=;i<=n;i++) ans[i]+=query(root[watch[i]-deep[i]+n*],,n,in[i],out[i]);
for(int i=;i<=n;i++) printf("%d ",ans[i]);
}

NOIP2016 天天爱跑步 正解的更多相关文章

  1. [Noip2016]天天爱跑步 LCA+DFS

    [Noip2016]天天爱跑步 Description 小c同学认为跑步非常有趣,于是决定制作一款叫做<天天爱跑步>的游戏.?天天爱跑步?是一个养成类游戏,需要玩家每天按时上线,完成打卡任 ...

  2. [NOIp2016]天天爱跑步 线段树合并

    [NOIp2016]天天爱跑步 LG传送门 作为一道被毒瘤出题人们玩坏了的NOIp经典题,我们先不看毒瘤的"动态爱跑步"和"天天爱仙人掌",回归一下本来的味道. ...

  3. 【LG1600】[NOIP2016]天天爱跑步

    [LG1600][NOIP2016]天天爱跑步 题面 洛谷 题解 考虑一条路径\(S\rightarrow T\)是如何给一个观测点\(x\)造成贡献的, 一种是从\(x\)的子树内出来,另外一种是从 ...

  4. BZOJ4719 [Noip2016]天天爱跑步

    本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000作者博客:http://www.cnblogs.com/ljh2000-jump/转 ...

  5. bzoj 4719: [Noip2016]天天爱跑步

    Description 小c同学认为跑步非常有趣,于是决定制作一款叫做<天天爱跑步>的游戏.?天天爱跑步?是一个养成类游戏,需要 玩家每天按时上线,完成打卡任务.这个游戏的地图可以看作一一 ...

  6. NOIP2016 天天爱跑步 80分暴力

    https://www.luogu.org/problem/show?pid=1600 题目描述 小c同学认为跑步非常有趣,于是决定制作一款叫做<天天爱跑步>的游戏.«天天爱跑步»是一个养 ...

  7. NOIP2016天天爱跑步 题解报告【lca+树上统计(桶)】

    题目描述 小c同学认为跑步非常有趣,于是决定制作一款叫做<天天爱跑步>的游戏.«天天爱跑步»是一个养成类游戏,需要玩家每天按时上线,完成打卡任务. 这个游戏的地图可以看作一一棵包含 nn个 ...

  8. ☆ [NOIp2016] 天天爱跑步 「树上差分」

    题目类型:LCA+思维 传送门:>Here< 题意:给出一棵树,有\(M\)个人在这棵树上跑步.每个人都从自己的起点\(s[i]\)跑到终点\(t[i]\),跑过一条边的时间为1秒.现在每 ...

  9. 4719: [Noip2016]天天爱跑步

    Time Limit: 40 Sec Memory Limit: 512 MB Submit: 1986 Solved: 752 [Submit][Status][Discuss] Descripti ...

随机推荐

  1. SVN的安装以及和eclipse的结合使用

    SVN概述 l 通常软件开发由多人协作开发,如果对代码文件.配置文件.文档等没有进行版本控制,将会出现很多问题: l 备份多个版本,占用磁盘空间大 l 解决代码冲突困难 l 容易引发BUG l 难于追 ...

  2. [THUWC 2017]在美妙的数学王国中畅游

    bzoj5020 \[答案误差只要小于 10^{-7}\] 题解 Taylor展开式: \[若f(x)的n阶导数在[a, b]内连续,则f(x)在x_{0}\in[a, b]可表示为\] \[f(x) ...

  3. WPF自学入门(五)WPF依赖属性

    在.NET中有事件也有属性,WPF中加入了路由事件,也加入了依赖属性.最近在写项目时还不知道WPF依赖属性是干什么用的,在使用依赖项属性的时候我都以为是在用.NET中的属性,但是确实上不是的,通过阅读 ...

  4. 1.4 如何在main()方法之前执行输出“hello world”

    public class Test{ static{ System.out.println("hello world"); } public static void main(St ...

  5. 关系型数据库工作原理-客户端连接管理器(翻译自Coding-Geek文章)

    本文翻译自Coding-Geek文章:< How does a relational database work>.原文链接:http://coding-geek.com/how-data ...

  6. ubuntu中eclipse 不支持gbk编码问题解决办法

    今天在ubuntu 下, 把Windows下工程导入Linux下Eclipse中,由于工程代码,是GBK编码,而Ubuntu默认不支持GBK编码,所以,要让Ubuntu支持GBK. 方法如下: 1.修 ...

  7. python3.5连接oracle数据及数据查询

    今天心血来潮研究下用python连接oracle数据库,看了一下demo,本以为很简单,从操作到成功还是有点坎坷,这里分享给大家,希望为后面学习的童鞋铺路. 一.首先按照cx_Oracle 二:在py ...

  8. 记kkpager分页控件的使用

    kkpager支持异步加载分页: 1.页面添加div标签和引用JS,默认标签为<div id="kkpager"></div> 引用JS和样式 <sc ...

  9. Jquery入门(初学者易懂)

    一.定义 jQuery是一个快速.简洁的JavaScript框架,是继Prototype之后又一个优秀的JavaScript代码库(或JavaScript框架).jQuery设计的宗旨是"w ...

  10. Mecanim之IK动画

    序言:IK动画全名是Inverse Kinematics 意思是逆向动力学,就是子骨骼节点带动父骨骼节点运动. 比如体操运动员,只靠手来带动身体各个部位的移动.手就是子骨骼,身体就是它的父骨骼,这时运 ...