bzoj1014: [JSOI2008]火星人prefix（splay+hash+二分）

　　题目大意：一个字符串三个操作：①求两个后缀的LCP②插入一个字符③修改一个字符。

　　前几天刚学了hash+二分求lcp，就看到这题。

　　原来splay还能这么用？！原来splay模板这么好写？我以前写的splay是假的吧woc

　　splay每个节点代表一个字符，并维护这个子树代表一个子串的哈希值。因为splay旋转不破坏树结构，所以不论怎么旋转这棵splay树都能代表这个字符串。

　　预处理处理少了调了半天呜呜呜

　　赶紧跑去更新自己的splay模板

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#define which(x) (son[fa[x]][1]==x)
#define ull unsigned long long
using namespace std;
const int maxn=;
int n,m,x,y,tot,root;
int fa[maxn],son[maxn][],size[maxn],s[maxn];
ull hs[maxn],mul[maxn];
char st[maxn],ch[maxn],ch2[maxn];
void read(int &k)
{
    int f=;k=;char c=getchar();
    while(c<''||c>'')c=='-'&&(f=-),c=getchar();
    while(c<=''&&c>='')k=k*+c-'',c=getchar();
    k*=f;
}
void update(int x)
{
    size[x]=size[son[x][]]+size[son[x][]]+;
    hs[x]=hs[son[x][]]+mul[size[son[x][]]]*s[x]+mul[size[son[x][]]+]*hs[son[x][]];
}
void rotate(int x)
{
    int f=fa[x];
    bool k=which(x);
    son[f][k]=son[x][!k];
    son[x][!k]=f;
    son[fa[f]][which(f)]=x;
    fa[son[f][k]]=f;
    fa[x]=fa[f];
    fa[f]=x;
    size[x]=size[f];
    hs[x]=hs[f];
    update(f);
}
void splay(int x,int g)
{
    while(fa[x]!=g)
    {
        int f=fa[x];
        if(fa[f]==g)
        {
            rotate(x);
            break;
        }
        if(which(x)^which(f))rotate(x);
        else rotate(f);
        rotate(x);
    }
    if(!g)root=x;
}
int rank(int x,int k)
{
    if(k<=size[son[x][]])return rank(son[x][],k);
    else if(k==(size[son[x][]]+))return x;
    else return rank(son[x][],k-size[son[x][]]-);
}
int build(int l,int r,int f)
{
    if(l>r)return ;
    int x=++tot,mid=(l+r)>>;
    fa[x]=f;s[x]=(int)(st[mid]-'a')+;
    son[x][]=build(l,mid-,x);
    son[x][]=build(mid+,r,x);
    update(x);
    return x;
}
void insert(int k,int ch)
{
    int x=rank(root,k),y=rank(root,k+);
    splay(x,);splay(y,x);
    s[++tot]=ch;
    fa[tot]=y;son[y][]=tot;
    update(tot);update(y);update(x);
}
void change(int k,int ch)
{
    int x=rank(root,k);
    splay(x,);
    s[x]=ch;
    update(x);
}
int lcp(int kx,int ky)
{
    int l=,r=n;
    while(l<r)
    {
        int mid=(l+r+)>>;
        if(ky+mid>n+)
        {
            r=mid-;
            continue;
        }
        int x=rank(root,kx-),y=rank(root,kx+mid);
        splay(x,);splay(y,x);
        ull haxi=hs[son[y][]];
        x=rank(root,ky-),y=rank(root,ky+mid);
        splay(x,);splay(y,x);
        if(haxi==hs[son[y][]])l=mid;
        else r=mid-;
    }
    return l;
}
int main()
{
    scanf("%s",st+);
    n=strlen(st+);
    mul[]=;
    for(int i=;i<=;i++)mul[i]=mul[i-]*;
    root=build(,n+,);
    read(m);
    for(int i=;i<=m;i++)
    {
        scanf("%s",ch);
        if(ch[]=='Q')
        {
            read(x);read(y);
            if(x>y)swap(x,y);
            if(x!=y)printf("%d\n",lcp(x+,y+));
            else printf("%d\n",n-x+);
        }
        if(ch[]=='R')
        {
            read(x);
            scanf("%s",ch2);
            change(x+,(int)(ch2[]-'a')+);
        }
        if(ch[]=='I')
        {
            read(x);
            scanf("%s",ch2);
            insert(x+,(int)(ch2[]-'a')+);
            n++;
        }
    }
    return ;
}