bzoj1588: [HNOI2002]营业额统计 平衡树 Splay

bzoj1588: [HNOI2002]营业额统计

Description

营业额统计 Tiger最近被公司升任为营业部经理，他上任后接受公司交给的第一项任务便是统计并分析公司成立以来的营业情况。 Tiger拿出了公司的账本，账本上记录了公司成立以来每天的营业额。分析营业情况是一项相当复杂的工作。由于节假日，大减价或者是其他情况的时候，营业额会出现一定的波动，当然一定的波动是能够接受的，但是在某些时候营业额突变得很高或是很低，这就证明公司此时的经营状况出现了问题。经济管理学上定义了一种最小波动值来衡量这种情况： 该天的最小波动值 当最小波动值越大时，就说明营业情况越不稳定。 而分析整个公司的从成立到现在营业情况是否稳定，只需要把每一天的最小波动值加起来就可以了。你的任务就是编写一个程序帮助Tiger来计算这一个值。 第一天的最小波动值为第一天的营业额。该天的最小波动值=min{|该天以前某一天的营业额-该天营业额|}。(这句bzoj上没有)

Input

第一行为正整数 ，表示该公司从成立一直到现在的天数，接下来的n行每行有一个整数(有可能有负数) ，表示第i
天公司的营业额。
天数n<=32767,
每天的营业额ai <= 1,000,000。
最后结果T<=2^31

Output

输出文件仅有一个正整数，即Sigma(每天最小的波动值) 。结果小于2^31 。

Sample Input

6
5
1
2
5
4
6

Sample Output

12

题解

仍然是一道平衡树裸题
就是找前缀后缀就好了
好像有双向链表的解法，而且跑得飞快，我不管我不管我就写平衡树
然而平衡树时间贼慢，luogu上没过
最后发现是我自己Splay写水了，栈爆了
我真 tmd 弱！！
注意一下吧，我之前写前驱后继的时候没有Splay在luogu上T/RE到死。
所以记住：
路见不平就Splay！！
各种操作都Splay！！
看它不爽/很爽都Splay！！
重要的事情说三遍。

代码

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#include<cmath>
using namespace std;
const int N = 1100000, inf = 0x3f3f3f3f;
int read() {
    char ch = getchar(); int x = 0, f = 1;
    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;
}
int fa[N], ch[N][2], val[N], root, sz, n, ans;
bool wh(int p) {return ch[fa[p]][1] == p;}
void Rotate(int p) {
    int f = fa[p], g = fa[f], c = wh(p);
    if(g) ch[g][wh(f)] = p; fa[p] = g;
    ch[f][c] = ch[p][c ^ 1]; if(ch[f][c]) fa[ch[f][c]] = f;
    ch[p][c ^ 1] = f; fa[f] = p;
}

void Splay(int p) {
    for(;fa[p]; Rotate(p))
        if(fa[fa[p]]) Rotate(wh(p) == wh(fa[p]) ? fa[p] : p);
    root = p;
}

int find_pre(int u) {
    int ret = 0, p, f;
    for(p = root; p; f = p, p = ch[p][val[p] <= u])
        if(val[p] <= u) ret = u - val[p];
    Splay(f);//这里一定要Splay！
    return ret;
}

int find_nxt(int u) {
    int ret = 0, p, f;
    for(p = root; p; f = p, p = ch[p][val[p] <= u])
        if(val[p] >= u) ret = val[p] - u;
    Splay(f);//这里一定要Splay！
    return ret;
}

void Ins(int u) {
    int p = root;
    for(; ch[p][u > val[p]] && u != val[p]; p = ch[p][u > val[p]]) ;
    if(u == val[p]) return; 
    fa[ch[p][u > val[p]] = ++sz] = p;
    val[ch[p][u > val[p]]] = u; 
    Splay(sz);
}

int main() {
    n = read(); root = 1; val[1] = inf;
    ch[1][0] = 2; fa[2] = 1; val[2] = -inf; sz = 2;
    for(int i = 1;i <= n; ++i) {
        int u = read();
        if(i == 1) Ins(u), ans += u;
        else ans += min(find_pre(u), find_nxt(u)), Ins(u);
    }
    printf("%d\n", ans);
    return 0;
}