【BZOJ3217】ALOEXT 分块

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3217

分块过掉辣!!!!!!$O(n^{1.5}+q\times \sqrt{n})$的分块过掉辣!!

而且速度贼快内存贼小啊(成功踩到b站第一)

由于此题需要强制在线的删除或者插入,所以我们基于块状链表分块,在每个块内存储指定区间内的所有数,以及该区间内的最大值和次大值,同时再维护一个由该区间内所有数组成的trie树。

对于修改一个数,首先在该块的trie数中删除该数(直接伪删),然后再插入,接着更新最大值和次大值。

对于插入一个数,直接在trie树中插入该数,随后在块状链表中插入,更新最大值和次大值。最后判断是否要将该块拆分(如果要察费直接暴力重构两个块即可)

对于删除一个数,trie树伪删除+块链删除+更新最大值和次大值。随后判断该块大小是否等于0(为零可以直接删掉这个块了)

对于查询最大值,首先查询出该区间内的次大值。对于完全包含的块,直接将该值丢到该块的trie树中求最大值,对于非完全包含的块,直接暴力求即可。

PS:经过测试,每个块的大小设置为3200最好(达到3200后切成两块1600的)。最初块大小为400,b站上跑了15.5s,把块大小改为3200后变成3.8s,直接rank1了。

UPD:前段时间被卡到第二,后来重新把块大小调整为2620,并且加上fread与fwrite,卡到3312ms,重回rank1。

(不知能否继续优化.....)

 #include<bits/stdc++.h>
#define M 5000000
#define N 100
using namespace std;
struct trie{
int a[],sum;
}p[M]; int uset=;
int add(int now,int x,int px){
if(!now) now=++uset;
int ret=now;
for(int i=;i>=;i--){
bool k=x&(<<i); p[now].sum+=px;
if(p[now].a[k]) now=p[now].a[k];
else uset++,p[now].a[k]=uset,now=uset;
}
p[now].sum+=px;
return ret;
}
int getmax(int now,int x){
int ans=;
for(int i=;i>=;i--){
bool k=x&(<<i);
if(p[p[now].a[-k]].sum)
now=p[now].a[-k],ans+=(<<i);
else now=p[now].a[k];
}
return ans;
}
struct kuai{
int rt,use,next;
int max1,max2,a[N*+];
kuai(){
rt=use=next=max1=max2=;
memset(a,,sizeof(a));
}
void getmax(){
max1=max2=;
for(int i=;i<=use;i++)
if(a[i]>max1) max2=max1,max1=a[i];
else if(a[i]>max2) max2=a[i];
}
}a[N*]; int use=;
void updata(int id,int zhi){
for(int i=;i;i=a[i].next)
if(a[i].use>=id){
int last=a[i].a[id]; a[i].a[id]=zhi;
a[i].rt=add(a[i].rt,last,-);
add(a[i].rt,zhi,);
a[i].getmax();
return;
}else id-=a[i].use;
}
void cut(int i){
int j=++use;
a[i].use=a[j].use=N*;
a[j].next=a[i].next; a[i].next=j;
for(int k=;k<=N*;k++)
a[j].a[k]=a[i].a[k+N*],a[i].a[k+N*]=;
a[i].max1=a[i].max2=;
a[i].getmax(); a[j].getmax();
a[i].rt=add(,a[i].a[],);
a[j].rt=add(,a[j].a[],);
for(int k=;k<=N*;k++){
add(a[i].rt,a[i].a[k],);
add(a[j].rt,a[j].a[k],);
}
}
void insert(int x,int id){
int i,pre;
for(i=;i;pre=i,i=a[i].next)
if(a[i].use<id) id-=a[i].use;
else break;
if(!i) i=pre,id+=a[i].use;
for(int j=a[i].use;j>=id;j--) a[i].a[j+]=a[i].a[j];
a[i].a[id]=x; a[i].use++;
a[i].rt=add(a[i].rt,x,);
if(x>a[i].max1) a[i].max2=a[i].max1,a[i].max1=x;
else if(x>a[i].max2) a[i].max2=x;
if(a[i].use>=N*) cut(i);
}
void del(int id){
int i,pre=,x;
for(i=;i;pre=i,i=a[i].next)
if(a[i].use<id) id-=a[i].use;
else break;
x=a[i].a[id]; add(a[i].rt,x,-);
for(int j=id;j<=a[i].use;j++) a[i].a[j]=a[i].a[j+];
a[i].use--;
if(a[i].use==) {a[pre].next=a[i].next; return;}
if(a[i].max2>x) return;
a[i].getmax();
}
int getmax2(int l,int r){
int max1=,max2=,sum=;
int i,pre;
for(i=;i;pre=i,i=a[i].next){
if(l<=sum+&&sum+a[i].use<=r){
if(a[i].max1>max1) max2=max1,max1=a[i].max1;
else if(a[i].max1>max2) max2=a[i].max1;
if(a[i].max2>max2) max2=a[i].max2;
}
else if(sum+<=l&&r<=sum+a[i].use){
for(int j=l-sum;j<=r-sum;j++)
if(a[i].a[j]>max1) max2=max1,max1=a[i].a[j];
else if(a[i].a[j]>max2) max2=a[i].a[j];
}
else if(sum+<=l&&l<=sum+a[i].use){
for(int j=l-sum;j<=a[i].use;j++)
if(a[i].a[j]>max1) max2=max1,max1=a[i].a[j];
else if(a[i].a[j]>max2) max2=a[i].a[j];
}
else if(sum+<=r&&r<=sum+a[i].use){
for(int j=;j<=r-sum;j++)
if(a[i].a[j]>max1) max2=max1,max1=a[i].a[j];
else if(a[i].a[j]>max2) max2=a[i].a[j];
}
sum+=a[i].use;
}
return max2;
}
int query(int l,int r){
int x=getmax2(l,r),maxn=,sum=;
int i,pre;
for(i=;i;pre=i,i=a[i].next){
if(l<=sum+&&sum+a[i].use<=r){
maxn=max(maxn,getmax(a[i].rt,x));
}
else if(sum+<=l&&r<=sum+a[i].use){
for(int j=l-sum;j<=r-sum;j++)
maxn=max(maxn,x^a[i].a[j]);
}
else if(sum+<=l&&l<=sum+a[i].use){
for(int j=l-sum;j<=a[i].use;j++)
maxn=max(maxn,x^a[i].a[j]);
}
else if(sum+<=r&&r<=sum+a[i].use){
for(int j=;j<=r-sum;j++)
maxn=max(maxn,x^a[i].a[j]);
}
sum+=a[i].use;
}
return maxn;
}
int main(){
int n,m,lastans=,n0;
scanf("%d%d",&n,&m); n0=n;
for(int i=;i<=n;i++){
int x; scanf("%d",&x);
insert(x,i);
}
while(m--){
char c[]; int x,y;
scanf("%s%d",c,&x);
x=(x+lastans)%n0+;
if(c[]=='D'){del(x); n0--;continue;}
scanf("%d",&y);
if(c[]=='F') y=(y+lastans)%n0+;
else y=(y+lastans)%;
if(c[]=='I') insert(y,x),n0++;
if(c[]=='C') updata(x,y);
if(c[]=='F') printf("%d\n",lastans=query(min(x,y),max(x,y)));
}
}
上一篇:LaTeX语法笔记


下一篇:JSTL标签之core标签的使用