BZOJ 1588: [HNOI2002]营业额统计

1588: [HNOI2002]营业额统计

Time Limit: 5 Sec  Memory Limit: 162 MB
Submit: 18547  Solved: 7748
[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

Source

#include<algorithm>
#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cmath>
using namespace std;
 
const int maxn=32767;
typedef long long ll;
 
struct node{
    int d,n,c,f,son[2];//d?a?μ,f?a???×μ?±ào?,c?a????μ??úμ?êy,n?aí??μμ??úμ???êy 
}t[41000];
int len,root;
 
inline void update(int x)//?üD?x?ù????μ??úμ?êy 
{
    int lc=t[x].son[0],rc=t[x].son[1];
    t[x].c=t[lc].c+t[rc].c+t[x].n;
}
 
inline void add(int d,int f)//ìí?ó?μ?adμ?μ?,è?f?a???×,í?ê±,fò2è??ü?ao￠×ó 
{
    len++;
    t[len].d=d;
    t[len].n=1;
    t[len].c=1;
    t[len].f=f;
    if(d<t[f].d)
        t[f].son[0]=len;
    else
        t[f].son[1]=len;
    t[len].son[0]=t[len].son[1]=0;
}
 
inline void rotate(int x,int w)//×óDy(x,0)?ò??óòDy(x,1) 
{
    int f=t[x].f,ff=t[f].f;//x?úDy×a???°,òaè·?¨xμ????×foíòˉòˉff
    //??à′?¨á￠1??μ 
    int r,R;//r′ú±í?ù±2,R±íê???±2 
    //óD4????é?:?òx,?òμ??ù×ó,?òμ????×,?òμ?òˉòˉ 
    r=t[x].son[w];R=f;//xμ??ù×ó->×?±?μ±D??ù×ó 
    t[R].son[1-w]=r;
    if(r!=0)
        t[r].f=R;
     
    r=x;R=ff;//x->×?±?μ±D??ù×ó 
    if(t[R].son[0]==f)
        t[R].son[0]=r;
    else   
        t[R].son[1]=r;
    t[r].f=R;
     
    r=f;R=x;//xμ????×->×?±?μ±D??ù×ó 
    t[R].son[w]=r;
    t[r].f=R;
     
    update(f);//?è?üD?′|óú??2?μ?μ?f 
    update(x);//?ù?üD?é?2?μ?x 
}
 
inline void splay(int x,int rt)//??oˉêy1|?üê??aá?è?x±?3értμ?o￠×ó(×óóò???éò?) 
{
    while(t[x].f!=rt)//è?1?xμ?òˉòˉê?rt,???′x??DèòaDy×aò?′?(?àμ±óúì?ò?2?) 
    {
        int f=t[x].f,ff=t[f].f;
        if(ff==rt)
        {
            if(t[f].son[0]==x)
                rotate(x,1);
            else
                rotate(x,0);
        }
        else
        {
            if(t[ff].son[0]==f&&t[f].son[0]==x)
            {
                rotate(f,1);rotate(x,1);
            }
            else if(t[ff].son[1]==f&&t[f].son[1]==x)
            {
                rotate(f,0);rotate(x,0);
            }
            else if(t[ff].son[0]==f&&t[f].son[1]==x)
            {
                rotate(x,0);rotate(x,1);
            }
            else if(t[ff].son[1]==f&&t[f].son[0]==x)
            {
                rotate(x,1);rotate(x,0);
            }
        }
    }
    if(rt==0)
        root=x;
}
 
inline int findip(int d)//?ò?μ?adμ??úμ?μ??·,213?:è?1?2?′??úd,óD?é?üê??ó?üdμ?(?ò′ó?òD?)
{
    int x=root;
    while(t[x].d!=d)
    {
        if(d<t[x].d)
        {
            if(t[x].son[0]==0)
                break;
            else
                x=t[x].son[0];
        }
        else
        {
            if(t[x].son[1]==0)
                break;
            else
                x=t[x].son[1];
        }
    }
    return x;
}
 
inline void insert(int d)//2?è?êy?μ?adμ?ò????úμ? 
{
    if(root==0)
    {
        add(d,0);
        root=len;
        return ;
    }
    int x=findip(d);
    if(t[x].d==d)
    {
        t[x].n++;
        update(x);
        splay(x,0);
    }
    else
    {
        add(d,x);
        update(x);
        splay(len,0);
    }
}
 
inline void del(int d)//é?3yêy?μ?adμ?ò????úμ? 
{
    int x=findip(d);
    splay(x,0);//?òè?,2￠?òè?yμ±ê÷?ù 
    if(t[x].n>1)//?à??éí·Y,?í2?ó?é?μ? 
    {
        t[x].n--;
        update(x);
        return ;
    }
    if(t[x].son[0]==0&&t[x].son[1]==0)
    {
        root=0;
        len=0;
        return ;
    }
    else if(t[x].son[0]==0&&t[x].son[1]!=0)
    {
        root=t[x].son[1];
        t[root].f=0;
    }
    else if(t[x].son[0]!=0&&t[x].son[1]==0)
    {
        root=t[x].son[0];
        t[root].f=0;
    }
    else
    {
        int p=t[x].son[0];
        while(t[p].son[1]!=0)
            p=t[p].son[1];
        splay(p,x);
         
        int r=t[x].son[1],R=p;
         
        t[R].son[1]=r;
        t[r].f=R;
         
        root=R;
        t[root].f=0;
        update(R);
    }
}
 
inline int findrank(int d)//?ò???? 
{
    int x=findip(d);
    splay(x,0);
    return t[t[x].son[0]].c+1;
}
 
inline int findshuzhi(int k)//?ò?????akμ??μ 
{
    int x=root;
    while(true)
    {
        int lc=t[x].son[0],rc=t[x].son[1];
        if(k<=t[lc].c) 
            x=lc;//è￥×óo￠×ó2é?ò 
        else if(k>t[lc].c+t[x].n)
        {
            k-=t[lc].c+t[x].n;
            x=rc;
        }//è￥óòo￠×ó2é?ò 
        else
            break;//?íê??? 
    }
    splay(x,0);
    return x;
}
 
inline int findqianqu(int d)//?ò?°?y 
{
    int x=findip(d);
    splay(x,0);
    if(d<=t[x].d&&t[x].son[0]!=0) 
    {
        x=t[x].son[0];
        while(t[x].son[1]!=0)
            x=t[x].son[1];
    }
    if(t[x].d>=d)//è?1?ê?if(t[x].d>d)?ò?òμ?μ?ê?:D?óúμèóúdμ??°?y 
        x=0;
    return x;
}
 
inline int findhouji(int d)//?òoó?ì 
{
    int x=findip(d);
    splay(x,0);
    if(t[x].d<=d&&t[x].son[1]!=0)
    {
        x=t[x].son[1];
        while(t[x].son[0]!=0)
            x=t[x].son[0];
    }
    if(t[x].d<=d)
        x=0;
    return x;
}
 
int main()
{
    int n;
    ll ans=0;
    scanf("%d",&n);
    len=0;
    root=0;
    for(int i=1;i<=n;i++)
    {
        int x;
        scanf("%d",&x);
        insert(x);
        if(i==1)
        {
            ans+=x;
            continue;
        }
        int qq=findqianqu(x),hj=findhouji(x),id=findip(x);
        if(t[id].n>1)continue;
        if(qq==0)
        {
            ans+=abs(x-t[hj].d);
        }
        else if(hj==0)
        {
            ans+=abs(x-t[qq].d);
        }
        else
        {
            ans+=min(abs(x-t[qq].d),abs(x-t[hj].d));
        }
    }
    printf("%lld
",ans);
    return 0;
}