NOIP2016 天天爱跑步 线段树合并_桶_思维题
Description
Input
Output
输出1行N 个整数,第个整数表示结点的观察员可以观察到多少人。
Code:
#include<bits/stdc++.h>
#define setIO(s) freopen(s".in","r",stdin)
#define maxn 400000
#define M 1000000
#define plus pl
#define minus mi
using namespace std;
vector<int>plus[M],minus[M];
int n,m;
int hd[maxn<<1],to[maxn<<1],nx[maxn<<1],edges;
int dep[maxn],F[22][maxn],tim[maxn],st[maxn],ed[maxn];
int answer[maxn];
void add(int u,int v){
nx[++edges]=hd[u],hd[u]=edges,to[edges]=v;
}
void dfs1(int u,int fa){
dep[u]=dep[fa]+1;
F[0][u]=fa;
for(int i=1;i<22;++i) F[i][u]=F[i-1][F[i-1][u]];
for(int v=hd[u];v;v=nx[v])
if(to[v]!=fa) dfs1(to[v],u);
}
int lca(int a,int b){
if(dep[a]>dep[b]) swap(a,b);
if(dep[a]!=dep[b]) for(int i=21;i>=0;--i) if(dep[F[i][b]]>=dep[a]) b=F[i][b];
if(a==b) return a;
for(int i=21;i>=0;--i) if(F[i][a]!=F[i][b]) a=F[i][a],b=F[i][b];
return F[0][a];
}
int tot,root[maxn];
struct Node{ int l,r,w; }node[maxn*11];
#define mid ((l+r)>>1)
void ins(int &o,int k,int delta,int l,int r){
if(!o) o=++tot;
node[o].w+=delta;
if(l==r) return;
else{
if(k<=mid) ins(node[o].l,k,delta,l,mid);
else ins(node[o].r,k,delta,mid+1,r);
}
}
int merge(int x,int y){
if(!x||!y) return x+y;
node[x].w+=node[y].w;
node[x].l=merge(node[x].l,node[y].l);
node[x].r=merge(node[x].r,node[y].r);
return x;
}
int query(int x,int l,int r,int k){
if(!x) return 0;
if(l==r) return node[x].w;
if(k<=mid) return query(node[x].l,l,mid,k);
else return query(node[x].r,mid+1,r,k);
}
void dfs2(int u){
for(int sz=plus[u].size(),i=0;i<sz;++i) ins(root[u],plus[u][i],1,1,M);
for(int sz=minus[u].size(),i=0;i<sz;++i) ins(root[u],minus[u][i],-1,1,M);
plus[u].clear(),minus[u].clear();
for(int v=hd[u];v;v=nx[v])
if(to[v]!=F[0][u]) dfs2(to[v]),root[u]=merge(root[u],root[to[v]]);
answer[u]+=query(root[u],1,M,tim[u]+dep[u]);
}
void dfs3(int u){
for(int sz=plus[u].size(),i=0;i<sz;++i) ins(root[u],plus[u][i],1,1,M);
for(int sz=minus[u].size(),i=0;i<sz;++i) ins(root[u],minus[u][i],-1,1,M);
for(int v=hd[u];v;v=nx[v])
if(to[v]!=F[0][u]) dfs3(to[v]),root[u]=merge(root[u],root[to[v]]);
answer[u]+=query(root[u],1,M,dep[u]-tim[u]+maxn);
}
void up(){
for(int i=1,c;i<=m;++i){
c=lca(st[i],ed[i]);
if(dep[c]<dep[st[i]]) {
plus[st[i]].push_back(dep[st[i]]);
minus[F[0][c]].push_back(dep[st[i]]);
}
}
dfs2(1);
memset(node,0,sizeof(node));
memset(root,0,sizeof(root));
tot=0;
}
void down(){
for(int i=1,c,pre;i<=m;++i) {
c=lca(st[i],ed[i]);
if(dep[c]<dep[ed[i]]) {
pre=dep[st[i]]-dep[c];
if(st[i]==c){
plus[ed[i]].push_back(dep[c]-pre+maxn);
minus[F[0][c]].push_back(dep[c]-pre+maxn);
}else {
plus[ed[i]].push_back(dep[c]-pre+maxn);
minus[c].push_back(dep[c]-pre+maxn);
}
}
}
dfs3(1);
}
int main(){
//setIO("input");
scanf("%d%d",&n,&m);
for(int i=1,a,b;i<n;++i) scanf("%d%d",&a,&b),add(a,b),add(b,a);
dfs1(1,0);
for(int i=1;i<=n;++i) scanf("%d",&tim[i]);
for(int i=1;i<=m;++i) {
scanf("%d%d",&st[i],&ed[i]);
if(st[i]==ed[i] && tim[st[i]]==0) ++answer[st[i]];
}
up(),down();
for(int i=1;i<=n;++i) printf("%d ",answer[i]);
return 0;
}
NOIP2016 天天爱跑步 线段树合并_桶_思维题的更多相关文章
- [NOIp2016]天天爱跑步 线段树合并
[NOIp2016]天天爱跑步 LG传送门 作为一道被毒瘤出题人们玩坏了的NOIp经典题,我们先不看毒瘤的"动态爱跑步"和"天天爱仙人掌",回归一下本来的味道. ...
- 2018.08.09 bzoj4719: [Noip2016]天天爱跑步(树链剖分)
传送门 话说开始上文化课之后写题时间好少啊. 这道题将一个人的跑步路线拆成s->lca,lca->t,然后对于第一段上坡路径要经过的点,当前这个人能对它产生贡献当且仅当dep[s]-dep ...
- [Vani有约会]雨天的尾巴 线段树合并
[Vani有约会]雨天的尾巴 LG传送门 线段树合并入门好题. 先别急着上线段树合并,考虑一下这题的暴力.一看就是树上差分,对于每一个节点统计每种救济粮的数量,再一遍dfs把差分的结果统计成答案.如果 ...
- CF700E Cool Slogans——SAM+线段树合并
RemoteJudge 又是一道用线段树合并来维护\(endpos\)的题,还有一道见我的博客CF666E 思路 先把\(SAM\)建出来 如果两个相邻的串\(s_i\)和\(s_{i+1}\)要满足 ...
- LOJ #2359. 「NOIP2016」天天爱跑步(倍增+线段树合并)
题意 LOJ #2359. 「NOIP2016」天天爱跑步 题解 考虑把一个玩家的路径 \((x, y)\) 拆成两条,一条是 \(x\) 到 \(lca\) ( \(x, y\) 最近公共祖先) 的 ...
- [NOIP2016 DAY1 T2]天天爱跑步-[差分+线段树合并][解题报告]
[NOIP2016 DAY1 T2]天天爱跑步 题面: B[NOIP2016 DAY1]天天爱跑步 时间限制 : - MS 空间限制 : 565536 KB 评测说明 : 2s Description ...
- [NOIP2016]天天爱跑步(树上差分+线段树合并)
将每个人跑步的路径拆分成x->lca,lca->y两条路径分别考虑: 对于在点i的观察点,这个人(s->t)能被观察到的充要条件为: 1.直向上的路径:w[i]=dep[s]-dep ...
- BZOJ4719 NOIP2016天天爱跑步(线段树合并)
线段树合并的话这个noip最难题就是个裸题了. 注意merge最后return x,以及如果需要区间查询的话这里还需要up,无数次死于这里. #include<iostream> #inc ...
- 【bzoj4719】[Noip2016]天天爱跑步 权值线段树合并
题目描述 给出一棵n个点的树,以及m次操作,每次操作从起点向终点以每秒一条边的速度移动(初始时刻为0),最后对于每个点询问有多少次操作在经过该点的时刻为某值. 输入 第一行有两个整数N和M .其中N代 ...
随机推荐
- 在ecshop模板使用自定义函数
https://blog.csdn.net/shaolinld/article/details/46400485 在ecshop模板使用自定义函数 可以增加自定义函数,在模板直接调用,例如: {$us ...
- C#关键字的个人理解与注释
C#关键字注释:abstract:抽象as:类型转换(返回转换结果)base:基类bool:布尔类型break:条件中断语句byte:字节case:条件语句catch:异常捕获后执行char:16 位 ...
- 在Win32 Application 环境下实现MFC窗口的创建
// Win32下MFC.cpp : Defines the entry point for the application.// #include "stdafx.h" clas ...
- 0208MySQL5.7之Group Replication
转自http://blog.itpub.net/29510932/viewspace-2055679/ MySQL Group Replication: Hello World! 对测试版(on la ...
- ExtJs之ExtJs.Model验证
没有搞好,但知道意思和配置,并且好像4和3的没有兼容性语法,所以网上找了新语法才出来了点东东. <!DOCTYPE html> <html> <head> < ...
- qsort快速排序
C库函数qsort七种使用方法示例 七种qsort排序方法<本文中排序都是采用的从小到大排序> 一.对int类型数组排序C++ / C 代码 int num[100]; Sample: i ...
- Column注解的的RetentionPolicy的属性值是RUTIME,这样注解处理器可以通过反射,获取到该注解的属性值,从而去做一些运行时的逻辑处理
1.Column注解的的RetentionPolicy的属性值是RUTIME,这样注解处理器可以通过反射,获取到该注解的属性值,从而去做一些运行时的逻辑处理 2. 自定义注解: 使用@interfac ...
- caffe Solve函数
下面来看Solver<Dtype>::Solve(const char* resume_file) solver.cpp template <typename Dtype> v ...
- UVA 10039 Railroads
这道题好吧,一开始便是拓扑排序的想法,搞了好久,试了多组测试数据,没错啊,可是没过...作孽啊,竟然忘了拓扑不能处理环,白浪费了一晚上... 只好用动态规划了.. DP[time][city]表示在t ...
- oracle强化练习题
编写简单的查询语句 1.显示部门表的所有信息 Select * from dept; 2.显示部门号码,部门名称 Selectdeptno,dname from dept; 3.显示下面字段及字符串的 ...