P4036 [JSOI2008]火星人(splay+hash+二分)

P4036 [JSOI2008]火星人

Splay维护hash,查询二分

$a[x].vl=a[lc].vl*ha[a[rc].sz+1]+a[x].w*ha[a[rc].sz]+a[rc].vl$

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef unsigned long long ull;
int read(){
    char c=getchar(); int x=0;
    while(c<'0'||c>'9') c=getchar();
    while('0'<=c&&c<='9') x=x*10+c-48,c=getchar();
    return x;
}
#define N 300005
const int bas=13331;
int n,Q; char q[N],q1[3]; ull ha[N];
struct node{int ch[2],sz,fa; ull vl,w;}a[N];
struct Splay{
    #define lc a[x].ch[0]
    #define rc a[x].ch[1]
    #define mid (l+r)/2
    int rt,u;
    inline void up(int x){
        a[x].sz=a[lc].sz+a[rc].sz+1,
        a[x].vl=a[lc].vl*ha[a[rc].sz+1]+a[x].w*ha[a[rc].sz]+a[rc].vl; 
    }
    inline int whi(int y,int x){return a[y].ch[1]==x;}
    void turn(int x,int &k){
        int y=a[x].fa,z=a[y].fa,l=whi(y,x),r=l^1;
        if(y==k) k=x;
        else a[z].ch[whi(z,y)]=x;
        a[a[x].ch[r]].fa=y; a[y].fa=x; a[x].fa=z;
        a[y].ch[l]=a[x].ch[r]; a[x].ch[r]=y;
        up(y); up(x);
    }
    void splay(int x,int &k){
        for(;x!=k;turn(x,k)){
            int y=a[x].fa,z=a[y].fa;
            if(y!=k) turn(whi(z,y)==whi(y,x)?y:x,k);
        }
    }
    int kth(int k){
        int x=rt;
        while(a[lc].sz+1!=k){
            if(a[lc].sz>=k) x=lc;
            else k=k-a[lc].sz-1,x=rc;
        }return x;
    }
    inline void rep(int x,ull k){a[x=kth(x)].w=k,splay(x,rt);}
    void ins(int x,ull k){
        x=kth(x); ++u;
        if(!rc) rc=u;
        else{for(x=rc;lc;x=lc); lc=u;}
        a[u].fa=x; a[u].w=k; up(u); splay(u,rt);
    }
    ull Fi(int x,int y){
        int L=kth(x-1),R=kth(y+1);
        splay(L,rt); splay(R,a[L].ch[1]);
        return a[a[R].ch[0]].vl;
    }
    int ask(int x,int y){
        if(x>y) swap(x,y);
        int l=1,r=u-y,re=0;
        while(l<=r){
            if(Fi(x,x+mid-1)!=Fi(y,y+mid-1)) r=mid-1;
            else re=mid,l=mid+1;
        }return re;
    }
    int build(int Fa,int l,int r){
        a[mid].w=(ull)q[mid-1]; a[mid].fa=Fa; ++u;
        if(l<mid) a[mid].ch[0]=build(mid,l,mid-1);
        if(r>mid) a[mid].ch[1]=build(mid,mid+1,r);
        up(mid); return mid;
    }
}S;
int main(){
    scanf("%s",q+1); n=strlen(q+1); ha[0]=1;
    for(int i=1;i<N;++i) ha[i]=ha[i-1]*bas;
    S.rt=S.build(0,1,n+2); Q=read();//建树时记得加入边界
    for(int i=1,w1,w2;i<=Q;++i){
        scanf("%s",q1);
        if(q1[0]=='Q') w1=read(),w2=read(),printf("%d\n",S.ask(w1+1,w2+1));
        if(q1[0]=='R') w1=read(),scanf("%s",q1),S.rep(w1+1,(ull)q1[0]);
        if(q1[0]=='I') w1=read(),scanf("%s",q1),S.ins(w1+1,(ull)q1[0]);
    }return 0;
}

 

转载于:https://www.cnblogs.com/kafuuchino/p/11443071.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值