BZOJ_4864_[BeiJing 2017 Wc]神秘物质_Splay

BZOJ4864_[BeiJing 2017 Wc]神秘物质_Splay

Description

21ZZ 年，冬。

小诚退休以后， 不知为何重新燃起了对物理学的兴趣。 他从研究所借了些实验仪器，整天研究各种微观粒子。这

一天， 小诚刚从研究所得到了一块奇异的陨石样本， 便迫不及待地开始观测。 在精密仪器的视野下，构成陨石

的每个原子都无比清晰。 小诚发现， 这些原子排成若干列， 每一列的结构具有高度相似性。于是，他决定对单

独一列原子进行测量和测试。被选中的这列共有 N 个顺序排列的原子。 最初， 第 i 个原子具有能量 Ei。 随着

时间推移和人为测试， 这列原子在观测上会产生两种变化：

merge x e 当前第 x 个原子和第 x+1 个原子合并，得到能量为 e 的新原子；

insert x e 在当前第 x 个原子和第 x+1 个原子之间插入一个能量为 e 的新原子。

对于一列原子，小诚关心的是相邻一段中能量最大和能量最小的两个原子的能量差值，

称为区间极差。 因此， 除了观测变化外，小诚还要经常统计这列原子的两类数据：

max x y 当前第 x 到第 y 个原子之间的任意子区间中区间极差的最大值；

min x y 当前第 x 到第 y 个原子之间的任意子区间中区间极差的最小值。

其中， 子区间指的是长度至少是 2 的子区间。

小诚坚信这项研究可以获得诺贝尔物理学奖。为了让小诚早日了结心愿，你能否帮助他实现上述的观测和测量呢？

Input

第一行， 两个整数 N， M， 分别表示最初的原子数目和事件总数。

第二行， N 个整数 E1, E2, …, EN， 由空格隔开。依次表示每个原子的能量。

接下来 M 行， 每行为一个字符串和两个整数， 描述一次事件，格式见题目描述。

N<=100,000,M<=100,000

1 ≤ e, Ei ≤ 109。 设 N’ 为当前时刻原子数目。

对于 merge 类事件， 1 ≤ x ≤ N’-1；

对于 insert 类事件， 1 ≤ x ≤ N’；

对于 max 和 min 类事件， 1 ≤ x < y ≤ N’。

任何时刻，保证 N’ ≥ 2。

Output

输出若干行， 按顺序依次表示每次 max 和 min 类事件的测量结果。

Sample Input

4 3
5 8 10 2
max 1 3
min 1 3
max 2 4

Sample Output

5 2 8

 

---

根据极差的定义，我们可以发现，一个子序列越长极差越大，子序列越短极差越小。

于是所求的最大值一定是整个区间的极差，最小值一定是区间内所有相邻两个数的极差的最小值。

由于操作有插入和删除，我们用Splay维护。

需要维护区间最小值，区间最大值，以及这个节点所代表的区间左、右端点，这三个都是可以由儿子推出来的。

 

代码：

#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define N 200050
#define ls ch[p][0]
#define rs ch[p][1]
#define get(x) (ch[f[x]][1]==x)
int ch[N][2],f[N],siz[N],mn[N],mx[N],rt,sz,mxrange[N],mnrange[N],val[N],tot;
int n,m,a[N],head[N],tail[N];
char opt[10];
int fabs(int x){return x>0?x:-x;}
void pushup(int p) {
    siz[p]=siz[ls]+siz[rs]+1;
    mx[p]=max(val[p],max(mx[ls],mx[rs]));
    mn[p]=min(val[p],min(mn[ls],mn[rs]));
    head[p]=ls?head[ls]:p;
    tail[p]=rs?tail[rs]:p;
    mnrange[p]=min(mnrange[ls],mnrange[rs]);
    if(ls) mnrange[p]=min(mnrange[p],fabs(val[p]-val[tail[ls]]));
    if(rs) mnrange[p]=min(mnrange[p],fabs(val[p]-val[head[rs]]));
}
void rotate(int x) {
    int y=f[x],z=f[y],k=get(x);
    ch[y][k]=ch[x][!k]; f[ch[y][k]]=y;
    ch[x][!k]=y; f[y]=x; f[x]=z;
    if(z) ch[z][ch[z][1]==y]=x;
    pushup(y); pushup(x);
    if(rt==y) rt=x;
}
void splay(int x,int y) {
    for(int fa;(fa=f[x])!=y;rotate(x))
        if(f[fa]!=y)
            rotate(get(fa)==get(x)?fa:x);
}
int find(int x) {
    int p=rt;
    while(1) {
        if(x<=siz[ls]) p=ls;
        else {
            x-=siz[ls]+1;
            if(!x) return p;
            p=rs;
        }
    }
}
void build(int fa,int l,int r) {
    if(l>r) return ;
    int mid=(l+r)>>1;
    ch[fa][mid>fa]=mid;
    f[mid]=fa;
    val[mid]=a[mid-1];
    build(mid,l,mid-1);
    build(mid,mid+1,r);
    pushup(mid);
}
void print() {
    int i;
    printf("tot=%d\n",tot);
    for(i=1;i<=tot;i++) {
        int p=find(i);
        printf("p=%d,val[p]=%d,siz[p]=%d\n",p,val[p],siz[p]);
    }
}
int main() {
    scanf("%d%d",&n,&m);
    int i,x,y,p;
    mx[0]=0; mn[0]=mnrange[0]=0x3fffffff;
    for(i=1;i<=n;i++) scanf("%d",&a[i]);
    build(0,1,n+2);
    rt=(n+3)>>1; sz=n+2; tot=n+2;
    while(m--) {
        scanf("%s%d%d",opt,&x,&y);
        //print();
        if(opt[1]=='e') {
            tot--;
            p=x+3; x=find(x); p=find(p); splay(x,0); splay(p,x);
            ch[ls][0]=ch[ls][1]=0; val[ls]=y;
            pushup(ls); pushup(p); pushup(x);
        }else if(opt[0]=='i') {
            tot++;
            x++; p=x+1; x=find(x); p=find(p); splay(x,0); splay(p,x);
            sz++; ls=sz; f[sz]=p; val[sz]=y;
            pushup(sz); pushup(p); pushup(x);
        }else if(opt[1]=='a') {
            x=find(x); p=find(y+2); splay(x,0); splay(p,x);
            printf("%d\n",mx[ls]-mn[ls]);
        }else {
            x=find(x); p=find(y+2); splay(x,0); splay(p,x);
            printf("%d\n",mnrange[ls]);
        }
    }
}