BZOJ_4867_[Ynoi2017]舌尖上的由乃_分块+dfs序

Description

由乃为了吃到最传统最纯净的美食,决定亲自开垦一片菜园。现有一片空地,由乃已经规划n个地点准备种上蔬菜。最新鲜的蔬菜需有最甘甜井水的灌溉,因此由乃将要打出两口井,分别记为井A、井B。现在问题来了,由乃可是一周目的神,为何要打井?是谁想出来的这些题面?由乃不善于搞事情,于是提出以下几个方法,再根据这些方法找出题人。方法如下:
1. 做完这个出题人出的所有题
2. 做完所有数据结构题
3. 出一道奇怪的数据结构题,而且卡常卡死所有做题的
至于为什么这样能找到出题人呢。。。我也不信她能找到我。。。(能找到岂不是更好?)
因为由乃是神,所有她有很多秒时间,她终于做完了世界上左右的数据结构题
于是该她出题了,她左思右想,出了一个简单数据结构题:
给你一个n个点的有根树,1为根,带边权,有m次操作。
1、求x的子树中第k小的深度的值,如果子树中没有k个点则输出-1;
2、将x与x父亲的边权加上k。
保证每次操作2的k以及原树的边权小于等于一个数len。
如果操作2中x为1,那么视为将x的基础深度加上了k。
由乃能不能找到出题人呢?

Input

第一行三个数n,m,len。
之后n - 1行每行两个数表示2~n每个点的父亲编号,以及他们到父亲的边权。。。
之后m行每行三个数 opt,x,k,opt表示操作种类,x,k意义如题所述。
n,m <= 100000

Output

对于每个操作1,输出一个数表示答案

Sample Input

3 5 3
1 3
2 3
1 1 3
2 3 3
1 1 3
2 1 2
1 1 3

Sample Output

6
9
11

HINT

对于所有数据,n,m <= 100000,len <= 10
因为出题人是sb,len <= 10,其实没有这个限制也可以解决这个问题

先用dfs序转化成区间问题。区间加,区间第K小。
可以二分答案转化为区间求比x小的数的个数,即每块维护一个有序二元组,表示数值和位置。
然后修改整块加标记,零散块暴力归并,时间复杂度$O(\sqrt n)$,询问对零散块暴力然后整块的二分,时间复杂度$O(\sqrt nlogn)$。
也就是说加上二分答案后时间复杂度为$O(n\sqrt nlognlogn)$,过不去此题,考虑修改块的大小。
设块大小为$size$,修改的复杂度是$O(size+n/size)$,查询的复杂度为$O((size+n/size)logn)$,如果我们把零散块预先合在一起可以优化到$O(size+n/size*logn)$。
算出$size=\sqrt nlogn$时最优,总时间复杂度为$O(n\sqrt nlogn)$。
然而复杂度对了写的丑还是过不去。。。加上register直接快了10秒多,然后因为我写的修改常数太大,块的大小应小一点更好。
手动测试是当块大小为3700时能过。
 
