bzoj2258 pku2758 Checking the Text 文本校对

http://www.elijahqi.win/archives/2840
Description

为了给Wind买生日礼物，Jiajia不得不找了一份检查文本的工作。这份工作很无聊：给你一段文本
要求比对从文本中某两个位置开始能匹配的最大长度是多少。但比无聊更糟糕的是，Jiajia的经理
还可能往文本里面插入一些字符。
Jiajia想用一个程序来解决这些繁琐的工作。这个程序的速度要足够快，因为Wind的生日就快要到了
Jiajia必须赚到足够多的钱，也就是处理足够多的文本。

Input

输入文件第一行是原始文本。
输入文件第二行是操作数n。此后n行，每行描述一条命令，命令有两种形式：
I ch p：表示将一个字符ch插入到当前文本的第p个字符之前，如果p大于当前文本长度则表示插入到当前文本末尾；
Q i j：表示询问当前文本从原始文本的第i个和第j个字符现在所在的位置开始能匹配的字符是多少。
你可以认为文本初始长度不超过50000，I命令最多200条，Q命令最多20000条。

Output

对于每条Q命令输出一行，为最长匹配长度。

Sample Input

abaab
5
Q 1 2
Q 1 3
I a 2
Q 1 2
Q 1 3

Sample Output

0
1
0
3

HINT

Source

看到这题一眼火星人prefix 可惜蒟蒻我人傻自带大常数 Poj过不去..

顺便复习下splay板子 火星人prefix见：http://www.elijahqi.win/archives/2344

#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define N 55000
#define g 27
#define ll unsigned int
using namespace std; 
inline int read(){
    int x=0;char ch=getchar();
    while(!isdigit(ch)) ch=getchar();
    while(isdigit(ch)) x=x*10+ch-'0',ch=getchar();
    return x;
}
char s[N];int fa[N],size[N],c[N][2],v[N],cnt,root;
ll hs[N],p[N];
inline void update(int x){
    int l=c[x][0],r=c[x][1];
    size[x]=size[l]+size[r]+1;
    hs[x]=hs[l]*p[size[r]+1]+v[x]*p[size[r]]+hs[r];
}
inline void rotate(int x,int &tar){
    int y=fa[x],z=fa[y],l=c[y][1]==x,r=l^1;
    if (y==tar) tar=x;else c[z][c[z][1]==y]=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 &tar){
    while(x!=tar){
        int y=fa[x],z=fa[y];
        if (y!=tar) {
            if (c[y][0]==x^c[z][0]==y) rotate(x,tar);else rotate(y,tar);
        }rotate(x,tar);
    }
}
inline void build(int f,int l,int r){
    if (l>r) return;int mid=l+r>>1;
    v[mid]=hs[mid]=s[mid]-'a'+1;size[mid]=1;
    c[f][mid>f]=mid;fa[mid]=f;if(l==r) return;
    build(mid,l,mid-1);build(mid,mid+1,r);
    update(mid);
}
inline int find(int x,int p){
    int l=c[x][0],r=c[x][1];
    if (size[l]+1==p) return x;
    if(p<=size[l]) return find(l,p);else return find(r,p-size[l]-1);
}
inline int split(int x,int y){
    int xx=find(root,x),yy=find(root,y);
    splay(yy,root);splay(xx,c[root][0]);return c[xx][1];
}
inline bool check(int l,int r,int mid){
    int x=split(l,l+mid+1);ll tmp1=hs[x];
    x=split(r,r+mid+1);ll tmp2=hs[x];return tmp1==tmp2;
}
inline void print(int x){
    if (c[x][0]) print(c[x][0]);
    printf("%c %llu\n",v[x]+'a'-1,hs[x]);
    if (c[x][1]) print(c[x][1]);
}
int main(){
    freopen("bzoj2258.in","r",stdin);
    scanf("%s",s+2);int n=strlen(s+2);p[0]=1;
    for (int i=1;i<N;++i) p[i]=p[i-1]*g;cnt=n+2;
    build(0,1,n+2);root=n+3>>1;int m=read();//print(root);
    while(m--){static char op[2];
        scanf("%s",op);
        if (op[0]=='Q'){
            int x=read(),y=read();splay(x+1,root);x=size[c[x+1][0]];
            splay(y+1,root);y=size[c[y+1][0]];
            int l=1,r=min(cnt-x-1,cnt-y-1);
            while(l<=r){
                int mid=l+r>>1;
                if (check(x,y,mid)) l=mid+1;else r=mid-1;
            }printf("%d\n",r);continue;
        }
        if (op[0]=='I'){
            static char ch[2];scanf("%s",ch);int p=read();
            if (p-2>size[root]) p=size[root]-1;int xx=find(root,p),yy=find(root,p+1);
            splay(yy,root);splay(xx,c[root][0]);v[++cnt]=ch[0]-'a'+1;hs[cnt]=ch[0]-'a'+1;fa[cnt]=xx;
            c[xx][1]=cnt;size[cnt]=1;update(xx);update(root);
            //print(root);
        }
    }
    return 0;
}