BZOJ3196: Tyvj 1730 二逼平衡树
主席树的常数蜜汁优越,在BZOJ上跑了rnk1。
做法很简单,主席树套BIT。
1-3做法很简单,第四个和第五个做法转换成前两个就行了。
//BZOJ 3196 //by Cydiater //2016.12.10 #include <iostream> #include <queue> #include <map> #include <ctime> #include <cmath> #include <cstring> #include <string> #include <iomanip> #include <algorithm> #include <cstdlib> #include <cstdio> #include <bitset> #include <set> #include <vector> using namespace std; #define ll long long #define up(i,j,n) for(int i=j;i<=n;i++) #define down(i,j,n) for(int i=j;i>=n;i--) #define cmax(a,b) a=max(a,b) #define cmin(a,b) a=min(a,b) #define FILE "psh" const int MAXN=5e4+5; const int oo=0x3f3f3f3f; inline int read(){ char ch=getchar();int x=0,f=1; while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } int N,M,root[MAXN<<1],arr[MAXN],fsort[MAXN<<1],rnum,cnt=0,top[2],preL[MAXN],preR[MAXN],tmpL[MAXN],tmpR[MAXN],TMP; struct Chair_man_Tree{ int son[2],sum; }t[MAXN<<6]; struct Query{ int opt,L,R,pos,K; }query[MAXN]; namespace solution{ inline int lowbit(int i){return ((i)&(-i));} int NewNode(int sum,int son0,int son1){ t[++cnt].sum=sum;t[cnt].son[0]=son0;t[cnt].son[1]=son1; return cnt; } void insert(int leftt,int rightt,int &Root,int last,int pos){ Root=NewNode(t[last].sum+1,t[last].son[0],t[last].son[1]); int mid=(leftt+rightt)>>1; if(leftt==rightt) return; if(pos<=mid) insert(leftt,mid,t[Root].son[0],t[last].son[0],pos); else insert(mid+1,rightt,t[Root].son[1],t[last].son[1],pos); } void PushTree(int L,int R){ top[0]=top[1]=0; preL[++top[0]]=root[(L-1==0?0:(N+L-1))]; preR[++top[1]]=root[R+N]; for(int j=L-1;j>=1;j-=lowbit(j)) preL[++top[0]]=root[j]; for(int j=R;j>=1;j-=lowbit(j)) preR[++top[1]]=root[j]; } void ChooseLeft(){ up(i,1,top[0])preL[i]=t[preL[i]].son[0]; up(i,1,top[1])preR[i]=t[preR[i]].son[0]; } void ChooseRight(){ up(i,1,top[0])preL[i]=t[preL[i]].son[1]; up(i,1,top[1])preR[i]=t[preR[i]].son[1]; } int Col(){ int sum=0; up(i,1,top[0])sum-=t[t[preL[i]].son[0]].sum; up(i,1,top[1])sum+=t[t[preR[i]].son[0]].sum; return sum; } void Prepare(){ rnum=N=read();M=read(); up(i,1,N)arr[i]=fsort[i]=read(); up(i,1,M){ query[i].opt=read(); if(query[i].opt==3){ query[i].pos=read();query[i].K=read(); fsort[++rnum]=query[i].K;continue; } query[i].L=read();query[i].R=read();query[i].K=read(); if(query[i].opt==4||query[i].opt==5) fsort[++rnum]=query[i].K; } sort(fsort+1,fsort+rnum+1); rnum=unique(fsort+1,fsort+rnum+1)-(fsort+1); up(i,1,N){ arr[i]=lower_bound(fsort+1,fsort+rnum+1,arr[i])-fsort; insert(1,rnum,root[i+N],root[i+N-1],arr[i]); } } int Rank(int leftt,int rightt,int pos){ if(leftt==rightt){ TMP=0; up(i,1,top[0])TMP-=t[preL[i]].sum; up(i,1,top[1])TMP+=t[preR[i]].sum; return 1; } int sum=0,mid=(leftt+rightt)>>1; sum=Col(); if(pos<=mid){ ChooseLeft(); return Rank(leftt,mid,pos); }else{ ChooseRight(); return sum+Rank(mid+1,rightt,pos); } } int Get(int leftt,int rightt,int rnk){ int sum,mid=(leftt+rightt)>>1; if(leftt==rightt) return fsort[leftt]; sum=Col(); if(rnk<=sum){ ChooseLeft(); return Get(leftt,mid,rnk); }else{ ChooseRight(); return Get(mid+1,rightt,rnk-sum); } } void Insert(int leftt,int rightt,int &Root,int pos,int tag){ if(!Root)Root=NewNode(0,0,0); t[Root].sum+=tag; if(leftt==rightt) return; int mid=(leftt+rightt)>>1; if(pos<=mid) Insert(leftt,mid,t[Root].son[0],pos,tag); else Insert(mid+1,rightt,t[Root].son[1],pos,tag); } void Slove(){ up(i,1,M){ int opt=query[i].opt; if(opt==1){ int L=query[i].L,R=query[i].R,K=lower_bound(fsort+1,fsort+rnum+1,query[i].K)-fsort; PushTree(L,R); printf("%d\n",Rank(1,rnum,K)); } if(opt==2){ int L=query[i].L,R=query[i].R,K=query[i].K; PushTree(L,R); printf("%d\n",Get(1,rnum,K)); } if(opt==3){ int Pos=query[i].pos,K=query[i].K; for(int j=Pos;j<=N;j+=lowbit(j)) Insert(1,rnum,root[j],arr[Pos],-1); arr[Pos]=lower_bound(fsort+1,fsort+rnum+1,K)-fsort; for(int j=Pos;j<=N;j+=lowbit(j)) Insert(1,rnum,root[j],arr[Pos],1); } if(opt==4){ int L=query[i].L,R=query[i].R,K=lower_bound(fsort+1,fsort+rnum+1,query[i].K)-fsort; PushTree(L,R); int rank=Rank(1,rnum,K); PushTree(L,R); printf("%d\n",Get(1,rnum,rank-1)); } if(opt==5){ int L=query[i].L,R=query[i].R,K=lower_bound(fsort+1,fsort+rnum+1,query[i].K)-fsort; PushTree(L,R); int rank=Rank(1,rnum,K); PushTree(L,R); printf("%d\n",Get(1,rnum,rank+TMP)); } } } } int main(){ using namespace solution; Prepare(); Slove(); return 0; }
BZOJ3196: Tyvj 1730 二逼平衡树的更多相关文章
- bzoj3196: Tyvj 1730 二逼平衡树 树套树
地址:http://www.lydsy.com/JudgeOnline/problem.php?id=3196 题目: 3196: Tyvj 1730 二逼平衡树 Time Limit: 10 Sec ...
- 【线段树套平衡树】【pb_ds】bzoj3196 Tyvj 1730 二逼平衡树
线段树套pb_ds里的平衡树,在洛谷OJ上测试,后三个测试点TLE #include<cstdio> #include<algorithm> #include<ext/p ...
- [bzoj3196]Tyvj 1730 二逼平衡树——线段树套平衡树
题目 Description 您需要写一种数据结构(可参考题目标题),来维护一个有序数列,其中需要提供以下操作: 1.查询k在区间内的排名 2.查询区间内排名为k的值 3.修改某一位值上的数值 4.查 ...
- [bzoj3196][Tyvj 1730][二逼平衡树] (线段树套treap)
Description 您需要写一种数据结构(可参考题目标题),来维护一个有序数列,其中需要提供以下操作: 1.查询k在区间内的排名 2.查询区间内排名为k的值 3.修改某一位值上的数值 4.查询k在 ...
- 【分块】bzoj3196 Tyvj 1730 二逼平衡树
分块 或 树套树. 在每个块中维护一个有序表,查询时各种二分,全都是分块的经典操作,就不详细说了. 块的大小定为sqrt(n*log2(n))比较快. #include<cstdio> # ...
- 【带修莫队】【权值分块】bzoj3196 Tyvj 1730 二逼平衡树
这题用了三种算法写: 分块+二分:O(n*sqrt(n*log(n)) 函数式权值分块:O(n*sqrt(n)) 带修莫队+权值分块:O(n5/3) 结果……复杂度越高的实际上跑得越快……最后这个竟然 ...
- 【函数式权值分块】【分块】bzoj3196 Tyvj 1730 二逼平衡树
#include<cstdio> #include<cmath> #include<algorithm> using namespace std; #define ...
- bzoj 3196 Tyvj 1730 二逼平衡树(线段树套名次树)
3196: Tyvj 1730 二逼平衡树 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1807 Solved: 772[Submit][Stat ...
- BZOJ 3196: Tyvj 1730 二逼平衡树( 树套树 )
这道题做法应该很多吧.... 我用了线段树套treap.... -------------------------------------------------------------------- ...
随机推荐
- 将SHP导入MySQL中
ogr2ogr -f MySQL MySQL:smfs,host=127.0.0.1,user=root,password=gis D:\spatialData\HB\HuBeiPicture\HuB ...
- iOS类中的属性设置背景色(统一)
unsigned int count; objc_property_t *properties = class_copyPropertyList([self class], &count); ...
- 关于bundle install 的一点补充
在第一次运行bundle install之后,生成了Gemfile.lock文件,里面记录gem的具体版本号,按照官方文档说明,以后运行bundle install就不会再依据Gemfile,而是根据 ...
- SQL Server 2012 新特性:新增和修改函数
转换函数 1.PARSE Parse是把字符串类型转化为想要的类型,看看和convert和cast的区别 SELECT PARSE ('2.111111' AS f ...
- 在虚拟机中安装CentOS
1.准备工具 我当时下载的是VMware9.0.2,之后升级即可. 2.安装VMware9.0.2,按照步骤安装即可,安装成功并运行 选择创建新的虚拟机,出现下图,选择"自定义"后 ...
- VS 插件ReSharper10 破解注册方法(转)
ReSharper 10.0.0.1 Ultimate 完美破解补丁使用方法,本资源来自互联网,感谢吾乐吧软件站的分享. ReSharper是一款由jetbrains开发的针对C#, VB.NET, ...
- WPF CheckBox样式 ScrollViewer样式 WrapPanel、StackPanel、Grid布局
本节讲述布局,顺带加点样式给大家看看~单纯学布局,肯定是枯燥的~哈哈 那如上界面,该如何设计呢? 1.一些布局元素经常用到.Grid StackPanel Canvas WrapPanel等.如上这种 ...
- 关于HTTP协议,一篇就够了
HTTP简介 HTTP协议是Hyper Text Transfer Protocol(超文本传输协议)的缩写,是用于从万维网(WWW:World Wide Web )服务器传输超文本到本地浏览器的传送 ...
- 【转】虚拟机VMware与主机共享文件介绍
from: http://www.cnblogs.com/kerrycode/p/3818095.html 写的比较详细,但是vm版本较旧. 2:通过共享文件夹功能 虚拟机VMware提供了在宿主机与 ...
- 时隔一年再读到the star
The Star Arthur C. Clarke It is three thousand light-years to the Vatican. Once, I believed that spa ...