【LOJ】#2116. 「HNOI2015」开店
题解
一道我觉得和二叉树没有关系的题……
因为直接上点分就过了,虽然很慢,而且代码很长
你需要记录一个点分树,对于每个点分树的重心,记录一下上一次进行分割时树的重心以及这个重心和上一次重心所连接的点以及连接的边的距离
然后计算这个重心和所在的树到上一个重心节点路径和的前缀和,还有节点个数和
处理每棵树的节点路径前缀和和节点个数前缀和
处理的时候枚举这个点经过的每个重心,统计过重心的路径和即可
代码
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
#define MAXN 150005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef long double db;
typedef unsigned int u32;
template<class T>
void read(T &res) {
res = 0;char c = getchar();T f = 1;
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 + c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {putchar('-');x = -x;}
if(x >= 10) out(x / 10);
putchar('0' + x % 10);
}
int N,Q,A,num[MAXN],fa[MAXN],siz[MAXN],son[MAXN],dis[MAXN];
vector<int> age[MAXN],aux[MAXN],dis_fa[MAXN],siztr[MAXN],siztr_fa[MAXN];
vector<int64> sum[MAXN],sum_fa[MAXN];
int que[MAXN],ql,qr,id[MAXN],dfn[MAXN],idx,pre[MAXN],pos[MAXN],pre_dis[MAXN];
bool vis[MAXN];
struct node {
int to,next;
int64 val;
}E[MAXN * 2];
int head[MAXN],sumE;
void add(int u,int v,int64 c) {
E[++sumE].next = head[u];
E[sumE].to = v;
E[sumE].val = c;
head[u] = sumE;
}
int CalcG(int st) {
que[ql = qr = 1] = st;
fa[st] = 0;
while(ql <= qr) {
int u = que[ql++];siz[u] = 1;son[u] = 0;
for(int i = head[u] ; i ; i = E[i].next) {
int v = E[i].to;
if(v != fa[u] && !vis[v]) {
que[++qr] = v;fa[v] = u;
}
}
}
int res = que[qr];
for(int i = qr ; i >= 1 ; --i) {
int u = que[i];
if(fa[u]) {
siz[fa[u]] += siz[u];
son[fa[u]] = max(son[fa[u]],siz[u]);
}
son[u] = max(son[u],qr - siz[u]);
if(son[u] < son[res]) res = u;
}
return res;
}
void Calc(int st,int64 val) {
que[ql = qr = 1] = st;fa[st] = 0;dis[st] = val;
while(ql <= qr) {
int u = que[ql++];
for(int i = head[u] ; i ; i = E[i].next) {
int v = E[i].to;
if(!vis[v] && v != fa[u]) {
dis[v] = dis[u] + E[i].val;
fa[v] = u;
que[++qr] = v;
}
}
}
}
int64 S(int u,int R) {
int siz = aux[u].size();
int64 res = 0;
for(int i = 0 ; i < siz - 1; ++i) {
int z = aux[u][i],y = aux[u][i + 1];
int t = upper_bound(age[z].begin(),age[z].end(),R) - age[z].begin() - 1;
res += 1LL * (siztr[z][t] - siztr_fa[y][t]) * dis_fa[u][i] + sum[z][t] - sum_fa[y][t];
}
res += sum[u][upper_bound(age[u].begin(),age[u].end(),R) - age[u].begin() - 1];
return res;
}
bool cmp(int a,int b) {
return num[a] < num[b];
}
void dfs_G(int u,int f,int64 c) {
int G = CalcG(u);
dfn[++idx] = G;
pre[G] = f;pos[G] = u;pre_dis[G] = c;vis[G] = 1;
for(int i = head[G] ; i; i = E[i].next) {
int v = E[i].to;
if(!vis[v]) dfs_G(v,G,E[i].val);
}
}
void Init() {
read(N);read(Q);read(A);
for(int i = 1 ; i <= N ; ++i) read(num[i]);
for(int i = 1 ; i < N ; ++i) {
int u,v;int64 c;
read(u);read(v);read(c);
add(u,v,c);add(v,u,c);
}
dfs_G(1,0,0);
memset(vis,0,sizeof(vis));
for(int i = 1 ; i <= idx ; ++i) {
int G = dfn[i];
if(pre[G]) {
Calc(pos[G],pre_dis[G]);
sum_fa[G].resize(sum[pre[G]].size());
siztr_fa[G].resize(sum[pre[G]].size());
for(int i = 1 ; i <= qr ; ++i) {
int u = que[i];
int t = lower_bound(age[pre[G]].begin(),age[pre[G]].end(),num[u]) - age[pre[G]].begin();
sum_fa[G][t] += dis[u];
siztr_fa[G][t]++;
}
int s = sum_fa[G].size();
for(int i = 1 ; i < s ; ++i) {sum_fa[G][i] += sum_fa[G][i - 1];siztr_fa[G][i] += siztr_fa[G][i - 1];}
}
Calc(G,0);
vis[G] = 1;
for(int i = 1 ; i <= qr ; ++i) id[i] = que[i];
sort(id + 1,id + qr + 1,cmp);
int s = 0;
sum[G].pb(0);age[G].pb(-1);siztr[G].pb(0);
for(int i = 1 ; i <= qr ; ++i) {
int u = id[i];
aux[u].pb(G);
dis_fa[u].pb(dis[u]);
if(i == 1 || num[u] != num[id[i - 1]]) {
age[G].pb(num[u]);
sum[G].pb(0);siztr[G].pb(0);++s;
}
sum[G][s] += dis[u];
siztr[G][s]++;
}
for(int i = 1 ; i <= s ; ++i) {sum[G][i] += sum[G][i - 1];siztr[G][i] += siztr[G][i - 1];}
age[G].pb(A + 1);sum[G].pb(sum[G][s]);
}
}
void Solve() {
int64 ans = 0;
int u,a,b;
for(int i = 1 ; i <= Q ; ++i) {
read(u);read(a);read(b);
a = (a + ans) % A;b = (b + ans) % A;
if(a > b) swap(a,b);
ans = S(u,b) - S(u,a - 1);
out(ans);enter;
}
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Init();
Solve();
return 0;
}
(谁来拯救一下我越写越慢越写越长的数据结构题= =)
【LOJ】#2116. 「HNOI2015」开店的更多相关文章
- loj #2116. 「HNOI2015」开店
#2116. 「HNOI2015」开店 题目描述 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的想法当然非 ...
- LOJ #2116 Luogu P3241「HNOI2015」开店
好久没写数据结构了 来补一发 果然写的时候思路极其混乱.... LOJ #2116 Luogu P3241 题意 $ Q$次询问,求树上点的颜色在$ [L,R]$中的所有点到询问点的距离 强制在线 询 ...
- [loj2116]「HNOI2015」开店 动态点分治
4012: [HNOI2015]开店 Time Limit: 70 Sec Memory Limit: 512 MBSubmit: 2452 Solved: 1089[Submit][Status ...
- 「HNOI2015」开店(树链剖分, 主席树)
/* 考虑将所求的值拆分 记每个点到根的路径长度为dis_i, 那么我们要求的就是\sum_{i = l} ^ r dis_i + dis[u] * (r - l + 1) - 2\sum_{i = ...
- LOJ #2135. 「ZJOI2015」幻想乡战略游戏(点分树)
题意 给你一颗 \(n\) 个点的树,每个点的度数不超过 \(20\) ,有 \(q\) 次修改点权的操作. 需要动态维护带权重心,也就是找到一个点 \(v\) 使得 \(\displaystyle ...
- Loj #2192. 「SHOI2014」概率充电器
Loj #2192. 「SHOI2014」概率充电器 题目描述 著名的电子产品品牌 SHOI 刚刚发布了引领世界潮流的下一代电子产品--概率充电器: 「采用全新纳米级加工技术,实现元件与导线能否通电完 ...
- Loj #3096. 「SNOI2019」数论
Loj #3096. 「SNOI2019」数论 题目描述 给出正整数 \(P, Q, T\),大小为 \(n\) 的整数集 \(A\) 和大小为 \(m\) 的整数集 \(B\),请你求出: \[ \ ...
- Loj #3093. 「BJOI2019」光线
Loj #3093. 「BJOI2019」光线 题目描述 当一束光打到一层玻璃上时,有一定比例的光会穿过这层玻璃,一定比例的光会被反射回去,剩下的光被玻璃吸收. 设对于任意 \(x\),有 \(x\t ...
- Loj #3089. 「BJOI2019」奥术神杖
Loj #3089. 「BJOI2019」奥术神杖 题目描述 Bezorath 大陆抵抗地灾军团入侵的战争进入了僵持的阶段,世世代代生活在 Bezorath 这片大陆的精灵们开始寻找远古时代诸神遗留的 ...
随机推荐
- CDN中,字体文件的跨域问题和解决
@font-face是CSS3中的一个特性,可以把自己定义的Web字体嵌入到网页中,随着@font-face,越来越多的网页采用字体图标作为网页中的小图形. 比如Bootstrap就采用了Glyphi ...
- Kubernetes 部署 1.9.7 高可用版
转载于https://codegreen.cn/2018/08/30/kubernetes-cluster-1.9.7/ 前言 在部署之前,首先感谢 手动搭建高可用的kubernetes 集群 博文的 ...
- Go_22: Golang 命令行 test 应用
1. 测试某一个包下的所有测试用例 cd /myGoProject/src/go-test/utils go test -v "-v" 参数 go test -v ... 表示无论 ...
- C语言复习---获取矩阵的对角和
#include <stdio.h> #include <stdlib.h> #include <math.h> #define MAX 5 int main() ...
- 拼接html a标签字符串,onClick传递两个字符串类型参数写法
在参数传递过程中字符串类型的参数要有引号,我一开始拼接的完成后,没有想到要加引号,后来想到了这一问题,可是怎么拼都不对,于是就搜了很多拼接的例子,发现并没有几个能借鉴的,最后终于在一个人的博客中看到, ...
- 2018年9月28日CCPC秦皇岛站参赛总结
day1: 被中间结果超出int范围给叉了,立刻意识到了自己的弱小以及校赛出题的时候是怎么叉别人的 day2: 签到签了40分钟,谨慎一些还是很好的,机子重启耽误了一些时间 一道暴力+LCS的简单题被 ...
- jdk1.8.0_45源码解读——Set接口和AbstractSet抽象类的实现
jdk1.8.0_45源码解读——Set接口和AbstractSet抽象类的实现 一. Set架构 如上图: (01) Set 是继承于Collection的接口.它是一个不允许有重复元素的集合.(0 ...
- MongoDB-3.4搭建副本集
搭建副本集 1:首先创建3台虚拟机作为配置环境 IP1:192.168.101.175 IP2:192.168.101.176 IP3:192.168.101.177 2.下载MongoDB 3.4版 ...
- iOS编码规范参考
目录 注释 1.1 多行注释 1.2 单行注释 1.3 函数的注释 命名 2.1 常量的命名 2.2 函数的命名 2.3 变量的命名 2.3.1 成员变量 2.3.2 公 ...
- python 多进程的启动和代码执行顺序
对照着廖雪峰的网站学习Python遇到些问题: 在进程中,父进程创建子进程时发现,显示不是按照顺序显示,疑问? 参照代码如下: from multiprocessing import Pool imp ...