5333. 【NOIP2017提高A组模拟8.23】大新闻 主席树/树套树

时间:2022-12-17 14:22:28

数据结构学傻选手。。。
题意:给出一个序列,每次可以在开头添加或删除,以及查询区间k大。

一看添加无脑树套树,成功fst,事实上只要一个主席树就可以了,倒过来每次覆盖掉。。

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
int n,m,tot,top,sz;
const int N=5e5+5;
const int maxn=2e5;
int v[N],num[N],hash[N];
int flag[N],A[N],B[N],K[N],root[N];
int sum[N*60],ls[N*60],rs[N*60],sta,end;
int L[100],R[100],a,b;
int lowbit(int x)
{return x&(-x);}
int find(int x)
{
int l=1,r=tot,mid;
while(l<=r)
{
int mid=(l+r)>>1;
if(hash[mid]<x)l=mid+1;
else r=mid-1;
}
return l;
}
void updata(int last,int l,int r,int &rt,int w,int x)
{
rt=++sz;
sum[rt]=sum[last]+x;ls[rt]=ls[last];rs[rt]=rs[last];
if(l==r) return;
int mid=(l+r)>>1;
if(w<=mid) updata(ls[last],l,mid,ls[rt],w,x);
else updata(rs[last],mid+1,r,rs[rt],w,x);
}
int query(int l,int r,int k)
{
if(l==r) return l;
int i,suml=0,sumr=0;
for(i=1;i<=a;i++) suml+=sum[ls[L[i]]];
for(i=1;i<=b;i++) sumr+=sum[ls[R[i]]];
int mid=(l+r)>>1;
if(sumr-suml>=k)
{
for(i=1;i<=a;i++) L[i]=ls[L[i]];
for(i=1;i<=b;i++) R[i]=ls[R[i]];
return query(l,mid,k);
}
else
{
for(i=1;i<=a;i++) L[i]=rs[L[i]];
for(i=1;i<=b;i++) R[i]=rs[R[i]];
return query(mid+1,r,k-(sumr-suml));
}
}
int main()
{
freopen("news.in","r",stdin);
freopen("news.out","w",stdout);
char s[3];
scanf("%d%d",&n,&m);
sta=1+maxn,end=n+maxn;
for(int i=sta;i<=end;i++)
{

scanf("%d",&v[i]);
num[++top]=v[i];
}
for(int i=1;i<=m;i++)
{
int op;
scanf("%d",&op);
if (op==1)A[i]=-1;
else if (op==2)scanf("%d",&B[i]),num[++top]=B[i];
else {scanf("%d%d%d",&A[i],&B[i],&K[i]);flag[i]=1;}
}
sort(num+1,num+top+1);
hash[++tot]=num[1];
for(int i=2;i<=top;i++)
if(num[i]!=num[i-1])
hash[++tot]=num[i];
for(int i=sta;i<=end;i++)
{
int t=find(v[i]);
for(int j=i;j<=end;j+=lowbit(j))
updata(root[j],1,tot,root[j],t,1);
}
for(int i=1;i<=m;i++)
if(flag[i])
{
A[i]+=sta-1;
B[i]+=sta-1;
a=0;b=0;A[i]--;
for(int j=A[i];j>=sta;j-=lowbit(j))
L[++a]=root[j];
for(int j=B[i];j>=sta;j-=lowbit(j))
R[++b]=root[j];
printf("%d\n",hash[query(1,tot,K[i])]);
}
else
{
if (A[i]==-1)
{
sta++;
continue;
}
sta--;int t;
if (v[sta])
{
t=find(v[sta]);
for(int j=sta;j<=end;j+=lowbit(j))
updata(root[j],1,tot,root[j],t,-1);
}
v[sta]=B[i];
t=find(B[i]);
for(int j=sta;j<=end;j+=lowbit(j))
updata(root[j],1,tot,root[j],t,1);
}

return 0;
}

主席树:

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

const int N=3e5+5;
const int inf=1e9;

int n,m,sz,a[N],root[N*2];
struct tree{int s,l,r;}t[N*2*30];

void ins(int &d,int p,int l,int r,int x)
{
d=++sz;t[d]=t[p];t[d].s++;
if (l==r) return;
int mid=(l+r)/2;
if (x<=mid) ins(t[d].l,t[p].l,l,mid,x);
else ins(t[d].r,t[p].r,mid+1,r,x);
}

int query(int d,int p,int l,int r,int k)
{
if (l==r) return l;
int mid=(l+r)/2;
if (t[t[d].l].s-t[t[p].l].s>=k) return query(t[d].l,t[p].l,l,mid,k);
else return query(t[d].r,t[p].r,mid+1,r,k-t[t[d].l].s+t[t[p].l].s);
}

int main()
{
freopen("news.in","r",stdin);
freopen("news.out","w",stdout);
scanf("%d%d",&n,&m);
for (int i=n;i>=1;i--) scanf("%d",&a[i]);
for (int i=1;i<=n;i++) ins(root[i],root[i-1],1,inf,a[i]);
while (m--)
{
int op;
scanf("%d",&op);
if (op==1) n--;
else if (op==2)
{
int x;
scanf("%d",&x);
n++;
ins(root[n],root[n-1],1,inf,x);
}
else
{
int l,r,k;
scanf("%d%d%d",&l,&r,&k);
printf("%d\n",query(root[n-l+1],root[n-r],1,inf,k));
}
}
return 0;
}