bzoj 3600: 没有人的算术 替罪羊树+线段树

题意

题目

分析

vfk出的题，强的一逼，这里是题解。
关于替罪羊树的介绍可以去看clj13年的集训队论文《重量平衡树和后缀平衡树在信息学奥赛中的应用》。
关于这题的话，大概就是开一棵平衡树来维护每一个数的相对大小。因为要给每个点赋一个rank值来O(1)比较大小，而带旋转的平衡树在旋转操作后需要花费大量时间来维护rank，所以这里采用通过暴力重构来维护平衡的替罪羊树。
具体维护rank的过程大概就是因为组成新数的两个数都已经出现过，所以可以很方便地比较大小。便可以直接找到其在平衡树中的对应位置。rank值的话，我们可以给每个节点赋一个区间(l,r)，根节点的区间就是(0,1)，那么每个节点的rank值就是(l+r)/2，其左儿子的区间就是(l,(l+r)/2)，右儿子的区间为((l+r)/2,r)。
为了避免精度问题，可以把区间变为long long。
然后开一棵线段树来维护最大值即可。

代码

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

typedef long long LL;

const int N=500005;
const LL inf=(LL)4611686018427387904;

int n,m,pos[N],rt,top,stack[N],R;
LL val[N],w1[N],w2[N];
struct data{int l,r;};
struct tree{int mx;}t[N*5];

bool operator < (data x,data y)
{
    return val[x.l]<val[y.l]||val[x.l]==val[y.l]&&val[x.r]<val[y.r];
}

bool operator == (data x,data y)
{
    return x.l==y.l&&x.r==y.r;
}

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

struct Scapegoat
{
    int sz;
    struct tree{int l,r,s,fa;}t[N];
    data bel[N];

    void dfs(int x)
    {
        if (t[x].l) dfs(t[x].l);
        stack[++top]=x;
        if (t[x].r) dfs(t[x].r);
    }

    void build(int &d,int l,int r,LL p,LL q)
    {
        if (l>r) {d=0;return;}
        LL mid=(p+q)/2;int m=(l+r)/2;d=stack[m];val[d]=mid;w1[d]=p;w2[d]=q;
        build(t[d].l,l,m-1,p,mid);if (t[d].l) t[t[d].l].fa=d;
        build(t[d].r,m+1,r,mid,q);if (t[d].r) t[t[d].r].fa=d;
        t[d].s=t[t[d].l].s+t[t[d].r].s+1;
    }

    void rebuild(int d,LL p,LL q)
    {
        top=0;int tmp=t[d].fa,op=(d==t[tmp].l?0:1);
        dfs(d);
        build(d,1,top,p,q);
        t[d].fa=tmp;
        if (!tmp) {rt=d;return;}
        if (!op) t[tmp].l=d;
        else t[tmp].r=d;
    }

    int ins(int &d,LL p,LL q,data w)
    {
        if (!d) {d=++sz;t[d].s=1;val[d]=(p+q)/2;w1[d]=p;w2[d]=q;bel[d]=w;return d;}
        if (bel[d]==w) return d;
        LL mid=(p+q)/2;int ans;
        if (w<bel[d]) ans=ins(t[d].l,p,mid,w),t[t[d].l].fa=d;
        else ans=ins(t[d].r,mid,q,w),t[t[d].r].fa=d;
        t[d].s=t[t[d].l].s+t[t[d].r].s+1;
        if (t[d].s*0.75<(max(t[t[d].l].s,t[t[d].r].s)+3)) R=d;
        return ans;
    }
}sc;

void modify(int d,int l,int r,int x)
{
    if (l==r) {t[d].mx=l;return;}
    int mid=(l+r)/2;
    if (x<=mid) modify(d*2,l,mid,x);
    else modify(d*2+1,mid+1,r,x);
    if (val[pos[t[d*2].mx]]>=val[pos[t[d*2+1].mx]]) t[d].mx=t[d*2].mx;
    else t[d].mx=t[d*2+1].mx;
}

int query(int d,int l,int r,int x,int y)
{
    if (x>y) return 0;
    if (l==x&&r==y) return t[d].mx;
    int mid=(l+r)/2;
    int u=query(d*2,l,mid,x,min(y,mid)),v=query(d*2+1,mid+1,r,max(x,mid+1),y);
    if (val[pos[u]]>=val[pos[v]]) return u;
    else return v;
}

int main()
{
    n=read();m=read();val[0]=-inf-1;
    sc.ins(rt,-inf,inf,(data){0,0});
    for (int i=1;i<=n;i++) pos[i]=1,modify(1,1,n,i);
    while (m--)
    {
        char ch[2];scanf("%s",ch);
        int l=read(),r=read();
        if (ch[0]=='C')
        {
            int k=read();
            pos[k]=sc.ins(rt,-inf,inf,(data){pos[l],pos[r]});
            if (R) sc.rebuild(R,w1[R],w2[R]),R=0;
            modify(1,1,n,k);
        }
        else printf("%d\n",query(1,1,n,l,r));
    }
    return 0;
}