luogu-P1110(treap)

题目描述
Q的妈妈是一个出纳，经常需要做一些统计报表的工作。今天是妈妈的生日，小Q希望可以帮妈妈分担一些工作，作为她的生日礼物之一。

经过仔细观察，小Q发现统计一张报表实际上是维护一个非负整数数列，并且进行一些查询操作。

在最开始的时候，有一个长度为N的整数序列，并且有以下三种操作：

INSERT i k：在原数列的第i个元素后面添加一个新元素kk；如果原数列的第ii个元素已经添加了若干元素，则添加在这些元素的最后（见下面的例子）
MIN_GAP：查询相邻两个元素的之间差值（绝对值）的最小值
MIN_SORT_GAP：查询所有元素中最接近的两个元素的差值（绝对值）
例如一开始的序列为5, 3, 15,3,1。

执行操作INSERT 2 9将得到：5, 3, 9, 15,3,9,1，此时MIN_GAP为22，MIN_SORT_GAP为22。

再执行操作INSERT 2 6将得到：5, 3, 9, 6, 15,3,9,6,1

注意这个时候原序列的第22个元素后面已经添加了一个99，此时添加的66应加在99的后面。这个时候MIN_GAP为22，MIN_SORT_GAP为11。

于是小Q写了一个程序，使得程序可以自动完成这些操作，但是他发现对于一些大的报表他的程序运行得很慢，你能帮助他改进程序么？

输入格式
第一行包含两个整数NN，MM，分别表示原数列的长度以及操作的次数。

第二行为NN个整数，为初始序列。

接下来的MM行每行一个操作，即INSERT i k，MIN_GAP，MIN_SORT_GAP 中的一种（无多余空格或者空行）。

输出格式
对于每一个MIN_GAP和MIN_SORT_GAP命令，输出一行答案即可。

输入输出样例
输入

3 5
5 3 1
INSERT 2 9
MIN_SORT_GAP
INSERT 2 6
MIN_GAP
MIN_SORT_GAP

输出

2
2
1

说明/提示
对于30%的数据，N ≤ 1000N≤1000 ,M ≤ 5000,M≤5000
对于100%的数据，N ,M ≤500000
对于所有的数据，序列内的整数不超过$5\times 10^8$
时限2s

这个题我们可以通过建立两棵平衡树来维护，对于MIN_SORT_GAP的答案，我们将每个数都加入平衡树，在加入时，我们找到加入的val的前驱和后继，而减数和被减数有val的减式中最小的值就是这两种情况之一，我们取最小值就好啦，而对于MIN_GAP的答案，我们需要用链表来模拟添加的情况，然后在平衡树中减去原有的相邻的减值，然后在加上加入这个数之后相邻的数相减的情况得到的值。
代码：

#pragma GCC optimize(2)
#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5 + 7;
std::mt19937 rnd(233);
struct treap
{
    int size, val, l, r, key;
} fhq1[maxn<<5], fhq2[maxn<<5];
int a[maxn << 4], b[maxn << 4], c[maxn << 4];
int cnt1, cnt2, root1, root2;
int ans = 1e9 + 10;
inline int newnode(treap *t,int val,int &cnt)
{
    t[++cnt] = treap{1, val, 0, 0, rnd()};
    return cnt;
}
inline void update(treap *t,int cnt)
{
    t[cnt].size = t[t[cnt].l].size + t[t[cnt].r].size + 1;
}
inline void split(treap *t,int now,int val,int &x,int &y)
{
    if(!now){
        x = y = 0;
        return;
    }
    if(t[now].val<=val){
        x = now;
        split(t, t[now].r, val, t[now].r, y);
    }
    else{
        y = now;
        split(t, t[now].l, val, x, t[now].l);
    }
    update(t, now);
}
inline int merge(treap*t,int x,int y)
{
    if(!x||!y){
        return x + y;
    }
    if(t[x].key>=t[y].key){
        t[x].r = merge(t, t[x].r, y);
        update(t, x);
        return x;
    }
    else{
        t[y].l = merge(t, x, t[y].l);
        update(t, y);
        return y;
    }
}
int x, y, z;
inline int pre(treap *t,int root)
{
    while(t[root].l)
        root = t[root].l;
    return t[root].val;
}
inline int nex(treap *t,int root)
{
    while(t[root].r)
        root = t[root].r;
    return t[root].val;
}
inline void insert(treap *t,int &root,int val,int &cnt)
{
    split(t, root, val, x, y);
    if(t==fhq1){
        if(t[x].size){
            ans = min(ans, abs(val - nex(t, x)));
        }
        if(t[y].size){
            ans = min(ans, abs(val - pre(t, y)));
        }
    }
    root = merge(t,merge(t,x, newnode(t, val, cnt)), y);
}
inline void del(treap *t,int &root,int val)
{
    split(t, root, val, x, z);
    split(t, x, val - 1, x, y);
    y = merge(t, t[y].l, t[y].r);
    root = merge(t, merge(t, x, y), z);
}
inline void read(int &val)
{
    int f = 1;
    val = 0;
    char c = getchar();
    while(c<'0' || c>'9'){
        if(c=='-')
            f = -1;
        c = getchar();
    }
    while(c>='0' && c<='9'){
        val = val * 10 + c - '0';
        c = getchar();
    }
    val *= f;
}
int len;
char arr[100];
int main()
{
    int n, m;
    read(n), read(m);
    for (register int i = 1; i <= n;i++){
        read(a[i]);
        insert(fhq1, root1, a[i], cnt1);
    }
    for (register int i = 2; i <= n;i++){
        insert(fhq2, root2, abs(a[i] - a[i - 1]), cnt2);
        b[i] = i - 1;
    }
    b[n + 1] = n;
    for (register int i = 1; i < n;i++){
        c[i] = i + 1;
    }
    len = n;
    while(m--)
    {
        int num, val;
        scanf("%s", arr);
        if(arr[0]=='I'){
            read(num), read(val);
            if(num<n){
                del(fhq2, root2, abs(a[num + 1] - a[b[num + 1]]));
            }
            insert(fhq2, root2, abs(val - a[b[num + 1]]), cnt2);
            if(num<n)
                insert(fhq2, root2, abs(val-a[num+1]), cnt2);
            a[++len] = val;
            b[len] = b[num + 1];
            c[len] = num + 1;
            c[b[num + 1]] = len, b[num + 1] = len;
            insert(fhq1, root1, val, cnt1);
        }
        else{
            if(arr[4]=='G')
                printf("%d\n", pre(fhq2, root2));
            else
                printf("%d\n", ans);
        }
    }
}

最后提一句，这道题会卡常，最后两个点一直过不去，一直T着，最后没办法，气急败坏之下开了氧气飘过[雾]。