代码:
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
using namespace std;
#define RR register
inline char nc() {
static char buf[100000],*p1,*p2;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int rd() {
RR int x=0;RR char s=nc();
while(s<'0'||s>'9') s=nc();
while(s>='0'&&s<='9') x=(x<<3)+(x<<1)+s-'0',s=nc();
return x;
}
#define N 100500
int head[N],to[N],nxt[N],val[N],cnt,dfn[N],son[N],L[N],R[N],block,dis[N],n,m,pos[N],mx,t[N],more[N],size,fa[N],tot;
struct A {
int v,id;
}a[N],b[N],c[N];
inline bool cmp1(const A &x,const A &y) {
return x.v<y.v;
}
inline void add(int u,int v,int w) {
to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt; val[cnt]=w;
}
void dfs(int x) {
RR int i; dfn[x]=++tot; a[tot].v=dis[x]; a[tot].id=tot;
for(i=head[x];i;i=nxt[i]) {
dis[to[i]]=dis[x]+val[i];
dfs(to[i]);
}
son[x]=tot;
}
inline void update(int l,int r,int v) {
RR int p=pos[l],q=pos[r],i,j,k;
RR int lb=0,lc=0;
if(p==q) {
for(i=L[p];i<=R[p];i++) {
if(a[i].id<=r&&a[i].id>=l) b[++lb]=a[i],b[lb].v+=v;
else c[++lc]=a[i];
}
i=1; j=1; k=L[p];
while(i<=lb&&j<=lc) {
if(b[i].v<=c[j].v) a[k++]=b[i++];
else a[k++]=c[j++];
}
while(i<=lb) a[k++]=b[i++];
while(j<=lc) a[k++]=c[j++];
}else {
for(i=p+1;i<q;i++) more[i]+=v; lb=0,lc=0;
for(i=L[p];i<=R[p];i++) {
if(a[i].id>=l) b[++lb]=a[i],b[lb].v+=v;
else c[++lc]=a[i];
}
i=1; j=1; k=L[p];
while(i<=lb&&j<=lc) {
if(b[i].v<=c[j].v) a[k++]=b[i++];
else a[k++]=c[j++];
}
while(i<=lb) a[k++]=b[i++];
while(j<=lc) a[k++]=c[j++]; lb=0; lc=0;
for(i=L[q];i<=R[q];i++) {
if(a[i].id<=r) b[++lb]=a[i],b[lb].v+=v;
else c[++lc]=a[i];
}
i=1; j=1; k=L[q];
while(i<=lb&&j<=lc) {
if(b[i].v<=c[j].v) a[k++]=b[i++];
else a[k++]=c[j++];
}
while(i<=lb) a[k++]=b[i++];
while(j<=lc) a[k++]=c[j++];
}
}
inline int search(int l,int r,int x) {
RR int tmp=l,p=pos[l]; r++;
x-=more[p];
while(l<r) {
int mid=(l+r)>>1;
if(a[mid].v<=x) l=mid+1;
else r=mid;
}
return l-tmp;
}
inline int findt(int x) {
RR int l=1,r=t[0]+1,mid;
while(l<r) {
mid=(l+r)>>1;
if(t[mid]<=x) l=mid+1;
else r=mid;
}
return l-1;
}
inline int query(int x,int y,int k) {
if(y-x+1<k) return -1;
RR int p=pos[x],q=pos[y],i,j;
if(p==q) {
t[0]=0;
for(i=L[p];i<=R[p];i++) {
if(a[i].id>=x&&a[i].id<=y) t[++t[0]]=a[i].v+more[p];
}
}else {
t[0]=0;
RR int lb=0,lc=0;
for(i=L[p];i<=R[p];i++) {
if(a[i].id>=x) b[++lb]=a[i],b[lb].v+=more[p];
}
for(i=L[q];i<=R[q];i++) {
if(a[i].id<=y) c[++lc]=a[i],c[lc].v+=more[q];
}
i=1; j=1;
while(i<=lb&&j<=lc) {
if(b[i].v<=c[j].v) t[++t[0]]=b[i++].v;
else t[++t[0]]=c[j++].v;
}
while(i<=lb) t[++t[0]]=b[i++].v;
while(j<=lc) t[++t[0]]=c[j++].v;
}
RR int l=0,r=mx+1,mid,re;
while(l<r) {
mid=(l+r)>>1; re=0;
if(t[0]) re+=findt(mid);
for(i=p+1;i<q;i++) re+=search(L[i],R[i],mid);
if(re>=k) r=mid;
else l=mid+1;
}
return l;
}
int main() {
//freopen("a.in","r",stdin);
//freopen("a.out","w",stdout);
RR int kitty,opt;
n=rd(); m=rd(); kitty=rd();
RR int i,x,y,j;
// int ln=n,lg=1;
// while(ln) lg++,ln>>=1;
for(i=2;i<=n;i++) {
x=rd(); y=rd(); add(x,i,y);
}
// size=((int)ceil(sqrt(n)))*(lg);
// printf("%d\n",size);return 0;
size=3700;
dfs(1);
block=n/size;
for(i=1;i<=block;i++) {
L[i]=R[i-1]+1; R[i]=i*size;
sort(a+L[i],a+R[i]+1,cmp1);
mx=max(mx,a[R[i]].v);
for(j=L[i];j<=R[i];j++) pos[j]=i;
}
if(R[block]!=n) {
block++; L[block]=R[block-1]+1; R[block]=n;
sort(a+L[block],a+n+1,cmp1);
for(i=L[block];i<=n;i++) pos[i]=block;
mx=max(mx,a[n].v);
}
while(m--) {
opt=rd(); x=rd(); y=rd();
if(opt==1) {
printf("%d\n",query(dfn[x],son[x],y));
}else {
update(dfn[x],son[x],y); mx+=y;
}
}
return 0;
}

BZOJ_4867_[Ynoi2017]舌尖上的由乃_分块+dfs序的更多相关文章

  1. BZOJ_4765_普通计算姬_分块+dfs序+树状数组

    BZOJ_4765_普通计算姬_分块 Description "奋战三星期,造台计算机".小G响应号召,花了三小时造了台普通计算姬.普通计算姬比普通计算机要厉害一些 .普通计算机能 ...

  2. BZOJ4867 : [Ynoi2017]舌尖上的由乃

    首先通过DFS序将原问题转化为序列上区间加.询问区间kth的问题. 考虑分块,设块大小为$K$,每块维护排序过后的$pair(值,编号)$. 对于修改,整块的部分可以直接打标记,而零碎的两块因为本来有 ...

  3. BZOJ4867 Ynoi2017舌尖上的由乃(dfs序+分块)

    容易想到用dfs序转化为序列上的问题.考虑分块,对每块排序,修改时对于整块打上标记,边界暴力重构排序数组,询问时二分答案,这样k=sqrt(nlogn)时取最优复杂度nsqrt(nlogn)logn, ...

  4. (原创)舌尖上的c++--相逢

    引子 前些时候,我在群里出了一道题目:将变参的类型连接在一起作为字符串并返回出来,要求只用函数实现,不能借助于结构体实现.用结构体来实现比较简单: template<typename... Ar ...

  5. 舌尖上的硬件:CPU/GPU芯片制造解析(高清)(组图)

    一沙一世界,一树一菩提,我们这个世界的深邃全部蕴藏于一个个普通的平凡当中.小小的厨房所容纳的不仅仅是人们对味道的情感,更有推动整个世界前进的动力.要想理解我们的世界,有的时候只需要细细品味一下我们所喜 ...

  6. BZOJ_3252_攻略_线段树+dfs序

    BZOJ_3252_攻略_线段树+dfs序 Description 题目简述:树版[k取方格数] 众所周知,桂木桂马是攻略之神,开启攻略之神模式后,他可以同时攻略k部游戏.今天他得到了一款新游戏< ...

  7. BZOJ_3729_Gty的游戏_博弈论+splay+dfs序

    BZOJ_3729_Gty的游戏_博弈论+splay+dfs序 Description 某一天gty在与他的妹子玩游戏. 妹子提出一个游戏,给定一棵有根树,每个节点有一些石子,每次可以将不多于L的石子 ...

  8. BZOJ_2821_作诗(Poetize)_分块

    BZOJ_2821_作诗(Poetize)_分块 Description 神犇SJY虐完HEOI之后给傻×LYD出了一题:SHY是T国的公主,平时的一大爱好是作诗.由于时间紧迫,SHY作完诗 之后还要 ...

  9. BZOJ_2434_[NOI2011]_阿狸的打字机_(AC自动机+dfs序+树状数组)

    描述 http://www.lydsy.com/JudgeOnline/problem.php?id=2434 给出\(n\)个字符串,\(m\)个询问,对于第\(i\)个询问,求第\(x_i\)个字 ...

随机推荐

  1. 十六进制颜色转换为iOS可以用的UIColor

    // //  UIColor+Transformation.h //  ContactApp // //  Created by 袁冬冬 on 15/9/11. //  Copyright (c) 2 ...

  2. ambari2.6.1汉化记录

    1.1测试机 Apache hadoop2.6Apache ambari 2.6.1集群规模:单节点操作系统 CentOS7以下所有操作均在root用户下执行 1.2安装环境 安装Maventar - ...

  3. 视频客观质量评价工具:MSU Video Quality Measurement Tool【ssim,psnr】

    MSU Video Quality Measurement Tool(msu vqmt)是莫斯科国立大学(Moscow State University)的Graphics and Media Lab ...

  4. Mac 下实现 pyenv/virtualenv 与 Anaconda 的兼容

    http://blog.csdn.net/vencent7/article/details/76849849 自己一直用的 pyenv 和 pyenv-virtualenv 管理不同的 python ...

  5. RA layer request failed

    新整的Eclipse环境出现这个问题,细化内容是不能connect,后来想起切换Eclipse底层库的事情,然后打开Eclipse的SVN设置.把SVN Client借口由JavaHL改为PureJa ...

  6. LR性能测试结果样例分析

    http://www.cnblogs.com/hyzhou/archive/2011/11/16/2251316.html   测试结果分析 LoadRunner性能测试结果分析是个复杂的过程,通常可 ...

  7. git - 简明指南

    助你入门 git 的简明指南,木有高深内容 ;) 作者:罗杰·杜德勒 感谢:@tfnico, @fhd 和 Namics如有纰漏,请在 github 提报问题 安装 下载 git OSX 版 下载 g ...

  8. echart css样式弄了几个月了,样式一直都没有好好解决,好激动

    legend 图例,每个图表最多仅有一个图例.try this » 名称 默认值 描述 {string} orient 'horizontal' 布局方式,默认为水平布局,可选为:'horizonta ...

  9. 写个批处理脚本来帮忙干活--遍历文件夹&字符串处理

    这次打算写几篇关于脚本方面的博客,主要是记录一下 Gradle 脚本和批处理脚本的一些写法,方便后续查阅. 前言 平常开发过程中,一些较为重复的手工性工作,如果能让脚本来帮忙处理,自然是最好的,刚好之 ...

  10. 找jar包的网站 还没用过2017.12.19

    http://www.findjar.com/index.x http://mvnrepository.com/这个比较好用