替罪羊树套Trie,Trie合并用线段树合并,注意常数优化。

顺便AC800题纪念~~~

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=200010,inf=19,M=32000010;
struct info{
int v1,v2;
info(){v1=v2=0;}
info(int _v1,int _v2){v1=_v1,v2=_v2;}
}maxv;
inline info merge(info a,info b){
if(b.v1>a.v1){
a.v2=b.v2>a.v1?b.v2:a.v1;
a.v1=b.v1;
}else if(b.v1>a.v2)a.v2=b.v1;
return a;
}
struct node{int l,r,v;node(){}node(int _l,int _r,int _v){l=_l,r=_r,v=_v;}}T[M];
int rub,getn[M];
const double A=0.8;
int size[N],real[N],ex[N],son[N][2],val[N],f[N],tot,root,id[N],cnt;
info mv[N];
int h[N],pool[N];
int ans,cntpool,cntone,one[N];
int n,q,x,y,i,Pow[20],Inv[20];char ch;
//Trie begin
inline void Max(int&a,int b){if(a<b)a=b;}
inline int Newtrienode(int l,int r,int v){
int x=getn[rub--];
T[x]=node(l,r,v);
return x;
}
void Ins(int&x,int a,int c,int p){
if(!x)x=Newtrienode(0,0,0);
T[x].v+=p;
if(~a)Ins(c&Pow[a]?T[x].r:T[x].l,a-1,c,p);
}
int Merge(int x,int y){
if(!x&&!y)return 0;
return Newtrienode(Merge(T[x].l,T[y].l),Merge(T[x].r,T[y].r),T[x].v+T[y].v);
}
inline int Ask(int x,int c){
int t=1048575,a;
for(a=inf;t>ans&&~a;a--)if(c&Pow[a]){
if(T[T[x].l].v)x=T[x].l;else t&=Inv[a],x=T[x].r;
}else{
if(T[T[x].r].v)x=T[x].r;else t&=Inv[a],x=T[x].l;
}
return t;
}
inline void Ask(int c){
int i;
for(i=1;i<=cntone;i++)Max(ans,one[i]^c);
for(i=1;i<=cntpool;i++)Max(ans,Ask(pool[i],c));
}
void deltree(int x){
if(!x)return;
deltree(T[x].l);deltree(T[x].r);
getn[++rub]=x;
}
//Trie end
//Scapegoat begin
inline void up(int x){
mv[x]=info(ex[x]?val[x]:0,0);
if(son[x][0])mv[x]=merge(mv[x],mv[son[x][0]]);
if(son[x][1])mv[x]=merge(mv[x],mv[son[x][1]]);
}
inline int newnode(int p,int fa){
int x=++tot;
f[x]=fa;size[x]=real[x]=ex[x]=1;son[x][0]=son[x][1]=0;
val[x]=p;
mv[x]=info(p,0);
Ins(h[x]=0,inf,p,1);
return x;
}
inline int Newnode(int x,int fa){
f[x]=fa;size[x]=1;real[x]=ex[x];son[x][0]=son[x][1]=0;
mv[x]=info(ex[x]?val[x]:0,0);
h[x]=0;
return x;
}
int ins(int x,int p,int b){
size[x]++;real[x]++;mv[x]=merge(mv[x],info(p,0));
Ins(h[x],inf,p,1);
if(!son[x][b])return son[x][b]=newnode(p,x);else return ins(son[x][b],p,b);
}
void dfs(int x){
if(son[x][0])dfs(son[x][0]);
id[++cnt]=x;
if(son[x][1])dfs(son[x][1]);
}
int build(int fa,int l,int r){
int mid=(l+r)>>1,x=Newnode(id[mid],fa);
if(l<mid)size[x]+=size[son[x][0]=build(x,l,mid-1)],real[x]+=real[son[x][0]];
if(r>mid)size[x]+=size[son[x][1]=build(x,mid+1,r)],real[x]+=real[son[x][1]];
h[x]=Merge(h[son[x][0]],h[son[x][1]]);
if(ex[x])Ins(h[x],inf,val[x],1);
up(x);
return x;
}
inline int rebuild(int x){
cnt=0;dfs(x);
for(int i=1;i<=cnt;i++)deltree(h[id[i]]);
return build(f[x],1,cnt);
}
inline int kth(int k,int p){
int x=root,rank;
while(1){
size[x]++;real[x]++;mv[x]=merge(mv[x],info(p,0));Ins(h[x],inf,p,1);
rank=real[son[x][0]]+1;
if(k==rank&&ex[x])return x;
if(k<rank)x=son[x][0];else k-=rank+ex[x]-1,x=son[x][1];
}
}
inline void kthchange(int k,int p){
int x=root,rank,del;
while(1){
rank=real[son[x][0]]+1;
if(k==rank&&ex[x])break;
if(k<rank)x=son[x][0];else k-=rank+ex[x]-1,x=son[x][1];
}
for(del=val[x],val[x]=p;x;x=f[x])Ins(h[x],inf,del,-1),Ins(h[x],inf,p,1),up(x);
}
inline void kthdel(int k){
int x=root,rank,del;
while(1){
rank=real[son[x][0]]+1;
if(k==rank&&ex[x])break;
if(k<rank)x=son[x][0];else k-=rank+ex[x]-1,x=son[x][1];
}
for(del=val[x],ex[x]=0;x;x=f[x])real[x]--,Ins(h[x],inf,del,-1),up(x);
}
inline void kthins(int k,int p){
if(!root){root=newnode(p,0);return;}
int x;
if(k==1)x=ins(root,p,0);
else if(k>n)x=ins(root,p,1);
else{
x=kth(k,p);
if(son[x][0])x=ins(son[x][0],p,1);else{
son[x][0]=newnode(p,x);
x=son[x][0];
}
}
int z=x,deep=0;while(f[z])z=f[z],deep++;
if(deep<log(tot)/log(1/A))return;
while((double)size[son[x][0]]<A*size[x]&&(double)size[son[x][1]]<A*size[x])x=f[x];
if(!x)return;
if(x==root){root=rebuild(x);return;}
int y=f[x],b=son[y][1]==x,now=rebuild(x);
son[y][b]=now;
}
inline void ask(int x,int a,int b,int c,int d){
if(!x)return;
if(c<=a&&b<=d){maxv=merge(maxv,mv[x]);if(h[x])pool[++cntpool]=h[x];return;}
int mid=a+real[son[x][0]];
if(ex[x]&&c<=mid&&mid<=d)maxv=merge(maxv,info(val[x],0)),one[++cntone]=val[x];
if(c<mid)ask(son[x][0],a,mid-1,c,d);
if(d>=mid+ex[x])ask(son[x][1],mid+ex[x],b,c,d);
}
//Scapegoat end
inline void read(int&a){
char ch;while(!(((ch=getchar())>='0')&&(ch<='9')));
a=ch-'0';while(((ch=getchar())>='0')&&(ch<='9'))(a*=10)+=ch-'0';
}
int main(){
for(Pow[0]=i=1;i<=inf;i++)Pow[i]=Pow[i-1]<<1;
for(i=0;i<=inf;i++)Inv[i]=Pow[i]^1048575;
for(i=1;i<M;i++)getn[++rub]=i;
read(n);read(q);
for(i=1;i<=n;i++)read(val[i]),id[i]=i,ex[i]=1;
root=build(0,1,tot=n);
while(q--){
while(!(((ch=getchar())=='D')||(ch=='F')||(ch=='C')||(ch=='I')));
read(x);
if(ch!='D')read(y);
if(ch=='F'){
x=(x+ans)%n,y=(y+ans)%n;
maxv=info();cntpool=cntone=0;
ask(root,1,n,x+1,y+1);
ans=0;
Ask(maxv.v2);
printf("%d\n",ans);
}else if(ch=='C')x=(x+ans)%n,y=(y+ans)&1048575,kthchange(x+1,y);
else if(ch=='I')x=(x+ans)%n,y=(y+ans)&1048575,kthins(x+1,y),n++;
else x=(x+ans)%n,kthdel(x+1),n--;
}
return 0;
}

  

