bzoj 4012: [HNOI2015]开店 主席树
Description
风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到
Input
第一行三个用空格分开的数 n、Q和A,表示树的大小、开店的方案个数和妖
Output
对于每个方案,输出一行表示方便值。
Sample Input
0 0 7 2 1 4 7 7 7 9
1 2 270
2 3 217
1 4 326
2 5 361
4 6 116
3 7 38
1 8 800
6 9 210
7 10 278
8 9 8
2 8 0
9 3 1
8 0 8
4 2 7
9 7 3
4 7 0
2 2 7
3 2 1
2 3 4
Sample Output
957
7161
9466
3232
5223
1879
1669
1282
0
HINT
满足 n<=150000,Q<=200000。对于所有数据,满足 A<=10^9
Source
哎,原来这是一道主席树的套路题啊,然而当年用动态点分治做码得要死。。。
对于每个u要求:
拆开一下,相当于是求:
前两项都可以通过dis数组的前缀和求解,问题的瓶颈在于后面一项。。。
这一项,1-lca的路径其实是1-v的路径和1-u的路径的交。。。
首先我们不考虑年龄限制,我们只需要把1--v路径上的点都打上标记,然后只要查询1--u路径上的和即可。。。
(跟LNOI的LCA类似,相当于是线段树上区间修改和区间求和。。。)
加上年龄限制的话,我们可以把点按照年龄大小依次加入,
然后可持久化一下,对与每个年龄都建一棵线段树(每个点都是在历史版本上修改log个区间)
然后查询的话直接利用前缀和的性质相减即可,lazy实现标记永久化。。。
(史上最强样例,没有之一。。。)
// MADE BY QT666
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<cstring>
#include<vector>
#define int long long
using namespace std;
typedef long long ll;
const int N=300050;
int head[N],to[N],nxt[N],v[N],cnt;
int size[N],top[N],son[N],deep[N],fa[N],dfn[N],tt,dis[N],d[N];
int sz,rt[N],ls[N*20],rs[N*20],sum[N*20],lazy[N*20],res[N],res2[N],num[N];
int hsh[N],a[N],tot,n,Q,A;
void dfs1(int x,int f){size[x]=1;deep[x]=deep[f]+1;for(int i=head[x];i;i=nxt[i]){int y=to[i];if(y!=f){d[y]=v[i];dis[y]=dis[x]+v[i];dfs1(y,x);fa[y]=x;size[x]+=size[y];if(size[y]>size[son[x]]) son[x]=y;}}
}
void dfs2(int x,int ff){top[x]=ff;dfn[x]=++tt;res[tt]=res[tt-1]+d[x];if(son[x]) dfs2(son[x],ff);for(int i=head[x];i;i=nxt[i]){int y=to[i];if(y!=fa[x]&&y!=son[x]) dfs2(y,y);}
}
void update(int x,int &y,int l,int r,int xl,int xr){y=++sz;ls[y]=ls[x],rs[y]=rs[x];lazy[y]=lazy[x];sum[y]=sum[x];if(xl<=l&&r<=xr){sum[y]+=(res[r]-res[l-1]);lazy[y]++;return;}int mid=(l+r)>>1;if(xr<=mid) update(ls[x],ls[y],l,mid,xl,xr);else if(xl>mid) update(rs[x],rs[y],mid+1,r,xl,xr);else update(ls[x],ls[y],l,mid,xl,mid),update(rs[x],rs[y],mid+1,r,mid+1,xr);sum[y]=sum[ls[y]]+sum[rs[y]]+lazy[y]*(res[r]-res[l-1]);
}
int query(int x,int y,int l,int r,int xl,int xr,int la){if(xl<=l&&r<=xr){return sum[y]-sum[x]+la*(res[r]-res[l-1]);}int mid=(l+r)>>1;la+=lazy[y]-lazy[x];if(xr<=mid) return query(ls[x],ls[y],l,mid,xl,xr,la);else if(xl>mid) return query(rs[x],rs[y],mid+1,r,xl,xr,la);else return query(ls[x],ls[y],l,mid,xl,mid,la)+query(rs[x],rs[y],mid+1,r,mid+1,xr,la);
}
vector<int> p[N];
void insert(int x,int age){while(x){update(rt[age],rt[age],1,n,dfn[top[x]],dfn[x]);x=fa[top[x]];}
}
int ask(int x,int l,int r){int ret=0;while(x){ret+=query(rt[l-1],rt[r],1,n,dfn[top[x]],dfn[x],0);x=fa[top[x]];}return ret;
}
void lnk(int x,int y,int z){to[++cnt]=y,nxt[cnt]=head[x],v[cnt]=z,head[x]=cnt;to[++cnt]=x,nxt[cnt]=head[y],v[cnt]=z,head[y]=cnt;
}
main(){scanf("%lld%lld%lld",&n,&Q,&A);for(int i=1;i<=n;i++) scanf("%lld",&a[i]),hsh[++tot]=a[i];sort(hsh+1,hsh+1+tot);tot=unique(hsh+1,hsh+1+tot)-hsh-1;for(int i=1;i<=n;i++) a[i]=lower_bound(hsh+1,hsh+1+tot,a[i])-hsh,p[a[i]].push_back(i);for(int i=1;i<n;i++){int x,y,z;scanf("%lld%lld%lld",&x,&y,&z);lnk(x,y,z);}dfs1(1,0);dfs2(1,1);for(int i=1;i<=tot;i++){rt[i]=rt[i-1];res2[i]=res2[i-1];num[i]=num[i-1];for(int j=0;j<p[i].size();j++) insert(p[i][j],i),res2[i]+=dis[p[i][j]],num[i]++;}int ans=0;while(Q--){int u,a,b;scanf("%lld%lld%lld",&u,&a,&b);int l=min((a+ans)%A,(b+ans)%A),r=max((a+ans)%A,(b+ans)%A);l=lower_bound(hsh+1,hsh+1+tot,l)-hsh,r=upper_bound(hsh+1,hsh+1+tot,r)-hsh-1;ans=res2[r]-res2[l-1]+(num[r]-num[l-1])*dis[u]-2*ask(u,l,r);printf("%lld\n",ans);}return 0;
}
转载于:https://www.cnblogs.com/qt666/p/7265067.html
bzoj 4012: [HNOI2015]开店 主席树相关推荐
- BZOJ 4012: [HNOI2015]开店 -- 动态树分治
4012: [HNOI2015]开店 Time Limit: 70 Sec Memory Limit: 512 MB Submit: 1463 Solved: 635 [Submit][Statu ...
- bzoj 4012: [HNOI2015]开店
Description 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现 ...
- 4012: [HNOI2015]开店
4012: [HNOI2015]开店 Time Limit: 70 Sec Memory Limit: 512 MB Submit: 1241 Solved: 536 [Submit][Statu ...
- [HNOI2015]开店(树剖+主席树+标记永久化)
听说正解点分树?我不会就对了 此题是 \([LNOI2014]LCA\) 强化版,也是差分一下,转化为区间加区间和 不过权值有大小要求,那么我们按照权值排序,依次加入主席树,询问的时候 \(lower ...
- BZOJ.3524.[POI2014]Couriers(主席树)
题目链接 //119964kb 5780ms //主席树裸题啊.. #include <cstdio> #include <cctype> #define gc() getch ...
- BZOJ 4009 HNOI2015 接水果 树套树
题目大意:给定一棵树和m<script type="math/tex" id="MathJax-Element-32">m</script&g ...
- BZOJ 3277 串 BZOJ 3473 字符串 (广义后缀自动机、时间复杂度分析、启发式合并、线段树合并、主席树)...
标签那么长是因为做法太多了... 题目链接: (bzoj 3277) https://www.lydsy.com/JudgeOnline/problem.php?id=3277 (bzoj 3473) ...
- BZOJ 3218 UOJ #77 A+B Problem (主席树、最小割)
BZOJ 3218 UOJ #77 A+B Problem (主席树.最小割) 大名鼎鼎的A+B Problem, 主席树优化最小割-- 调题死活调不对,一怒之下改了一种写法交上去A了,但是改写法之后 ...
- BZOJ 2588 Spoj 10628 Count on a tree | 树上主席树
BZOJ 2588 Count on a tree 题面 求树上两点之间路径上第k大的点权. 题解 一开始看到这道题觉得是树剖,然后又听说是主席树,然后以为是主席树+树剖,差点没吓死-- 然后发现,如 ...
最新文章
- 作为程序员,你是否曾经想过写一本书?
- 1370亿参数、接近人类水平,谷歌重磅推出对话AI模型LaMDA
- #define、#undef、#ifdef、#ifndef、#if、#elif、#else、#endif、defined解释
- 用tcc自己编译自己
- TP5 连接oracle数据库
- 自我学习的技巧和建议
- 使用DynamoDBMapper插入DynamoDB项目
- 雄安大学真的来了!占地总面积2550亩!
- [转载] python3基础:异常处理及python常见异常类型总结
- 数学建模——遗传算法(GA)
- 20145235 《Java程序设计》第6周学习总结
- 操作系统课程设计-文件管理系统(JAVA)
- 电子书,世界名著,中国名著
- Relational Knowledge Distillation解读
- 大学生php实训总结_php实训报告.doc
- Windows安装ElasticSearch 7.3.0
- Linux下用脚本命令打开文档、表格、PPT
- 明日之后什么服务器物品最便宜,明日之后:玩家晒物价最低的服务器,看到uzi标价,观众直接酸了...
- JS---------------网页版的消灭星星
- ABAQUS如何输出应力应变曲线(XY曲线)