bzoj2827 千山鸟飞绝(平衡树+map)

47 篇文章 0 订阅

把坐标用map离散化掉。我们考虑把一个点x加到集合y中的影响:
首先是对x点的影响:用集合y中的最大值更新ans1,用集合大小更新ans2.
然后是对集合y的影响:每个点的ans1都可以由v[x]更新,ans2都可以由集合大小更新。
因此我们就要维护集合最大值,集合大小,还要能够打标记,支持插入删除。
用splay即可。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <utility>
#include <map>
using namespace std;
#define ll long long
#define N 30010
#define M 300010
#define pa pair<int,int>
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,v[N],sz[N],fa[N],c[N][2],mx[N],tag1[N],tag2[N],rt[M],owo=0,pos[N];
map<pa,int> mp;int ans1[N],ans2[N],q[N];
inline void pushdown(int p){
    int l=c[p][0],r=c[p][1];
    if(tag1[p]){
        ans1[p]=max(ans1[p],tag1[p]);
        tag1[l]=max(tag1[l],tag1[p]);
        tag1[r]=max(tag1[r],tag1[p]);tag1[p]=0;
    }if(tag2[p]){
        ans2[p]=max(ans2[p],tag2[p]);
        tag2[l]=max(tag2[l],tag2[p]);
        tag2[r]=max(tag2[r],tag2[p]);tag2[p]=0;
    }
}
inline void update(int p){
    int l=c[p][0],r=c[p][1];
    sz[p]=sz[l]+sz[r]+1;
    mx[p]=max(max(mx[l],mx[r]),v[p]);
}
inline void rotate(int x,int &k){
    int y=fa[x],z=fa[y],l=x==c[y][1],r=l^1;
    if(y==k) k=x;
    else c[z][y==c[z][1]]=x;
    fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
    c[y][l]=c[x][r];c[x][r]=y;update(y);update(x);
}
inline void splay(int x,int &k){
    int top=0;for(int xx=x;xx;xx=fa[xx]) q[++top]=xx;
    while(top) pushdown(q[top--]);
    while(x!=k){
        int y=fa[x],z=fa[y];
        if(y!=k){
            if(x==c[y][1]^y==c[z][1]) rotate(x,k);
            else rotate(y,k);
        }rotate(x,k);
    }
}
inline void ins(int &p,int x,int Fa){
    if(!p){p=x;sz[p]=1;fa[p]=Fa;c[p][0]=c[p][1]=0;tag1[p]=tag2[p]=0;mx[p]=v[p];return;}
    pushdown(p);
    if(!c[p][0]) ins(c[p][0],x,p);
    else ins(c[p][1],x,p);
}
inline void del(int &rt,int x){
    splay(x,rt);
    if(c[x][0]*c[x][1]==0){
        rt=c[x][0]+c[x][1];fa[rt]=0;c[x][0]=c[x][1]=0;return;
    }int pre=c[x][0],succ=c[x][1];
    while(c[pre][1]) pre=c[pre][1];
    while(c[succ][0]) succ=c[succ][0];
    splay(pre,rt);splay(succ,c[rt][1]);c[succ][0]=fa[x]=0;
    update(succ);update(pre);
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();for(int i=1;i<=n;++i){
        v[i]=read();int x=read(),y=read();
        if(!mp[make_pair(x,y)]) mp[make_pair(x,y)]=++owo;
        pos[i]=mp[make_pair(x,y)];
        ans1[i]=max(ans1[i],mx[rt[pos[i]]]);
        ans2[i]=max(ans2[i],sz[rt[pos[i]]]);
        tag1[rt[pos[i]]]=v[i];tag2[rt[pos[i]]]=sz[rt[pos[i]]];
        ins(rt[pos[i]],i,0);splay(i,rt[pos[i]]);
    }int m=read();while(m--){
        int id=read(),x=read(),y=read();
        if(!mp[make_pair(x,y)]) mp[make_pair(x,y)]=++owo;
        del(rt[pos[id]],id);pos[id]=mp[make_pair(x,y)];
        ans1[id]=max(ans1[id],mx[rt[pos[id]]]);
        ans2[id]=max(ans2[id],sz[rt[pos[id]]]);
        tag1[rt[pos[id]]]=max(tag1[rt[pos[id]]],v[id]);
        tag2[rt[pos[id]]]=max(tag2[rt[pos[id]]],sz[rt[pos[id]]]);
        ins(rt[pos[id]],id,0);splay(id,rt[pos[id]]);
    }for(int i=1;i<=n;++i) splay(i,rt[pos[i]]);
    for(int i=1;i<=n;++i) printf("%lld\n",(ll)ans1[i]*ans2[i]);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值