BZOJ3217 : ALOEXT的更多相关文章

  1. 【BZOJ3217】ALOEXT 替罪羊树+Trie树

    [BZOJ3217]ALOEXT Description taorunz平时最喜欢的东西就是可移动存储器了……只要看到别人的可移动存储器,他总是用尽一切办法把它里面的东西弄到手. 突然有一天,taor ...

  2. 【BZOJ3217】ALOEXT 分块

    题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3217 分块过掉辣!!!!!!$O(n^{1.5}+q\times \sqrt{n})$的 ...

  3. 【bzoj3217】ALOEXT 替罪羊树套Trie树

    题目描述 taorunz平时最喜欢的东西就是可移动存储器了……只要看到别人的可移动存储器,他总是用尽一切办法把它里面的东西弄到手. 突然有一天,taorunz来到了一个密室,里面放着一排可移动存储器, ...

  4. bzoj 3217: ALOEXT

    将此神题作为博客园的第一篇文章,至此,数据结构基本学完了(或者说数据结构轮流虐了我一次!) 人生第一道7K代码题! 没什么,就是treap套个trie,然后tle是因为一定要用指针当时p党谁会用那么丑 ...

  5. BZOJ 3217: ALOEXT (块状链表套trie)

    第一次写块状链表,发现还挺好写的,但是一点地方写错加上强制在线就会各种姿势WA/TLE/RE爆- 想法就是分块后,在每一个块上维护最大值和次大值,还在每一个块上维护一棵trie树来求异或最大值.散块直 ...

  6. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

