https://www.luogu.org/problemnew/show/P4848

我的数据结构大概已经废了。

外层权值线段树内层kdtree,外层线段树上二分答案。

码数据结构一时爽,码完debug火葬场。

要rebuild时少写了个else什么的

插入不upd下去的时候没把值更新完比如sz什么的

比较Int的时候把外层数组套多套少什么的

拿什么拯救你,我的数据结构

吸氧才能过的代码,据说没人能不吸氧过。

//Achen
#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);i++)
#define Rep(i,a,b) for(int i=(a);i>=(b);i--)
#define Formylove return 0
const int N=100007,UP=1e9;
typedef long long LL;
typedef double db;
using namespace std;
int n,q,lastans;
int xl,yl,xr,yr,Qrs; template<typename T> void read(T &x) {
char ch=getchar(); x=0; T f=1;
while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
if(ch=='-') f=-1,ch=getchar();
for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0'; x*=f;
} #define lc ch[x][0]
#define rc ch[x][1]
#define mid ((l+r)>>1) struct pt {
int d[2],v;
}p[N*33];
int rt[N*33],cnt,dt[N*33][2][2],ch[N*33][2],sz[N*33],f[N*33],D; void upd(int x) {
dt[x][0][0]=dt[x][0][1]=p[x].d[0];
dt[x][1][0]=dt[x][1][1]=p[x].d[1];
if(lc) {
dt[x][0][0]=min(dt[x][0][0],dt[lc][0][0]);
dt[x][0][1]=max(dt[x][0][1],dt[lc][0][1]);
dt[x][1][0]=min(dt[x][1][0],dt[lc][1][0]);
dt[x][1][1]=max(dt[x][1][1],dt[lc][1][1]);
}
if(rc) {
dt[x][0][0]=min(dt[x][0][0],dt[rc][0][0]);
dt[x][0][1]=max(dt[x][0][1],dt[rc][0][1]);
dt[x][1][0]=min(dt[x][1][0],dt[rc][1][0]);
dt[x][1][1]=max(dt[x][1][1],dt[rc][1][1]);
}
sz[x]=sz[lc]+sz[rc]+1;
} int sta[N],top;
void tra(int x) {
if(lc) tra(lc);
sta[++top]=x;
if(rc) tra(rc);
} bool cmp(const int &A,const int &B) { return p[A].d[D]<p[B].d[D]; } int build(int l,int r,int dd) {
D=dd;
if(l>r) return 0;
nth_element(sta+l,sta+mid,sta+r+1,cmp);
int x=sta[mid];
lc=build(l,mid-1,dd^1); if(lc) f[lc]=x;
rc=build(mid+1,r,dd^1); if(rc) f[rc]=x;
upd(x); return x;
} int rebuild(int x,int F,int dd) {
top=0; tra(x);
int rs=build(1,top,dd);
f[rs]=F; return rs;
} int bal(int x) { return sz[x]*0.85>=max(sz[lc],sz[rc]); } void insert(int &RT,pt A) {
int x=RT;
for(int fa=0,l=0,dd=0;;dd^=1) {
if(!x) {
x=++cnt; p[x]=A;
upd(x);
f[x]=fa;
if(fa) ch[fa][l]=x;
else RT=x;
break ;
}
dt[x][0][0]=min(dt[x][0][0],A.d[0]);
dt[x][0][1]=max(dt[x][0][1],A.d[0]);
dt[x][1][0]=min(dt[x][1][0],A.d[1]);
dt[x][1][1]=max(dt[x][1][1],A.d[1]); sz[x]++;
if(A.d[dd]<p[x].d[dd]) fa=x,x=lc,l=0;
else fa=x,x=rc,l=1;
}
int ls=-1,dd=0,lsd;
for(int tp=x;;tp=f[tp],dd^=1) {
if(!bal(tp)) ls=tp,lsd=dd;
if(tp==RT) break;
}
if(ls!=-1) {
if(ls==RT) RT=rebuild(ls,0,0);
else ch[f[ls]][ls==ch[f[ls]][1]]=rebuild(ls,f[ls],lsd^dd);
}
} void qry(int x) {
if(!x) return ;
if(dt[x][0][0]>xr||dt[x][0][1]<xl||dt[x][1][0]>yr||dt[x][1][1]<yl) return;
if(dt[x][0][0]>=xl&&dt[x][0][1]<=xr&&dt[x][1][0]>=yl&&dt[x][1][1]<=yr) {
Qrs+=sz[x]; return ;
}
if(p[x].d[0]>=xl&&p[x].d[0]<=xr&&p[x].d[1]>=yl&&p[x].d[1]<=yr) Qrs++;
if(lc) qry(lc); if(rc) qry(rc);
} int sgrt,tot,lson[N*33],rson[N*33];
void upd(int &x,int l,int r,int pos,pt A) {
if(!x) x=++tot; insert(rt[x],A);
if(l==r) return ;
if(pos<=mid) upd(lson[x],l,mid,pos,A);
else upd(rson[x],mid+1,r,pos,A);
} int qry(int x,int l,int r,int k) {
Qrs=0;
if(l==r) {
qry(rt[x]);
if(Qrs>=k) return l;
else return 0;
}
qry(rt[rson[x]]);
if(Qrs>=k) return qry(rson[x],mid+1,r,k);
else return qry(lson[x],l,mid,k-Qrs);
} int main() {
//freopen("4605.in","r",stdin);
//freopen("4605.out","w",stdout);
read(n); read(q);
For(i,1,q) {
int o;
read(o);
if(o==1) {
int x,y,v;
read(x); read(y); read(v);
x^=lastans; y^=lastans; v^=lastans;
upd(sgrt,1,UP,v,(pt){x,y,v});
}
else {
int k;
read(xl); read(yl); read(xr); read(yr);
xl^=lastans; yl^=lastans;
xr^=lastans; yr^=lastans;
read(k); k^=lastans;
lastans=qry(sgrt,1,UP,k);
if(!lastans) puts("NAIVE!ORZzyz.");
else printf("%d\n",lastans);
}
}
Formylove;
}
//Achen
#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);i++)
#define Rep(i,a,b) for(int i=(a);i>=(b);i--)
#define Formylove return 0
const int N=;
typedef long long LL;
typedef double db;
using namespace std;
int a[N][N]; template<typename T> void read(T &x) {
char ch=getchar(); T f=; x=;
while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
if(ch=='-') f=-,ch=getchar();
for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
} int main() {
//freopen("1.in","w",stdout);
srand(time());
int n=,m=rand()%+;
printf("%d %d\n",n,m);
For(i,,m) {
int o;
o=rand()%+;
//if(i<=m/2) o=1; else o=2;
if(i<=) {
int x=rand()%n+,y=rand()%n+,v=rand()%+;
printf("1 %d %d %d\n",x,y,v);
}
else {
int xl,yl,xr,yr,k;
xl=rand()%n+,xr=rand()%n+;
yl=rand()%n+,yr=rand()%n+;
if(xl>xr) swap(xl,xr);
if(yl>yr) swap(yl,yr);
k=rand()%n+;
printf("2 %d %d %d %d %d\n",xl,yl,xr,yr,k);
}
}
Formylove;
}

