【权值分块】bzoj1588 [HNOI2002]营业额统计

权值分块就是快……Rank5……

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define N 38001
#define INF 2147483647
#define min(a,b) (((a)<(b))?(a):(b))
struct Point{int v,p;}t[N];
bool operator < (const Point &a,const Point &b){return a.v<b.v;}
int n,ma[N],a[N],en,sum,sz,l[201],r[201],sumv[201],num[N],b[N],ans;
inline void R(int &x){
    char c=0;int f=1;
    for(;c<'0'||c>'9';c=getchar())if(c=='-')f=-1;
    for(x=0;c>='0'&&c<='9';c=getchar())(x*=10)+=(c-'0');
    x*=f;
}
void makeblock()
{
    sz=sqrt(en); if(!sz) sz=1;
    for(sum=1;sum*sz<en;sum++)
      {
          l[sum]=r[sum-1]+1;
          r[sum]=sz*sum;
          for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
      }
    l[sum]=r[sum-1]+1;
    r[sum]=en;
    for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
}
inline void Insert(const int &x){b[x]++; sumv[num[x]]++;}
inline int Next(const int &x)
{
    if(b[x]>1) return x;
    for(int i=x+1;i<=r[num[x]];i++) if(b[i]) return i;
    for(int i=num[x]+1;i<=sum;i++) if(sumv[i])
      for(int j=l[i];;j++)
        if(b[j]) return j;
    return INF;
}
inline int Pre(const int &x)
{
    if(b[x]>1) return x;
    for(int i=x-1;i>=l[num[x]];i--) if(b[i]) return i;
    for(int i=num[x]-1;i>=1;i--) if(sumv[i])
      for(int j=r[i];;j--)
        if(b[j]) return j;
    return INF;
}
int main()
{
    R(n); for(int i=1;i<n;i++)
      {
          R(t[i].v);
          t[i].p=i;
      }
    if(scanf("%d",&t[n].v)==EOF) t[n].v=0; t[n].p=n;
    sort(t+1,t+n+1);
    ma[a[t[1].p]=++en]=t[1].v;
    for(int i=2;i<=n;i++)
      {
          if(t[i].v!=t[i-1].v) en++;
          ma[a[t[i].p]=en]=t[i].v;
      } makeblock();
    Insert(a[1]); ans+=ma[a[1]];
    for(int i=2;i<=n;i++)
      {
          Insert(a[i]);
          int Up=Next(a[i]),Down=Pre(a[i]);
          if(Up==INF) ans+=(ma[a[i]]-ma[Down]);
          else if(Down==INF) ans+=(ma[Up]-ma[a[i]]);
        else ans+=min(ma[Up]-ma[a[i]],ma[a[i]]-ma[Down]);
      } printf("%d
",ans);
    return 0;
}