随机推荐

  1. UITapGestureRecognizer

    UITapGestureRecognizer IOS的手势非常多, 但是特别容易其他视图起冲突的手势,要数UITapGestureRecognizer 于是有了gestureRecognizerSho ...

  2. Linux常用热键(持续更新)

    (这些文章都是从我的个人主页上粘贴过来的,大家也可以访问我的主页 www.iwangzheng.com) --圣诞节怎么过, --略过. 今天装ubuntu的时候把windows覆盖了, 凌乱,TX童 ...

  3. [BZOJ1177][Apio2009]Oil

    [BZOJ1177][Apio2009]Oil 试题描述 采油区域 Siruseri政府决定将石油资源丰富的Navalur省的土地拍卖给私人承包商以建立油井.被拍卖的整块土地为一个矩形区域,被划分为M ...

  4. django 架构点点滴滴

    前言: 零星发现一些,零星记录一些,因此可能整体比较混乱,因为显然不是一气呵成写的. 关于CBV(Class Based View): 首先吐槽下,cbv的整体继承结构,可真的不是很优美,可以查看这里 ...

  5. django-cms 代码研究(八)app hooks

    app钩子,啥玩意呢? 就是把现有的app,集成到cms的一种手段. 有两种实现方式: 1) 定义cms_app.py,如下: from cms.app_base import CMSApp from ...

  6. TCP/IP详解学习笔记(5)-- ICMP:internet 控制报文协议

    1.概述      ICMP是(Internet Control Message Protocol)Internet控制报文协议.它是TCP/IP协议族的一个子协议,用于在IP主机.路由器之间传递控制 ...

  7. iOS tableView 选中某个cell时 标准的处理方法

    以前选中cell时,常常判断选中的行数,但是当cell的顺序发生变化时,就要改动处理函数,特别是行数比较多的时候,很麻烦. 之后运用cell的title的内容判断,但是这种判断与现实的内容密切相关,如 ...

  8. Java for LeetCode 066 Plus One

    Given a non-negative number represented as an array of digits, plus one to the number. The digits ar ...

  9. 14.python笔记之paramiko

    作者:刘耀 博客:www.liuyao.me 博客园:www.cnblogs.com/liu-yao 转载请注明 一,介绍 1.使用paramiko可以很好的解决以上问题,比起前面的方法,它仅需要在本 ...

  10. Gym 100851E Easy Problemset (模拟题)

    Problem E. Easy ProblemsetInput file: easy.in Output file: easy.outPerhaps one of the hardest problems ...