bzoj3600: 没有人的算术

http://www.lydsy.com/JudgeOnline/problem.php?id=3600

题意大概是维护一个由二叉树组成的序列a，初始每个位置为单个点的树，支持单点修改（把a[l]和a[r]分别作为根的左右子树，放在a[x]），查询区间最大值（按先序遍历字典序比较，空<非空）的位置。

用动态标号法维护一下当前已出现过的所有二叉树的大小顺序，线段树支持区间查询。

#include<bits/stdc++.h>
typedef long long i64;
const int P=1844677,MX=1<<30;
char ib[10000007],*ip=ib;
int _(){
    int x=0;
    while(*ip<48)++ip;
    while(*ip>47)x=x*10+*ip++-48;
    return x;
}
int n,m,tr[277777],a[133333],mx,idp=0;
struct node{
    node*c[2];
    int rnd,M,ls,rs;
    void rb(int L,int R){
        if(!this)return;
        M=L+R;
        c[0]->rb(L,M>>1);
        c[1]->rb(M>>1,R);
    }
}ns[500007],*np=ns,*rt;
bool cmp(int a,int b){return ns[a].M<ns[b].M;}
void up(int&w,int l,int r){w=cmp(a[l],a[r])?r:l;}
bool cmp(node*a,node*b){return a==ns||(a->ls!=b->ls?cmp(a->ls,b->ls):cmp(a->rs,b->rs));}
int h[P][4];
#define $(a,b,c) (*a=b,a=&c,b=*a)
void ins(int x,int ls,int rs){
    node**w=&rt;
    *np=(node){0,0,rand(),0,ls,rs};
    int L=0,R=MX;
    for(;;){
        node*a=*w;
        if(!a||np->rnd>a->rnd){
            *w=np;
            np->M=L+R;
            node**l=np->c,**r=l+1;
            while(a)cmp(a,np)?$(l,a,a->c[1]):$(r,a,a->c[0]);
            *l=*r=0;
            return np++->rb(L,R);
        }
        int d=cmp(a,np);
        w=a->c+d;
        (d?L:R)=(L+R)>>1;
    }
}
int getid(int a,int b){
    int w=(a*149+b*293)%P;
    while(h[w][0]){
        if(h[w][1]==a&&h[w][2]==b)return h[w][3];
        w=(w+13999)%P;
    }
    h[w][0]=1;
    h[w][1]=a;
    h[w][2]=b;
    h[w][3]=++idp;
    ins(idp,a,b);
    return idp;
}
int main(){
    fread(ib,1,sizeof(ib),stdin);
    n=_();m=_();
    srand(n+m+131);
    for(mx=2;mx<=n+2;mx<<=1);
    for(int i=1;i<=n;++i)tr[mx+i]=i;
    for(int i=mx-1;i;--i)tr[i]=tr[i*2];
    ins(0,-1,-1);
    while(m--){
        int o=_(),l=_(),r=_();
        if(o+48=='C'){
            int x=_();
            a[x]=getid(a[l],a[r]);
            for(int w=mx+x>>1;w;w>>=1)up(tr[w],tr[w*2],tr[w*2+1]);
        }else{
            int ml=l,mr=r;
            for(l+=mx,r+=mx;r-l>1;l>>=1,r>>=1){
                if(~l&1)up(ml,ml,tr[l+1]);
                if(r&1)up(mr,tr[r-1],mr);
            }
            up(ml,ml,mr);
            printf("%d
",ml);
        }
    }
    return 0;
}