【题解】SDOI-2008 郁闷的小J

话说这题bzoj上没有,丢一个洛谷的连接 点我

这题不是很难,对于每一种书建立一个有序的位置数组,就是这种书出现的位置集合,然后就是维护这个序列有序了

蒟蒻只能想到两种:set和平衡树,set就二分查找,SBT就find

但这题最变态的一点是卡空间!!!

蒟蒻的SBT结构体在数据极端情况下欲哭无泪
把里面和外面的数组改来改去都A不了

于是蒟蒻使出了浑身解数卡空间:

队列回收空间

vector动态数组

剪除冗杂

而且是在结构体里边外边都回收空间,开上vector

最后发现直接用一个大数组维护森林即可

改来改去的代码:

#include<bits/stdc++.h>
using namespace std;

#define rg register
#define cl(x) memset(x,0,sizeof(x))
#define clm(x) memset(x,0x7fffffff,sizeof(x))
#define oo 0x7fffffff

#define L ch[p][0]
#define R ch[p][1]
#define LL ch[ch[p][0]][0]
#define LR ch[ch[p][0]][1]
#define RL ch[ch[p][1]][0]
#define RR ch[ch[p][1]][1]

template <typename _Tp> inline void read(_Tp&x){
    char c11=getchar();x=0;bool booo=0;
    while(c11!='-'&&!isdigit(c11))c11=getchar();if(c11=='-'){c11=getchar();booo=1;}
    while(isdigit(c11)){x=x*10+c11-'0';c11=getchar();}if(booo)x=-x;return ;
}

const int limit1=100001;
const int limit2=50;

const int N=500000;
int sz;
int data[N],ch[N][2],size[N],rt[N];
void up(int p){
        size[p]=size[L]+size[R]+1;
        return ;
    }
    void clear(){
        cl(data),cl(ch);
        return ;
    }
    void zig(int&p){
        int k=ch[p][1];
        ch[p][1]=ch[k][0];
        ch[k][0]=p;
        up(p);up(k);
        p=k;
        return ;
    }
    void zag(int&p){
        int k=ch[p][0];
        ch[p][0]=ch[k][1];
        ch[k][1]=p;
        up(p);up(k);
        p=k;
        return ;
    }
    void maintain(int&p,bool flag){
        if(!flag)
            if(size[LL]>size[R])zag(p);
            else
                if(size[LR]>size[R]){zig(L);zag(p);}
                else return ;
        else
            if(size[RR]>size[L])zig(p);
            else
                if(size[RL]>size[L]){zag(R);zig(p);}
                else return ;
        maintain(L,0);maintain(R,1);
        maintain(p,0);maintain(p,1);
        return ;
    }
    void ins(int&p,int x){
        if(!p){p=++sz;data[p]=data[sz]=x,size[sz]=1;return ;}
        if(data[p]==x)return ;
        ++size[p];
        ins(ch[p][data[p]<x],x);
        up(p);
        maintain(p,x>=data[p]);
        return ;
    }
    int erase(int &p,int x){
        --size[p];int tmp;
        if(x==data[p]||(x<data[p]&&!ch[p][0])||(data[p]<x&&!ch[p][1])){
            tmp=data[p];
            if(!ch[p][0]||!ch[p][1])p=ch[p][0]+ch[p][1];
            else data[p]=erase(ch[p][0],data[p]+1);
            return tmp;
        }
        tmp=erase(ch[p][data[p]<=x],x);
        return tmp;
    }
    int rank(int t,int x){
        int p=t,tmp=0;
        while(1){
            if(data[p]==x)return tmp+size[L]+1;
            if(x<data[p])
                if(L)p=L;
                else return tmp;
            else{
                tmp+=size[L]+1;
                if(R)p=R;
                else return tmp;
            }
        }
        return tmp;
    }
    int build(int l,int r){
        if(l>r)return 0;
        int mid=(l+r)>>1;
        int p=++sz;
        L=build(l,mid-1);
        read(data[p]);size[p]=1;
        R=build(mid+1,r);
        up(p);
        return p;
    }

int a[N],z=0;
map <int,int> mp;
queue <int> q;

int main(){
    freopen("depressedj.in","r",stdin);
    freopen("depressedj.out","w",stdout);
    int n,m,tmp;
    read(n);read(m);
    for(rg int i=1;i<=n;++i){
        read(a[i]);
        tmp=mp[a[i]];
        if(!tmp){mp[a[i]]=++z,tmp=z;}
        ins(rt[tmp],i);
    }
    char c[3];int l,r,x;
    while(m--){
        scanf("%s",c);
        if(c[0]=='Q'){
            read(l);read(r);read(x);
            int p=mp[x];
            if(!p||rt[p]==0){puts("0");continue;}
            printf("%d\n",rank(rt[p],r)-rank(rt[p],l-1));
        }
        else {
            read(l);read(r);
            int p1=mp[a[l]],p2=mp[r];
            if(!p2)mp[r]=p2=++z;
            erase(rt[p1],l);
            ins(rt[p2],l);a[l]=r;
        }
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值