BZOJ 1588:[HNOI2002]营业额统计

裸的平衡树
辣鸡BZOJ
昨天我的splay莫名其妙一个 void v o i d 打成 int i n t 就特么狂T
(黑人问号)
今天又抽风一个下午&&一个晚上
边上那神犇去问管理员,丫说在修复等俩天
mdzz m d z z
还好昨儿向他要了数据。。
写了个替罪羊树,自己手测了一发,大约是能过的。。。
splay改的不堪入目了而且风格和某神犇极(yi)其(mao)相(yi)似(yang)所以不贴了。。
感觉替罪羊树很好写呀,为什么没人学呢。。
明明比treap好写呀,而且比treap快多了
替罪羊树就是当某个节点不够平衡时,暴力重构整棵子树
具体如何不够平衡嘿嘿嘿看你是否够强迫症了
我看到的版本是一边子树大小超过 0.7 0.7 或者 0.75 0.75 算不平衡,如果强迫症末期的患者也许可以试试 0.65 0.65 ,但是我想(只是我想)重构次数太多反而会拖慢程序。。
关于重构,先求出 dfs d f s 序,然后和二分有点类似,具体看代码吧,其实很好懂。。
然后找前驱和后继和其他平衡树一样的
插入的话也差不多,就是要判断是不是不平衡。。
删除操作下次再写。。看这篇BZOJ 3224:Tyvj 1728 普通平衡树


#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define g getchar()
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
inline ll read(){
    ll x=0,f=1;char ch=g;
    for(;ch<'0'||ch>'9';ch=g)if(ch=='-')f=-1;
    for(;ch>='0'&&ch<='9';ch=g)x=x*10+ch-'0';
    return x*f;
}
inline void out(ll x){
    int a[25],t=0;
    if(x<0)putchar('-'),x=-x;
    for(;x;x/=10)a[++t]=x%10;
    for(int i=t;i;--i)putchar('0'+a[i]);
    if(t==0)putchar('0');
    putchar('\n');
}
int sz,size[50005],t[50005][2],fa[50005],v[50005],id,q[50005],top,rt,ans1,ans2,n,ans;
int ins(int &k,int x,int f){
    if(!k){k=++sz;size[k]=1;fa[k]=f;return 0;}
    if(x<v[k])ins(t[k][0],x,k);else ins(t[k][1],x,k);
    if(max(size[t[k][0]],size[t[k][1]])>size[k]*0.7)id=k;//如果不平衡,记录要重构的子树的根
    size[k]=size[t[k][0]]+size[t[k][1]]+1;//size记录节点的重量
}
int build(int &k,int l,int r){
    if(l>r)return 0;
    int mid=(l+r)>>1;
    k=q[mid];
    build(t[k][0],l,mid-1);
    build(t[k][1],mid+1,r);
    fa[t[k][0]]=fa[t[k][1]]=k;
    size[k]=size[t[k][0]]+size[t[k][1]]+1;
}//构建
void dfs(int x){
    if(!x)return;
    dfs(t[x][0]);
    q[++top]=x;
    dfs(t[x][1]);
}//求dfs序
int re_build(int x){
    int pre=fa[x];
    dfs(x);
    for(int i=1;i<=top;++i)
    t[q[i]][0]=t[q[i]][1]=fa[q[i]]=size[q[i]]=0;
    if(!pre){build(rt,1,top);}//如果x是总树的根,重构总树
    else
    if(t[pre][0]==x)build(t[pre][0],1,top);
    else build(t[pre][1],1,top);
    fa[t[pre][0]]=fa[t[pre][1]]=pre;
    top=id=0;
}//重构以x为根的树
void lower(int k,int x){
    if(!k)return;if(v[k]<=x){ans1=v[k];lower(t[k][1],x);}else lower(t[k][0],x);
}
void higher(int k,int x){
    if(!k)return;if(v[k]>=x){ans2=v[k];higher(t[k][0],x);}else higher(t[k][1],x);
}

int main(){
//  freopen("turnover0.in","r",stdin);
//  freopen("","w",stdout);
    n=read();
    for(int i=1;i<=n;++i){
        int x;
        if(scanf("%d",&x)==EOF)x=0;
        ans1=-inf;ans2=inf;
        lower(rt,x);
        higher(rt,x);
        ans+=i==1?x:min(x-ans1,ans2-x);
        v[i]=x;
        ins(rt,x,0);
        if(id)re_build(id);
    }
    out(ans);
    return 0;
}
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值