传送门

题目分析:

一年前还是个傻子的时候居然直接放弃了这题。

首先列出两个方程:如果i节点的观察员能够观察到由s->t的那个人,那么:

\[dep[s] - dep[i] = w[i], dep[t] - dep[i] = len - w[i]
\]

整理得到:$$dep[s] = w[i] + dep[i] (1), dep[t] - len = dep[i] - w[i] (2)$$

也就是说:只要起点满足方程(1)和终点满足方程(2)的都能被i看到,下面考虑下差分:因为树的dfs是连续的,所以以节点i为根节点的子树的答案就是:扫描子树后的 - 扫描子树前的,对于一个人从s->t,在s节点+dep[s],在t节点+dep[t] - len,在lca节点-dep[s], 在fa[lca]节点-(dep[t] - len)(s+, t+, lca-, fa[lca]-,这样才能做到补充不漏,自行脑补),每次处理该点的差分标记更新cnt1和cnt2数组(cnt1表示关于起点的差分值的个数,cnt2表示关于终点的差分值的个数,对应上面的方程),节点i的答案就是:$$(扫描后的cnt1[w[i] + dep[i]] + 扫描后的cnt2[dep[i] - w[i]]) - (扫描前的cnt1[w[i] + dep[i]] + 扫描前的cnt2[dep[i] - w[i]])$$

最后还要注意cnt数组记录的差分值可能出现负数,所以要给作差的差分值+offset,来保证下标是合法的。

code

#include<bits/stdc++.h>
using namespace std; typedef long long ll; namespace IO{
inline ll read(){
ll i = 0, f = 1; char ch = getchar();
for(; (ch < '0' || ch > '9') && ch != '-'; ch = getchar());
if(ch == '-') f = -1, ch = getchar();
for(; ch >= '0' && ch <= '9'; ch = getchar()) i = (i << 3) + (i << 1) + (ch - '0');
return i * f;
}
inline void wr(ll x){
if(x < 0) putchar('-'), x = -x;
if(x > 9) wr(x / 10);
putchar(x % 10 + '0');
}
}using namespace IO; const int N = 3e5 + 5, M = 3e5 + 5;
int n, m, ecnt, adj[N], nxt[M << 1], go[M << 1];
int cnt1[N * 3], cnt2[N * 3], w[N], offset;
int fa[N][25], dep[N], ans[N];
struct node{int type, val, flag;};
vector<node> tag[N]; inline void addEdge(int u, int v){nxt[++ecnt] = adj[u], adj[u] = ecnt, go[ecnt] = v;} inline void pre(int u, int f){
dep[u] = dep[f] + 1;
fa[u][0] = f;
for(int i = 1; i <= 20; i++) fa[u][i] = fa[fa[u][i - 1]][i - 1];
for(int v, e = adj[u]; e; e = nxt[e]){
if((v = go[e]) == f) continue;
pre(v, u);
}
} inline int getLca(int u, int v){
if(dep[u] < dep[v]) swap(u, v);
int delta = dep[u] - dep[v];
for(int i = 20; i >= 0; i--) if((1 << i) & delta) u = fa[u][i];
if(u == v) return u;
for(int i = 20; i >= 0; i--) if(fa[u][i] != fa[v][i]) u = fa[u][i], v = fa[v][i];
return fa[u][0];
} inline void addTag(int u, int v){
int lca = getLca(u, v), len = dep[u] + dep[v] - 2 * dep[lca];
tag[u].push_back((node){1, dep[u], 1}), tag[v].push_back((node){2, dep[v] - len + offset, 1});
tag[lca].push_back((node){1, dep[u], -1}), tag[fa[lca][0]].push_back((node){2, dep[v] - len + offset, -1});
} inline void dfs(int u, int f){
int cur = cnt1[w[u] + dep[u]] + cnt2[dep[u] - w[u] + offset];
for(int i = 0; i < tag[u].size(); i++){
if(tag[u][i].type == 1) cnt1[tag[u][i].val] += tag[u][i].flag;
else cnt2[tag[u][i].val] += tag[u][i].flag;
}
for(int v, e = adj[u]; e; e = nxt[e]){
if((v = go[e]) == f) continue;
dfs(v, u);
}
ans[u] = -cur + cnt1[w[u] + dep[u]] + cnt2[dep[u] - w[u] + offset];
} int main(){
freopen("h.in", "r", stdin);
n = read(), m = read(); offset = n * 3;
for(int i = 1; i < n; i++){int x = read(), y = read(); addEdge(x, y), addEdge(y, x);}
pre(1, 0);
for(int i = 1; i <= n; i++) w[i] = read();
for(int i = 1; i <= m; i++){int u = read(), v = read(); addTag(u, v);}
dfs(1, 0);
for(int i = 1; i <= n; i++) wr(ans[i]), putchar(' ');
return 0;
}

