崂山白花蛇草水 权值线段树套KDtree
崂山白花蛇草水 权值线段树套KDtree
Description
Input
Output
Code:
// luogu-judger-enable-o2
// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define setIO(s) freopen(s".in","r",stdin) //,freopen(s".out","w",stdout)
struct Data
{int ch[2],w,minv[2],maxv[2],p[2],siz,sum;
}node[20000000], arr[20000000];
namespace KDtree
{#define maxn 20000000#define Min(a,b) (a = a>b?b:a) #define Max(a,b) (a = b>a?b:a) int d,c_arr; std::queue<int>Q; void Init(){for(int i=1;i<maxn;++i) Q.push(i); } int newnode(){ int q = Q.front(); Q.pop(); return q; }bool cmp(Data i,Data j){ return i.p[d]==j.p[d]?i.p[d^1]<j.p[d^1]:i.p[d]<j.p[d]; }bool isout(int k,int x1,int y1,int x2,int y2){if(node[k].maxv[0]<x1||node[k].minv[0]>x2||node[k].maxv[1]<y1||node[k].minv[1]>y2) return 1; return 0; }bool isin(int k,int x1,int y1,int x2,int y2){if(node[k].maxv[0]<=x2 && node[k].minv[0]>=x1&&node[k].maxv[1]<=y2&&node[k].minv[1]>=y1) return 1; return 0; }void pushup(int x,Data p){node[x].sum+=p.sum;node[x].siz+=p.siz; Min(node[x].minv[0],p.minv[0]); Min(node[x].minv[1],p.minv[1]); Max(node[x].maxv[0],p.maxv[0]); Max(node[x].maxv[1],p.maxv[1]); }int query(int x,int x1,int y1,int x2,int y2){if(!x||isout(x,x1,y1,x2,y2)) return 0; if(isin(x,x1,y1,x2,y2)) {return node[x].sum; }int ans=0; if(node[x].p[0]>=x1&&node[x].p[0]<=x2&&node[x].p[1]>=y1&&node[x].p[1]<=y2) ans=1; ans += query(node[x].ch[0],x1,y1,x2,y2)+query(node[x].ch[1],x1,y1,x2,y2); return ans; } void dfs(int &o){if(!o) return; dfs(node[o].ch[0]);arr[++c_arr] = node[o];arr[c_arr].sum=arr[c_arr].w=1; arr[c_arr].maxv[0]=arr[c_arr].minv[0]=arr[c_arr].p[0]; arr[c_arr].maxv[1]=arr[c_arr].minv[1]=arr[c_arr].p[1]; arr[c_arr].siz=1; Q.push(o);dfs(node[o].ch[1]); o=0; } void rebuild(int &x,int l,int r,int o){int mid=(l+r)>>1; d=o,std::nth_element(arr+l,arr+mid,arr+1+r,cmp); x=newnode(); node[x]=arr[mid]; node[x].minv[0]=node[x].maxv[0]=node[x].p[0]; node[x].minv[1]=node[x].maxv[1]=node[x].p[1]; node[x].sum=node[x].w=1; node[x].ch[0]=node[x].ch[1]=0;node[x].siz=1; if(l<mid) rebuild(node[x].ch[0],l,mid-1,o^1), pushup(x, node[node[x].ch[0]]); if(r>mid) rebuild(node[x].ch[1],mid+1,r,o^1), pushup(x, node[node[x].ch[1]]); } void Reconstruct(int &u){c_arr=0; dfs(u); rebuild(u,1,c_arr,d=0); } bool check(int son,int x){ if(son*10>x*9) return true; return false; }void insert(int &x,Data a,int de,int tag){if(!x){x=newnode(); node[x].p[0]=node[x].maxv[0]=node[x].minv[0]=a.p[0]; node[x].p[1]=node[x].maxv[1]=node[x].minv[1]=a.p[1]; node[x].sum=node[x].w=1;node[x].siz=1,node[x].ch[0]=node[x].ch[1]=0; return; }d = de;int is=0; if(cmp(a,node[x]) == 1) { is = check(node[node[x].ch[0]].siz+1,node[x].siz+1);insert(node[x].ch[0],a,de^1,is||tag);pushup(x,a); }//a is bigger (rson) else {is = check(node[node[x].ch[1]].siz+1,node[x].siz+1);insert(node[x].ch[1],a,de^1,is||tag);pushup(x,a); } if(tag && is) Reconstruct(x); }
};
#define y1 opopopop
int lson[2000000],rson[2000000],rt[2000000];
int x1,y1,x2,y2, tot=0,root=0;
void modify(int &x,int l,int r,int k,Data i)
{ if(!x) x = ++tot; KDtree::insert(rt[x],i,0,0); if(l==r) return; int mid = (l+r)>>1; if(k <= mid) modify(lson[x],l,mid,k,i); else modify(rson[x],mid+1,r,k,i);
}
int query(int x,int l,int r,int kth)
{ if(l==r) return l; int u = KDtree::query(rt[rson[x]],x1,y1,x2,y2); int mid = (l + r) >> 1; if(u < kth) return query(lson[x],l,mid,kth-u); else return query(rson[x],mid+1,r,kth);
}
int Query(int kth)
{ return query(root,0,1000000000,kth);
}
int n,q,lastans=0;
int main()
{// setIO("input"); KDtree::Init(); scanf("%d%d",&n,&q); for(int cc=1;cc<=q;++cc) {int opt,a,b,c,d; scanf("%d",&opt); if(opt==1){ scanf("%d%d%d",&a,&b,&c); a^=lastans,b^=lastans,c^=lastans; Data f; f.p[0]=f.maxv[0]=f.minv[0]=a;f.p[1]=f.maxv[1]=f.minv[1]=b; f.ch[0]=f.ch[1]=0; f.siz=1; f.sum=f.w=1; modify(root,0,1000000000,c,f); }if(opt==2){scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&a); x1^=lastans,y1^=lastans,x2^=lastans,y2^=lastans,a^=lastans; lastans = Query(a); if(lastans==0) printf("NAIVE!ORZzyz.\n"); else printf("%d\n",lastans); }}return 0;
}
崂山白花蛇草水 权值线段树套KDtree相关推荐
- bzoj4605: 崂山白花蛇草水 权值线段树套KDtree
bzoj4605: 崂山白花蛇草水 链接 bzoj loj 思路 强制在线,那就权值线段树套KDtree好了,没啥好讲的. KDtree要加平衡因子来重构.另外,那水真难喝. 错误 树套树一边写过了, ...
- 【bzoj4605】崂山白花蛇草水 权值线段树套KD-tree
题目描述 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实力,他轻松地进了山东省省队,现在便是他履行诺言的时候了.蒟蒻Bob ...
- bzoj4605 崂山白花蛇草水 权值线段树套kd树
Description Q次操作,要求资瓷 在(x,y)处放一个数字x 查询(x1,y1)到(x2,y2)矩形内第k大 Solution 非常裸的权值线段树套kd树,为了保证复杂度可以定期重构也可以平 ...
- 【BZOJ4605】崂山白花蛇草水 权值线段树+kd-tree
[BZOJ4605]崂山白花蛇草水 Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实力,他轻松地进了 ...
- BZOJ 4605 崂山白花蛇草水 权值线段树+K-D树
Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实 力,他轻松地进了山东省省队,现在便是他履行诺言的时 ...
- BZOJ 4605: 崂山白花蛇草水 树套树 权值线段树套kdtree
4605: 崂山白花蛇草水 Time Limit: 80 Sec Memory Limit: 512 MB Submit: 527 Solved: 153 [Submit][Status][Dis ...
- bzoj4605: 崂山白花蛇草水(权值线段树套k-dtree)
题目 #include<bits/stdc++.h> using namespace std; #define mid (l+r>>1) const int N=100002, ...
- 刷题总结——序列操作(权值线段树套树状数组)
题目: 题目描述 给出序列 a1,a2,-,an(0≤ai≤109),有关序列的两种操作. 1. ai(1≤i≤n)变成 x(0≤x≤109). 2. 求 al,al+1,-,ar(1≤l≤r≤n)第 ...
- 【BZOJ3110】【codevs1616】K大数查询,权值线段树套普通线段树
Time:2016.05.09 Author:xiaoyimi 转载注明出处谢谢 传送门1 传送门2 思路: 之前没怎么接触过权值线段树(非主席树),这次就当学习了一下吧.一开始还把题意理解错了,我的 ...
最新文章
- getservbyname和getservbyport
- 万能makefile深入浅出 - 第三篇
- [js高手之路]从原型链开始图解继承到组合继承的产生
- Linux网络编程服务器模型选择之IO复用循环并发服务器
- 编译内核_将驱动编译进内核(Kernel)的步骤记录
- 【渝粤教育】国家开放大学2018年秋季 0680-22T会计基础知识 参考试题
- c 语言从大到小排序算法,10 大经典排序算法(动图演示+ C 语言代码)
- c++可达矩阵_测试分析人员必备知识—需求管理和可追溯性矩阵
- python 学习之 WINDOWS 注册表学习(一)
- yml配置mysql自动重连_数据库和Doctrine(转载自http://www111cnnet/phper/332/85987htm) - Doctrine...
- i2c-tools 使用集锦
- Office下载记录
- SQL按每5分钟或者10分钟分组统计数据
- 后勤事务繁杂低效?三步骤解决企业行政管理难题
- 墨天轮国产数据库沙龙 | 胡彦军:华为GaussDB迁移工具解密
- 百分点科技启动上市:消息称其国内业务持续亏损,苏萌为实控人
- java获取手机通讯录权限_android获取手机通讯录
- Java中角度和弧度的转换,三角函数,反三角函数
- php外边框样式,CSS的边框样式
- 《HelloGitHub》第 49 期