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

数据结构学傻选手。。。
题意：给出一个序列，每次可以在开头添加或删除，以及查询区间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;
}