[bzoj1014][JSOI2008]火星人prefix

Description

我们定义LCQ(x,y)为x这个后缀和y这个后缀的LCP(最长公共前缀)
给出一个字符串，可以在中间插入某一个字符，或者修改某一个字符，并且维护LCQ.
Q<=1.5*10^5，|S|<=10^5，字符集为所有的小写字母，询问数<=10^4.

Solution

一眼看到还以为是裸SA（智商堪忧）
后来发现有动态插入，立马想到splay。
维护什么呢？Hash值！
然后每次询问二分答案判定是否合法。
打一个splay需要20min，但用splayA一题可能需要2h+
这道题常数写大了，一直TLE

然后一怒之下把LL扔掉了，然后就。。。

踩线过啦！（各位4ms的爷请收下我的膝盖Orz）
所以，当你看到你跑多了100ms的话，LL神教向你招手。

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define N 100010
using namespace std;
typedef long long ll;
const int mo=9875321;
int f[N],t[N][2],size[N],h[N],mi[N],key[N],tot,n,m,x,y,root;
char s[N],ch[2];
void updata(int x) {
    int l=t[x][0],r=t[x][1];
    size[x]=size[l]+size[r]+1;
    h[x]=h[l]+(ll)key[x]*mi[size[l]]%mo+(ll)h[r]*mi[size[l]+1]%mo;
    h[x]%=mo;
}
int build(int l,int r,int fa) {
    if (l>r) return 0;
    int m=(l+r)/2;
    f[m]=fa;size[m]=1;
    h[m]=key[m]=s[m]-'a';
    if (l==r) return m;
    t[m][0]=build(l,m-1,m);
    t[m][1]=build(m+1,r,m);
    updata(m);return m;
}
int son(int x) {
    if (t[f[x]][0]==x) return 0;else return 1;
}
void rotate(int x) {
    int y=f[x],z=son(x);f[x]=f[y];
    if (f[x]) t[f[x]][son(y)]=x;
    if (t[x][1-z]) f[t[x][1-z]]=y;
    t[y][z]=t[x][1-z];t[x][1-z]=y;f[y]=x;
    updata(y);updata(x);
}
void splay(int x,int y) {
    while (f[x]!=y) {
        if (f[f[x]]!=y) 
            if (son(x)==son(f[x])) rotate(f[x]);
            else rotate(x);
        rotate(x);
    }
    if (!y) root=x;
}
int find(int x,int id) {
    if (!x) return 0;
    int l=t[x][0],r=t[x][1];
    if (size[l]+1==id) return x;
    if (size[l]+1>id) return find(l,id);
    else return find(r,id-size[l]-1);
}
void insert(int x,int c) {
    int y=find(root,x);splay(y,0);
    int z=find(root,x+1);splay(z,y);
    t[z][0]=++tot;f[tot]=z;key[tot]=h[tot]=c;size[tot]=1;
    updata(z);updata(y);
}
void change(int x,int c) {
    int y=find(root,x);splay(y,0);
    key[y]=c;updata(y);
}
int query(int l,int r) {
    int x=find(root,l-1);splay(x,0);
    int y=find(root,r+1);splay(y,x);
    return h[t[y][0]];
}
int main() {
    scanf("%s",s+2);int n=strlen(s+2);
    mi[0]=1;fo(i,1,N-1) mi[i]=mi[i-1]*27%mo;
    root=build(1,n+2,0);tot=n+2;
    for(scanf("%d",&m);m;m--) {
        scanf("%s",s);
        if (s[0]=='Q') {
            scanf("%d%d",&x,&y);
            int l=1,r=min(tot-x,tot-y);x++;y++;
            while (l<r) {
                int mid=(l+r)/2;
                if (query(x,x+mid-1)==query(y,y+mid-1)) l=mid+1;
                else r=mid;
            }
            printf("%d\n",l-1);
        } else if (s[0]=='R') {
            scanf("%d%s",&x,ch);
            change(++x,ch[0]-'a');
        } else {
            scanf("%d%s",&x,ch);
            insert(++x,ch[0]-'a');
        }
    }
}