luogu P4848 崂山白花蛇草水的更多相关文章

  1. 洛谷P4848 崂山白花蛇草水 权值线段树+KDtree

    题目描述 神犇 \(Aleph\) 在 \(SDOI\ Round2\) 前立了一个 \(flag\):如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇 \(Aleph\) 的实力,他轻松地进了山 ...

  2. P4848 崂山白花蛇草水

    题意:支持修改的矩形第 \(k\) 大. 题解:动态开点权值线段树 套 Kd-tree. 然后也没什么难的但就是写不对...调了两天才调出来然后发现跑的巨慢,于是又%了一发Claris'题解,跑的真快 ...

  3. bzoj 4605: 崂山白花蛇草水

    Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实 力,他轻松地进了山东省省队,现在便是他履行诺言的时 ...

  4. LG4848 崂山白花蛇草水

    题意 神犇 Aleph 在 SDOI Round2 前立了一个 flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇 Aleph 的实力,他轻松地进了山东省省队,现在便是他履行诺言的时候了. ...

  5. BZOJ4605:崂山白花蛇草水

    浅谈\(K-D\) \(Tree\):https://www.cnblogs.com/AKMer/p/10387266.html 题目传送门:https://lydsy.com/JudgeOnline ...

  6. 【BZOJ4605】崂山白花蛇草水 权值线段树+kd-tree

    [BZOJ4605]崂山白花蛇草水 Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实力,他轻松地进了 ...

  7. 【bzoj4605】崂山白花蛇草水 权值线段树套KD-tree

    题目描述 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实力,他轻松地进了山东省省队,现在便是他履行诺言的时候了.蒟蒻Bob ...

  8. 崂山白花蛇草水 权值线段树套KDtree

    Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实 力,他轻松地进了山东省省队,现在便是他履行诺言的时 ...

  9. bzoj4605: 崂山白花蛇草水 权值线段树套KDtree

    bzoj4605: 崂山白花蛇草水 链接 bzoj loj 思路 强制在线,那就权值线段树套KDtree好了,没啥好讲的. KDtree要加平衡因子来重构.另外,那水真难喝. 错误 树套树一边写过了, ...

