bzoj 2827 千山鸟飞绝 平衡树

原创 2017年01月04日 08:10:06

splay和我的常数加到一起简直是灾难。。。
坐标没有用,把它离散掉。
然后splay维护一个最大值和次大值和最大size的标记。
插入的时候查询一下,打个标记。

#include <bits/stdc++.h>
using namespace std;
#define N 410000
#define M 810000
#define ll long long
#define ls(x) ch[x][0]
#define rs(x) ch[x][1]
#define which(x) (ch[fa[x]][1]==x)
#define icmp(x,y) (x^y ? v[x]^v[y] ? v[x]>v[y] ? 1:-1 : x>y ? 1:-1 :0)
const int inf=~0U>>1;
int n,m,cnt,tot;
int v[M],a[N],av[M],as[M],rem[N],root[N];
int size[M],L[N],R[N],l,r;
ll val[N],pos[N],org[N];
int ch[M][2],fa[M],mx[M],sx[M],sum[M];
char getc()
{
    static const int LEN = 4096;
    static char buf[LEN],*S=buf,*T=buf;
    if(S == T)
    {
        T = (S=buf)+fread(buf,1,LEN,stdin);
        if(S == T)return EOF;
    }
    return *S++;
}
int read()
{
    static char ch;
    static int D;
    while(!isdigit(ch=getc()));
    for(D=ch-'0'; isdigit(ch=getc());)
        D=(D<<3)+(D<<1)+(ch-'0');
    return D;
}
void upd(int x)
{
    as[x]=max(as[x],sum[x]);
    int t=mx[x]==x ? sx[x]:mx[x];
    if(v[t]>v[av[x]])av[x]=t;
}
void cal(int x,int y)
{
    if(!x||!y||mx[x]==y||sx[x]==y)return;
    if(v[y]>v[mx[x]])sx[x]=mx[x],mx[x]=y;
    else if(v[y]>v[sx[x]])sx[x]=y;
}
void pushdown(int x)
{
    if(sum[x])
    {
        sum[ls(x)]=max(sum[ls(x)],sum[x]);
        sum[rs(x)]=max(sum[rs(x)],sum[x]);
    }
    cal(ls(x),mx[x]);cal(rs(x),mx[x]);
    cal(ls(x),sx[x]);cal(rs(x),sx[x]);
    upd(ls(x));upd(rs(x));
    mx[x]=sx[x]=0;sum[x]=0;
}
void pushup(int x)
{size[x]=size[ls(x)]+size[rs(x)]+1;}
void down(int x)
{
    if(fa[x])down(fa[x]);
    pushdown(x);
}
void rotate(int x)
{
    int y=fa[x],k=which(x);
    ch[y][k]=ch[x][k^1];
    ch[x][k^1]=y;
    ch[fa[y]][which(y)]=x;
    fa[x]=fa[y];fa[y]=x;
    fa[ch[y][k]]=y;

    pushup(y);pushup(x);
}
void splay(int x,int tar)
{
    down(x);
    while(fa[x]!=tar)
    {
        int y=fa[x];
        if(fa[y]==tar)rotate(x);
        else
        {
            if(which(x)^which(y))rotate(x);
            else rotate(y);
            rotate(x);
        }
    }
}
void pre(int x,int y)
{
    if(!x)return;
    if(icmp(x,y)<0)
    {
        if(icmp(x,l)>0)l=x;
        pre(ch[x][1],y);
    }
    else pre(ch[x][0],y);
}
void nex(int x,int y)
{
    if(!x)return;
    if(icmp(x,y)>0)
    {
        if(icmp(x,r)<0)r=x;
        nex(ch[x][0],y);
    }
    else nex(ch[x][1],y);
}
int find(int x,int K)
{
    int t=size[ch[x][0]];
    if(t+1==K)return x;
    if(t>=K)return find(ch[x][0],K);
    return find(ch[x][1],K-t-1);
}
void build(int x,int y,int pos)
{
    root[pos]=x;
    L[pos]=x;R[pos]=y;
    ch[x][1]=y;v[x]=-1;v[y]=inf;
    fa[y]=x;size[x]=2;size[y]=1;
}
void del(int pos,int x)
{
    int &rt=root[pos];
    l=L[pos];r=R[pos];
    pre(rt,x);nex(rt,x);
    splay(l,0);rt=l;splay(r,l);
    size[l]--;size[r]--;
    ch[r][0]=fa[x]=0;
}
void ins(int pos,int x)
{
    int &rt=root[pos];
    if(size[rt]>2)
    {
        int t=find(rt,size[rt]-1);
        if(icmp(t,av[x])>0)av[x]=t;
    }
    l=L[pos];r=R[pos];
    pre(rt,x);nex(rt,x);
    splay(l,0);rt=l;splay(r,l);
    size[l]++;size[r]++;
    ch[r][0]=x;fa[x]=r;size[x]=1;
    cal(rt,x);sum[rt]=max(sum[rt],size[rt]-3);
    upd(rt);
}
int main()
{
    n=read();
    for(int i=1,x,y;i<=n;i++)
    {
        v[i]=read();x=read();y=read();
        org[i]=val[++cnt]=(ll)x*inf+y;
    }
    m=read();
    for(int i=1,x,y;i<=m;i++)
    {
        a[i]=read();x=read();y=read();
        pos[i]=val[++cnt]=(ll)x*inf+y;
    }
    sort(val+1,val+1+cnt);
    cnt=unique(val+1,val+1+cnt)-val-1;
    tot=n;
    for(int i=1;i<=cnt;i++)
        build(tot+1,tot+2,i),tot+=2;
    for(int i=1;i<=n;i++)
    {
        org[i]=lower_bound(val+1,val+1+cnt,org[i])-val;
        ins(org[i],i);
    }
    for(int i=1,t;i<=m;i++)
    {
        t=a[i];del(org[t],t);
        pos[i]=lower_bound(val+1,val+1+cnt,pos[i])-val;
        org[t]=pos[i];ins(org[t],t);
    }
    for(int i=1;i<=n;i++)
    {
        down(i);
        printf("%lld\n",(ll)v[av[i]]*as[i]);
    }
    return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

【BZOJ2827】千山鸟飞绝 离散化+splay

题解: 首先先把坐标离散化一下, 然后对于每个坐标点我们建一棵平衡树,每次插入操作后给整颗平衡树下传一下需求的两个标记。 注意: splay有的人(比如我)习惯每棵都先建-inf、inf两...
  • Vmurder
  • Vmurder
  • 2015年05月14日 15:17
  • 1523

[BZOJ1503][NOI2004]郁闷的出纳员(平衡树splay)

海水是咸的,因为它包含了人世所有的泪水;泪水是咸的,因为它承载了太多的欢乐和苦悲。...

【bzoj3262】【陌上花开】【树状数组套平衡树】

Description 有n朵花,每朵花有三个属性:花形(s)、颜色(c)、气味(m),又三个整数表示。现要对每朵花评级,一朵花的级别是它拥有的美丽能超过的花的数量。定义一朵花A比另一朵花B要美...

BZOJ3224:普通平衡树(含SBT、Treap、Splay模板)

传送门 这道题就是裸的平衡树的题目,操作都是模板,也就不多说了,重在代码。SBT:#include struct SBT { int v, sz, ch[2], cnt; } t[30000...

bzoj 2733【hnoi2012】永无乡(平衡树)

时间限制:1秒 内存限制:64M 【问题描述】   永无乡包含 n 座岛,编号从 1 到 n,每座岛都有自己的独一无二的重要度,按照重要度可 以将这 n 座岛排名,名次用 1 到 n 来表示。某些...

BZOJ 1146 CTSC 2008 网络管理 Network 树链剖分+二分答案+平衡树

题目大意:有n个路由器,他们由n-1条边连接(形成一棵树)。每一个路由器有一个延时值。有两种操作: 1.查询树上x,y两点之间的路径上第k大的权值是多少 2.修改x位置的权值为y 思路:当我大概想...

BZOJ 2733 HNOI 2012 永无乡 平衡树启发式合并

题目大意:有一些岛屿,一开始由一些无向边连接。后来也有不断的无向边加入,每一个岛屿有个一独一无二的重要度,问任意时刻的与一个岛屿联通的所有岛中重要度第k大的岛的编号是什么。 思路:首先连通性...

bzoj 3223: Tyvj 1729 文艺平衡树(splay 模板题 区间翻转)

3223: Tyvj 1729 文艺平衡树 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 3321  Solved: 1887 [Submit...

【BZOJ 3224】普通平衡树

BZOJ 3224题意有n个操作,一共有6种。 编号 操作 含义 1 insert(x) 插入一个数x 2 remove(x) 删去一个数x,多个删去一个 3 ...
  • Ripped
  • Ripped
  • 2017年02月09日 14:10
  • 151

BZOJ3224普通平衡树splay,SBT代码

splay#include #include #include using namespace std; int v[100010],fa[100010],ch[100010][2],sz[10001...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:bzoj 2827 千山鸟飞绝 平衡树
举报原因:
原因补充:

(最多只允许输入30个字)