bzoj4605: 崂山白花蛇草水 权值线段树套KDtree
bzoj4605: 崂山白花蛇草水
链接
思路
强制在线,那就权值线段树套KDtree好了,没啥好讲的。
KDtree要加平衡因子来重构。另外,那水真难喝。
错误
树套树一边写过了,然后是各种傻屌错误。
我居然离散化了权值,要被gzy嘲笑了。
我一开始还笑话那些写动态开点的慢,
然后就被啪啪脸了。
我回收内存不知道咋的多了个if。
然后就逼近内存上限半个G。
总之就是很多傻吊错误。
代码
#include <bits/stdc++.h>
#define cmin(a,b) (a>b?a=b:a)
#define cmax(a,b) (a>b?a:a=b)
using namespace std;
const int N=1e5+7;
const double alpha=0.8;
int read() {
int x=0,f=1;char s=getchar();
for(;s>'9'||s<'0';s=getchar()) if(s=='-') f=-1;
for(;s>='0'&&s<='9';s=getchar()) x=x*10+s-'0';
return x*f;
}
int Wn;
struct point {int x[2];} p[N];
bool cmp(point a,point b) {
return a.x[Wn]<b.x[Wn];
}
struct kdnode {
int ch[2],ma[2],mi[2],siz;
point a;
}t[N*30];
int stak[N*30],top,num;
int newnode() {
if(top) return stak[top--];
return ++num;
}
void update(int x,int y) {
cmin(t[x].mi[0],t[y].mi[0]);
cmin(t[x].mi[1],t[y].mi[1]);
cmax(t[x].ma[0],t[y].ma[0]);
cmax(t[x].ma[1],t[y].ma[1]);
}
void up(int u) {
t[u].siz=t[t[u].ch[0]].siz+t[t[u].ch[1]].siz+1;
t[u].mi[0]=t[u].ma[0]=t[u].a.x[0];
t[u].mi[1]=t[u].ma[1]=t[u].a.x[1];
if(t[u].ch[0]) update(u,t[u].ch[0]);
if(t[u].ch[1]) update(u,t[u].ch[1]);
}
int build(int l,int r,int wn) {
if(l>r) return 0;
int mid=(l+r)>>1;
Wn=wn;
nth_element(p+l,p+mid,p+r+1,cmp);
int u=newnode();
t[u].a=p[mid];
t[u].ch[0]=build(l,mid-1,wn^1);
t[u].ch[1]=build(mid+1,r,wn^1);
up(u);
return u;
}
void pia(int u,int num) {
if(t[u].ch[0]) pia(t[u].ch[0],num);
p[num+t[t[u].ch[0]].siz+1]=t[u].a;
stak[++top]=u;
if(t[u].ch[1]) pia(t[u].ch[1],num+t[t[u].ch[0]].siz+1);
}
void check(int &x,int wn) {
if(t[x].siz*alpha<t[t[x].ch[0]].siz||t[x].siz*alpha<t[t[x].ch[1]].siz)
pia(x,0),x=build(1,t[x].siz,wn);
}
void insert(point a,int &u,int wn) {
if(!u) {
u=newnode();
t[u].a=a;
t[u].ch[0]=t[u].ch[1]=0;
up(u);
return;
}
if(a.x[wn]<t[u].a.x[wn]) insert(a,t[u].ch[0],wn^1);
else insert(a,t[u].ch[1],wn^1);
up(u),check(u,wn);
}
int pd(point a,point b,int k) {
if(a.x[0]<=t[k].mi[0]&&t[k].ma[0]<=b.x[0]&&
a.x[1]<=t[k].mi[1]&&t[k].ma[1]<=b.x[1]) return 2;
if(a.x[0]>t[k].ma[0]||b.x[0]<t[k].mi[0]) return 0;
if(a.x[1]>t[k].ma[1]||b.x[1]<t[k].mi[1]) return 0;
return 1;
}
int query(point a,point b,int u) {
if(!u||!pd(a,b,u)) return 0;
if(pd(a,b,u)==2) return t[u].siz;
int ans=0;
if(a.x[0]<=t[u].a.x[0]&&t[u].a.x[0]<=b.x[0]&&a.x[1]<=t[u].a.x[1]&&t[u].a.x[1]<=b.x[1]) ans++;
ans+=query(a,b,t[u].ch[0]);
ans+=query(a,b,t[u].ch[1]);
return ans;
}
namespace seg {
struct dsr_sb {
int ch[2],rt;
}t[N*30];
int tot;
void insert(int l,int r,int &u,point a,int val) {
if(!u) u=++tot;
insert(a,t[u].rt,0);
if(l==r) return;
int mid=(l+r)>>1;
if(val<=mid) insert(l,mid,t[u].ch[0],a,val);
else insert(mid+1,r,t[u].ch[1],a,val);
}
int ask(int l,int r,int u,point a,point b,int k) {
if(l==r) return l;
int mid=(l+r)>>1,tmp=query(a,b,t[t[u].ch[1]].rt);
if(k>tmp) return ask(l,mid,t[u].ch[0],a,b,k-tmp);
else return ask(mid+1,r,t[u].ch[1],a,b,k);
}
}
int main() {
int n=read(),q=read(),lastans=0,rt=0;
n+=n;
for(int i=1;i<=q;++i) {
int opt=read(),k;
point a,b;
a.x[0]=read()^lastans,a.x[1]=read()^lastans;
if(opt==1) {
k=read()^lastans;
seg::insert(1,1e9,rt,a,k);
} else {
b.x[0]=read()^lastans,b.x[1]=read()^lastans;
k=read()^lastans;
int tmp=query(a,b,seg::t[rt].rt);
if(k>tmp) {
puts("NAIVE!ORZzyz.");
lastans=0;
} else {
lastans=seg::ask(1,1e9,rt,a,b,k);
printf("%d\n",lastans);
}
}
}
return 0;
}
bzoj4605: 崂山白花蛇草水 权值线段树套KDtree的更多相关文章
- 【bzoj4605】崂山白花蛇草水 权值线段树套KD-tree
题目描述 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实力,他轻松地进了山东省省队,现在便是他履行诺言的时候了.蒟蒻Bob ...
- 崂山白花蛇草水 权值线段树套KDtree
Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实 力,他轻松地进了山东省省队,现在便是他履行诺言的时 ...
- 【BZOJ4605】崂山白花蛇草水 权值线段树+kd-tree
[BZOJ4605]崂山白花蛇草水 Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实力,他轻松地进了 ...
- 洛谷P4848 崂山白花蛇草水 权值线段树+KDtree
题目描述 神犇 \(Aleph\) 在 \(SDOI\ Round2\) 前立了一个 \(flag\):如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇 \(Aleph\) 的实力,他轻松地进了山 ...
- [bzoj3196][Tyvj1730]二逼平衡树_树套树_位置线段树套非旋转Treap/树状数组套主席树/权值线段树套位置线段树
二逼平衡树 bzoj-3196 Tyvj-1730 题目大意:请写出一个维护序列的数据结构支持:查询给定权值排名:查询区间k小值:单点修改:查询区间内定值前驱:查询区间内定值后继. 注释:$1\le ...
- [BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树)
[BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树) 题面 原题面有点歧义,不过从样例可以看出来真正的意思 有n个位置,每个位置可以看做一个集合. ...
- BZOJ 3110 ZJOI 2013 K大数查询 树套树(权值线段树套区间线段树)
题目大意:有一些位置.这些位置上能够放若干个数字. 如今有两种操作. 1.在区间l到r上加入一个数字x 2.求出l到r上的第k大的数字是什么 思路:这样的题一看就是树套树,关键是怎么套,怎么写.(话说 ...
- BZOJ3110[Zjoi2013]K大数查询——权值线段树套线段树
题目描述 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c如果是2 a b c形式,表示询问从第a个位置到第b个位置,第C大的数是 ...
- 【bzoj3110】[Zjoi2013]K大数查询 权值线段树套区间线段树
题目描述 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c.如果是2 a b c形式,表示询问从第a个位置到第b个位置,第C大的数 ...
随机推荐
- python基础流程控制
流程控制主要分为三大类: 1.if 判断语句 2.while 循坏语句 3.for 循坏语句 下面以举例说明: if 判断语句: user1 = 'seven' user2 = 'alex' pass ...
- Java自学-I/O File类
Java 的File类,以及常用方法 文件和文件夹都是用File代表 步骤 1 : 创建一个文件对象 使用绝对路径或者相对路径创建File对象 package file; import java.io ...
- webpack 里的 import, exports 实现原理
在使用 webpack 对脚本进行打包, 在开发中, 每个文件中都会使用 import 语句来导入一些功能,又会使用 export 语句导出一些功能,为了研究 import 和 export 原理,研 ...
- java web的转发与重定向
java web的转发与重定向 原创 2012年12月15日 11:51:39 19312 转发和重定向的区别 一.使用方法 转发:request.getRequestDispatcher(" ...
- Golang: 解析JSON数据之一
JSON 作为目前最流行的数据传输格式, 相信每个程序员都跟它打过交道吧.使用 Go 语言时,也不可避免的要操作 JSON 数据,令人惊喜的是,Go 内置了序列化和反序列化 JSON 的功能,今天就来 ...
- 【Spring Cloud】Spring Cloud之Zipkin server搭建以及HTTP收集,分布式服务跟踪(2)
一.搭建步骤 1)新建Spring Boot项目,引入pom坐标 <parent> <groupId>org.springframework.boot</groupId& ...
- springboot+jndi+tomcat配置多数据源
1.在application.properties中,添加jndi配置,如下图 2.新建dataSourceConfig类 3.dataSourceConfig类详细代码,这里只贴出其中一个,多个数据 ...
- springboot java.awt.AWTException: headless environment
修改启动类即可 @SpringBootApplicationpublic class SystemctlApplication { public static void main(String[] a ...
- selenium+python自动化100-centos上搭建selenium启动chrome浏览器headless无界面模式
环境准备 前言 selenium在windows机器上运行,每次会启动界面,运行很不稳定.于是想到用chrome来了的headless无界面模式,确实方便了不少. 为了提高自动化运行的效率和稳定性,于 ...
- version GLIBCXX3.4.21 not defined in file libstdc++.so.6 with link time reference
问题:在运行C++程序时,输入 ./a.out 输出:symbol _ZNSt7__cxx1112basic_stringIcSt11char_traitsIcESaIcEEC1Ev, version ...