BZOJ 4719 [Noip2016]天天爱跑步 ——树链剖分
一直以为自己当时是TLE了,但是再看发现居然WA?
然后把数组扩大一倍,就A掉了。QaQ
没什么好说的。一段路径分成两段考虑,上升的一段深度+时间是定值,下降的一段深度-时间是定值,然后打标记统计即可。
发现大概是统计数组因为深度+时间太大炸掉了。
现在想想,当时没有对拍,真是后怕。
#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
int read()
{
int ret=0,f=1; char ch=getchar();
while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
while (ch>='0'&&ch<='9'){ret*=10;ret+=ch-'0';ch=getchar();}
return ret*f;
}
int n,m;
int h[600005],to[600002],ne[600002],en=0;
void add(int a,int b)
{
to[en]=b;
ne[en]=h[a];
h[a]=en++;
}
int root=1,dep[600005],fa[600005],son[600005],siz[600005];
vector <int> st[600005];
vector <int> ed[600005];
void dfs1(int k)
{
siz[k]=1;
for (int i=h[k];i>=0;i=ne[i])
if (to[i]!=fa[k]){
fa[to[i]]=k;
dep[to[i]]=dep[k]+1;
dfs1(to[i]);
siz[k]+=siz[to[i]];
if (siz[son[k]]<siz[to[i]]) son[k]=to[i];
}
}
int pos[600005],top[600005],li[600005],cnt=0;
void dfs2(int k,int tp)
{
top[k]=tp;
cnt++; pos[k]=cnt; li[cnt]=k;
if (!son[k]) return ;
if (son[k]) dfs2(son[k],tp);
for (int i=h[k];i>=0;i=ne[i])
if (to[i]!=son[k]&&to[i]!=fa[k]) dfs2(to[i],to[i]);
return ;
}
int w[600005];
vector <int> que[600005];
vector <int> bel[600005];
struct node{int s,t,lca;}q[600005];
int vis[600005],f[600005];
int gf(int k)
{
if (f[k]==k) return k;
else return f[k]=gf(f[k]);
}
void dfs3(int k)
{
f[k]=k;
vis[k]=1;
for (int i=0;i<que[k].size();++i)
if (vis[que[k][i]]) q[bel[k][i]].lca=gf(que[k][i]);
for (int i=h[k];i>=0;i=ne[i])
if (fa[k]!=to[i]){
dfs3(to[i]);
f[to[i]]=k;
}
return ;
}
void lca(){dfs3(1);}
int s[600005],ans[600005];
void add(int a,int b,int fn)
{
if (dep[a]<dep[b]) swap(a,b);
while (top[a]!=top[b])
{
if (dep[top[a]]<dep[top[b]]) swap(a,b);
st[pos[top[a]]].push_back(fn);
ed[pos[a]+1].push_back(fn);
a=fa[top[a]];
}
if (dep[a]<dep[b]) swap(a,b);
st[pos[b]].push_back(fn);
ed[pos[a]+1].push_back(fn);
}
void cal1()
{
for (int i=1;i<=n;++i) w[i]-=dep[i];
for (int i=1;i<=m;++i)
{
int at;
at=dep[q[i].s]-dep[q[i].lca];
if (w[q[i].lca]==at-dep[q[i].lca]) ans[q[i].lca]--;
add(q[i].lca,q[i].t,at-dep[q[i].lca]);
}
for (int i=1;i<=n;++i)
{
for (int j=0;j<st[i].size();++j) s[st[i][j]]++;
for (int j=0;j<ed[i].size();++j) s[ed[i][j]]--;
ans[li[i]]+=s[w[li[i]]];
}
}
void cal2()
{
for (int i=0;i<=n+5;++i) st[i].clear(),ed[i].clear();
memset(s,0,sizeof s);
for (int i=1;i<=n;++i) w[i]+=2*dep[i];
for (int i=1;i<=m;++i)
{
int at;
at=dep[q[i].s]-dep[q[i].lca];
add(q[i].lca,q[i].s,at+dep[q[i].lca]);
}
for (int i=1;i<=n;++i)
{
for (int j=0;j<st[i].size();++j) s[st[i][j]]++;
for (int j=0;j<ed[i].size();++j) s[ed[i][j]]--;
ans[li[i]]+=s[w[li[i]]];
}
}
void out()
{
for (int i=1;i<n;++i)printf("%d ",ans[i]); printf("%d",ans[n]);
}
int main()
{
memset(h,-1,sizeof h);
n=read(); m=read();
for (int i=1;i<n;++i)
{
int a,b;
a=read(); b=read();
add(a,b); add(b,a);
}
for (int i=1;i<=n;++i)
w[i]=read();
for (int i=1;i<=m;++i)
{
int a,b;
a=read(); b=read();
que[a].push_back(b);
que[b].push_back(a);
bel[a].push_back(i);
bel[b].push_back(i);
q[i].s=a;q[i].t=b;
}
dfs1(root);
dfs2(root,root);
lca();
cal1();
cal2();
out();
return 0;
}
BZOJ 4719 [Noip2016]天天爱跑步 ——树链剖分的更多相关文章
- NOIP2016提高组Day1T2 天天爱跑步 树链剖分 LCA 倍增 差分
原文链接https://www.cnblogs.com/zhouzhendong/p/9275606.html 题目传送门 - 洛谷P1600 题目传送门 - LOJ#2359 题目传送门 - Vij ...
- bzoj 4719: [Noip2016]天天爱跑步
Description 小c同学认为跑步非常有趣,于是决定制作一款叫做<天天爱跑步>的游戏.?天天爱跑步?是一个养成类游戏,需要 玩家每天按时上线,完成打卡任务.这个游戏的地图可以看作一一 ...
- bzoj 4719: [Noip2016]天天爱跑步【树上差分+dfs】
长久以来的心理阴影?但是其实非常简单-- 预处理出deep和每组st的lca,在这里我简单粗暴的拿树剖爆算了 然后考虑对于一组s t lca来说,被这组贡献的观察员x当且仅当: x在s到lca的路径上 ...
- [BZOJ 2243] [SDOI 2011] 染色 【树链剖分】
题目链接:BZOJ - 2243 题目分析 树链剖分...写了200+行...Debug了整整一天+... 静态读代码读了 5 遍 ,没发现错误,自己做小数据也过了. 提交之后全 WA . ————— ...
- BZOJ 4568 [Scoi2016]幸运数字(树链剖分 + 异或线性基)
题目链接 BZOJ 4568 考虑树链剖分+线段树维护每一段区域的异或线性基 对于每个询问,求出该点集的异或线性基.然后求一下这个线性基里面能异或出的最大值即可. #include <bits ...
- Bzoj 4196: [Noi2015]软件包管理器 树链剖分
4196: [Noi2015]软件包管理器 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 721 Solved: 419[Submit][Statu ...
- BZOJ 1969: [Ahoi2005]LANE 航线规划( 树链剖分 )
首先我们要时光倒流, 倒着做, 变成加边操作维护关键边. 先随意搞出一颗树, 树上每条边都是关键边(因为是树, 去掉就不连通了)....然后加边(u, v)时, 路径(u, v)上的所有边都变成非关键 ...
- bzoj 4196 [Noi2015]软件包管理器 (树链剖分+线段树)
4196: [Noi2015]软件包管理器 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 2852 Solved: 1668[Submit][Sta ...
- BZOJ - 3631 松鼠的新家 (树链剖分)
题目链接 树链剖分基础题,路径权值修改+差分 #include<bits/stdc++.h> using namespace std; typedef long long ll; ; in ...
随机推荐
- webpack前端构建工具学习总结(四)之自动化生成项目中的html页面
接续上文:webpack前端构建工具学习总结(三)之webpack.config.js配置文件 插件的介绍文档:https://www.npmjs.com/package/html-webpack-p ...
- 洛谷 P1330 封锁阳光大学
题目描述 曹是一只爱刷街的老曹,暑假期间,他每天都欢快地在阳光大学的校园里刷街.河蟹看到欢快的曹,感到不爽.河蟹决定封锁阳光大学,不让曹刷街. 阳光大学的校园是一张由N个点构成的无向图,N个点之间由M ...
- EL1008E: Property or field 'timestamp' cannot be found on object of type 'java.util.HashMap
2018-06-22 09:50:19.488 INFO 20096 --- [nio-8081-exec-2] o.a.c.c.C.[Tomcat].[localhost].[/] : ...
- python基础面试题整理---从零开始 每天十题(01)
最近在弄flask的东西,好久没写博客的,感觉少了点什么,感觉被别人落下好多,可能渐渐的养成了写博客的习惯吧.也是自己想学的东西太多了(说白了就是基础太差了,只是know how,不能做到konw w ...
- AC自动机讲解+[HDU2222]:Keywords Search(AC自动机)
首先,有这样一道题: 给你一个单词W和一个文章T,问W在T中出现了几次(原题见POJ3461). OK,so easy~ HASH or KMP 轻松解决. 那么还有一道例题: 给定n个长度不超过50 ...
- Java第十二次作业:什么是一维数组?什么是对象数组?吃金币游戏2.0版 新增炸弹功能 新增游戏倒计时功能 新增胜利失败检测功能 使用如鹏游戏引擎制作窗体 一维数组设置金币
什么是数组? 数组的定义:是用统一的名字代表这批数据,用序号来区分各个数据.数组是无序的数据元素按有序的下标组成的集合,分配固定空间大小的一种容器. 如何理解:其实就是一个同时放很多数据的变量. a= ...
- Mac 电源管理
在安装BatteryManager后,可以删除NullPowerMananger,AppleIntelPowerMananger, AppleIntelPowerClientMananger三个kex ...
- visibilitychange 标签可见性
var pageVisibility = document.visibilityState;// 监听 visibility change 事件document.addEventListener('v ...
- js基本语法之 值类型(数据类型)(变量类型)
一.不可改变的原始值(栈数据)(五个) 数字(number),字符串(string),布尔值(boolean),undefined,null 其中;undefined是未定义的意思,而null是空的意 ...
- python面向对象编程(OOP)
python作为一种解释性语言,其主要的编程方式就是面向对象,而且python的框架django也是主要面向对象的编程. 类(class)和对象(object) 类(class)是用来描述具有相同属性 ...