随机推荐

  1. php的异常和处理

    常见错误处理类型 语法错误 环境错误 逻辑错误 常见错误级别 Deprecated 最低级别的错误 不推荐,不建议,使用一些过期函数的时候会出现,程序继续执行 Notice 通知级别的错误 使用一些未 ...

  2. Python学习进程(15)常用内置函数

        本节介绍Python的一些常用的内置函数.     (1)cmp(x, y): cmp()函数比较 x 和 y 两个对象,并根据比较结果返回一个整数,如果 x<y,则返回-1:如果x&g ...

  3. 【leetcode刷题笔记】Text Justification

    Given an array of words and a length L, format the text such that each line has exactly L characters ...

  4. 0731 #Django rest framework

    FBV:Function base viewsdef index(request):    if request.method == 'POST':        return HTTPrespons ...

  5. blast+学习之search tools

    search tools:blastn, blastp, blastx, tblastx, tblastn, psiblast, rpsblast, and rpstblastn 1.blastn: ...

  6. 20145230《Java程序设计》第4周学习总结

    20145230<Java程序设计>第4周学习总结 教材学习内容总结 继承共同行为 本周学习的首先是Java语言中的继承与多态.何为我们的继承呢?在我们面向对象中,子类继承父类,避免重复的 ...

  7. Ubuntu 15.10环境下安装Hive

    1. 安装MySQL sudo tar -xzvf mysql-5.7.13-linux-glibc2.5-x86_64.tar.gz sudo mv mysql-5.7.13-linux-glibc ...

  8. avaweb学习总结(八)——HttpServletResponse对象(二)

    一.HttpServletResponse常见应用——生成验证码 1.1.生成随机图片用作验证码 生成图片主要用到了一个BufferedImage类,

  9. Kubernetes pod网络解析

    在Kubernetes中,会为每一个pod分配一个IP地址,pod内的所有容器都共享这个pod的network namespace,彼此之间使用localhost通信. 那么pod内所有容器间的网络是 ...

  10. XXL-Job高可用集群搭建

    如果XXL-Job admin挂掉就完蛋了,所有任务无法执行 调度中心:管理任务的触发 调度中心如何实现集群? XXL-Job如何实现集群? 底层已经实现好了!文档里面有的 如果想实现Job集群:   ...