NOIP2016 天天爱跑步 - 树上差分的更多相关文章

  1. 洛谷 1600 (NOIp2016) 天天爱跑步——树上差分

    题目:https://www.luogu.org/problemnew/show/P1600 看TJ:https://blog.csdn.net/clove_unique/article/detail ...

  2. NOIP2016 天天爱跑步 (树上差分+dfs)

    题目大意:给你一颗树,树上每个点都有一个观察员,他们仅会在 w[i] 时刻出现,观察正在跑步的玩家 一共有m个玩家,他们分别从节点 s[i] 同时出发,以每秒跑一条边的速度,沿着到 t[i] 的唯一路 ...

  3. [NOIP2016]天天爱跑步(树上差分+线段树合并)

    将每个人跑步的路径拆分成x->lca,lca->y两条路径分别考虑: 对于在点i的观察点,这个人(s->t)能被观察到的充要条件为: 1.直向上的路径:w[i]=dep[s]-dep ...

  4. NOIP2016 Day1 T2 天天爱跑步(树上差分,LCA)

    原文链接 原题链接 题目描述 小c同学认为跑步非常有趣,于是决定制作一款叫做<天天爱跑步>的游戏.<天天爱跑步>是一个养成类游戏,需要玩家每天按时上线,完成打卡任务. 这个游戏 ...

  5. 洛谷P1600 天天爱跑步——树上差分

    题目:https://www.luogu.org/problemnew/show/P1600 看博客:https://blog.csdn.net/clove_unique/article/detail ...

  6. 洛谷$P1600$ 天天爱跑步 树上差分

    正解:树上差分 解题报告: 传送门$QwQ$! 这题还挺妙的,,,我想了半天才会$kk$ 首先对一条链$S-T$,考虑先将它拆成$S-LCA$和$LCA-T$,分别做.因为总体上来说差不多接下来我就只 ...

  7. BZOJ 4719--天天爱跑步(LCA&差分)

    4719: [Noip2016]天天爱跑步 Time Limit: 40 Sec  Memory Limit: 512 MBSubmit: 1464  Solved: 490[Submit][Stat ...

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

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

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

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

随机推荐

  1. debian 9 安装后需做的几件事

    debian 9 安装后需做的几件事 安装环境:X86 >> Debian 9 Linux/GNU apt源更新 注意连上有线网络 刚安装好的debian系统中,/etc/apt/sour ...

  2. [Angular] Protect The Session Id with https and http only

    For the whole signup process. we need to Hash the password to create a password digest Store the use ...

  3. amazeui学习笔记--css(常用组件6)--图标Icon

    amazeui学习笔记--css(常用组件6)--图标Icon 一.总结 1.关注用法即可:在 HTML 上添加添加 am-icon-{图标名称} class. <span class=&quo ...

  4. Js里面的arguments

    了解这个对象之前先来认识一下javascript的一些功能: 其实Javascript并没有重载函数的功能,但是Arguments对象能够模拟重载.Javascrip中国每个函数都会有一个Argume ...

  5. 【hdu 6214】Smallest Minimum Cut

    [链接] 我是链接,点我呀:) [题意] 求最小割中最少的边数. [题解] 模板题 [代码] const int INF = 1e9; const int maxn = 1e3 + 7; const ...

  6. vue中的select框的值动态绑定

    <--这两种写法效果一样--> 1: <select v-model="wxStatus"> <option label="已添加" ...

  7. 【习题 3-2 UVA - 1586】Molar mass

    [链接] 我是链接,点我呀:) [题意] 在这里输入题意 [题解] 模拟. 主要是找那个数字. [代码] #include <bits/stdc++.h> using namespace ...

  8. TOP全异步模式

    Top全异步方式调用技术方案 背景:目前top通过servlet3.0技术结合异步管道化框架做到半异步调用,半异步调用采用异步线程同步调用后端的方式来做api call @飞不起的奥特曼 的部分文档) ...

  9. 关于spring获取webApplication.getBean多种途径和简单解释

    ApplicationContext ac1 = new FileSystemXmlApplicationContext("com/spark/system/applicationConte ...

  10. IWorkSpace与IWorkSpaceFactory与IWorkSpaceEdit

    樱木 原文 IWorkSpace与IWorkSpaceFactory与IWorkSpaceEdit 1.IWorkSpace是一个容器,里面存放着各种空间数据和非空间数据,比如Featureclass ...