BZOJ 1588 HNOI2002 营业额统计

1588: [HNOI2002]营业额统计

Time Limit: 5 Sec  Memory Limit: 162 MB
Submit: 17427  Solved: 7104
[Submit][Status][Discuss]
Description

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

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

HINT

结果说明：5+|1-5|+|2-1|+|5-5|+|4-5|+|6-5|=5+4+1+0+1+1=12


该题数据bug已修复.----2016.5.15

刚开始看道题的时候看不懂

这道题目的大意就是对于每一个i找到i之前离a[i]最接近的数，他们俩的差值作为波动，最后累加

维护区间查询离某个值最近当然Treep更好用一些，getprev和getnext函数查询大和小最近的

一比较即可

代码如下：

#include <bits/stdc++.h>
#define ll long long 
using namespace std;
inline int read(){
    int x=0;int f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int MAXN=1e6+10;
namespace zhangenming{
    int n,a[MAXN],root=0,cnt=0;
    struct tree{
        int leftt,rightt,v,fix;
    }T[MAXN<<1];
    inline void left_rota(int &root){
        int k=T[root].rightt;
        T[root].rightt=T[k].leftt;
        T[k].leftt=root;
        root=k;
    }
    inline void right_rota(int &root){
        int k=T[root].leftt;
        T[root].leftt=T[k].rightt;
        T[k].rightt=root;
        root=k;
    }
    inline void myinsert(int &root,int x){
        //cout<<root<<' '<<x<<endl; 
        if(root==0){
            root=++cnt;
            T[root].v=x;
            T[root].fix=rand()&1007;
            return;
        }
        if(x==T[root].v){
            return;
        }
        else{
            if(x>T[root].v){
                myinsert(T[root].rightt,x);
                if(T[T[root].rightt].fix>T[root].fix) left_rota(root); 
            }
            else{
                myinsert(T[root].leftt,x);
                if(T[T[root].leftt].fix>T[root].fix) right_rota(root);
            }
        }
    }
    inline int getnext(int root,int x){
        if(root==0) return 100000000;
        if(T[root].v<x) return getnext(T[root].rightt,x);
        else {
            return min(T[root].v,getnext(T[root].leftt,x));
        }
    }
    inline int getprev(int root,int x){
        if(root==0) return -1000000;
        if(T[root].v>x){
            return getprev(T[root].leftt,x);
        }
        else return max(T[root].v,getprev(T[root].rightt,x));
    }
    void init(){
        n=read();
        for(int i=1;i<=n;i++){
            a[i]=read();
        }
    }
    void solve(){
        ll sum=a[1];
        memset(T,0,sizeof(T));
        myinsert(root,a[1]);
        for(int i=2;i<=n;i++){
            int k=getprev(root,a[i]);
            //cout<<k<<' ';
            int l=getnext(root,a[i]);
            //cout<<l<<endl;
            sum+=min(abs(k-a[i]),abs(l-a[i]));
            //cout<<sum<<endl;
            myinsert(root,a[i]);
        }
        
        cout<<sum<<endl;
    }
}
int main(){
    //freopen("All.in","r",stdin);
    //freopen("a.out","w",stdout);
    srand(int(time(NULL)));
    using namespace zhangenming;
    init();
    solve();
    return 0;
}

码了半个小时QAQ