bzoj1014 [JSOI2008]火星人prefix（二分答案+哈希+平衡树）

最新推荐文章于 2019-08-15 21:12:00 发布

Icefox_zhx

最新推荐文章于 2019-08-15 21:12:00 发布

阅读量436

点赞数

CC 4.0 BY-SA版权

分类专栏： bzoj 平衡树 --------字符串-------- Hash 二分答案

本文链接：https://blog.youkuaiyun.com/icefox_zhx/article/details/79209246

bzoj 同时被 3 个专栏收录

764 篇文章

订阅专栏

--------字符串--------

114 篇文章

订阅专栏

二分答案

66 篇文章

订阅专栏

本文介绍了一种使用Splay树和Hash算法来高效解决最长公共前缀(LCP)问题的方法。通过构建Splay树并利用Hash进行节点值的更新与查询，实现了快速的插入、修改和查询操作。代码示例展示了具体的实现细节。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

我每次询问，用二分答案+hash去求lcp。
splay维护一下插入和修改。点p的hash值就是
$hs[l]*(sz[r]+1)+hs[p]*sz[r]+hs[r]$

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define inf 0x3f3f3f3f
#define N 100010
#define ll long long
#define uint unsigned int 
#define base 29
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,fa[N],sz[N],c[N][2],rt,v[N],cnt;
char s[N];uint bin[N],hs[N];
inline void update(int p){
    int l=c[p][0],r=c[p][1];
    sz[p]=sz[l]+sz[r]+1;
    hs[p]=hs[l]*bin[sz[r]+1]+v[p]*bin[sz[r]]+hs[r];
}
inline void build(int &p,int l,int r){
    int mid=l+r>>1;p=mid;v[p]=s[mid-1]-'a';
    if(l<mid) build(c[p][0],l,mid-1),fa[c[p][0]]=p;
    if(r>mid) build(c[p][1],mid+1,r),fa[c[p][1]]=p;
    update(p);
}
inline int find(int p,int k){
    if(sz[c[p][0]]+1==k) return p;
    if(k<=sz[c[p][0]]) return find(c[p][0],k);
    return find(c[p][1],k-sz[c[p][0]]-1);
}
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(fa[y]);
}
inline void splay(int x,int &k){
    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);
            rotate(y,k);
        }rotate(x,k);
    }
}
inline int split(int xx,int yy){
    int x=find(rt,xx),y=find(rt,yy);
    splay(x,rt);splay(y,c[x][1]);return c[y][0];
}
inline void change(int pos,int val){
    int x=split(pos-1,pos+1),y=fa[x];
    v[x]=val;update(x);update(y);update(fa[y]);
}
inline void ins(int pos,int val){
    int xx=pos,yy=pos+1;
    int x=find(rt,xx),y=find(rt,yy);
    splay(x,rt);splay(y,c[x][1]);c[y][0]=++cnt;fa[cnt]=y;
    v[cnt]=val;update(cnt);update(y);update(fa[y]);
}
inline uint hash(int l,int r){
    int x=split(l-1,r+1);return hs[x];
}
inline bool jud(int x,int y,int len){
    return hash(x,x+len-1)==hash(y,y+len-1);
}
inline void solve(int x,int y){
    int l=0,r=min(cnt-2-x+1,cnt-2-y+1);x++;y++;
    while(l<=r){
        int mid=l+r>>1;
        if(jud(x,y,mid)) l=mid+1;
        else r=mid-1;
    }printf("%d\n",l-1);
}
int main(){
//  freopen("a.in","r",stdin);
    scanf("%s",s+1);n=strlen(s+1);bin[0]=1;
    for(int i=1;i<=N-1;++i) bin[i]=bin[i-1]*base;
    build(rt,1,n+2);m=read();cnt=n+2;
    while(m--){
        char op[5];scanf("%s",op+1);int x=read();
        if(op[1]=='Q'){solve(x,read());continue;}
        scanf("%s",s+1);
        if(op[1]=='R') change(x+1,s[1]-'a');
        else ins(x+1,s[1]-'a');
    }return 0;
}