bzoj2555 SubString（SAM+lct）

支持在末尾加入一个字符串
查询一个字符串的出现次数
强制在线

我们建出SAM，找到这个字符串，如果找不到直接返回0.这个节点的|Right|就是答案。因此我们需要动态维护|Right|.考虑我们构造的过程，把np连到parent树上之后，np到根的路径上的所有点都+1了。而用nq代替q时，需要把v[nq]变成v[q]，然后还是np到根的路径上的所有点+1.

我们可以用有根树lct来维护这个路径加法。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 1200010
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=0,m,son[N][26],mask=0,mx[N],par[N],rt,last;
int fa[N],c[N][2],v[N],tag[N],q[N];
char s[3000010];
inline bool isroot(int x){return x!=c[fa[x]][0]&&x!=c[fa[x]][1];}
inline void add(int p,int val){
    if(!p) return;tag[p]+=val;v[p]+=val;
}
inline void pushdown(int p){
    if(!tag[p]) return;
    add(c[p][0],tag[p]);add(c[p][1],tag[p]);tag[p]=0;
}
inline void rotate(int x){
    int y=fa[x],z=fa[y],l=x==c[y][1],r=l^1;
    if(!isroot(y)) 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;
}
inline void splay(int x){
    int top=0;q[++top]=x;
    for(int xx=x;!isroot(xx);xx=fa[xx]) q[++top]=fa[xx];
    while(top) pushdown(q[top--]);
    while(!isroot(x)){
        int y=fa[x],z=fa[y];
        if(!isroot(y)){
            if(x==c[y][1]^y==c[z][1]) rotate(x);
            else rotate(y);
        }rotate(x);
    }
}
inline void access(int x){
    int y=0;while(x){splay(x);c[x][1]=y;y=x;x=fa[x];}
}
inline void link(int x,int y){
    fa[x]=y;
}
inline void cut(int x){
    access(x);splay(x);fa[c[x][0]]=0;c[x][0]=0;
}
inline void update(int p){
    access(p);splay(p);add(p,1);
}
inline void ins(int ch){
    int p=last,np=++n;last=np;mx[np]=mx[p]+1;
    for(;p&&!son[p][ch];p=par[p]) son[p][ch]=np;
    if(!p){par[np]=rt;link(np,rt);update(np);return;}
    int q=son[p][ch];
    if(mx[q]==mx[p]+1){par[np]=q;link(np,q);update(np);return;}
    int nq=++n;memcpy(son[nq],son[q],sizeof(son[q]));cut(q);
    mx[nq]=mx[p]+1;v[nq]=v[q];par[nq]=par[q];link(nq,par[q]);
    par[q]=par[np]=nq;link(q,nq);link(np,nq);
    for(;p&&son[p][ch]==q;p=par[p]) son[p][ch]=nq;update(np);
}
inline void decode(int mask){
    for(int i=1;i<=m;++i){
        mask=(mask*131+i-1)%m;
        swap(s[i],s[mask+1]);
    }
}
inline int query(){
    int p=rt;
    for(int i=1;i<=m;++i){
        if(!son[p][s[i]-'A']) return 0;
        p=son[p][s[i]-'A'];
    }splay(p);return v[p];
}
int main(){
//  freopen("a.in","r",stdin);
    int owo=read();scanf("%s",s+1);m=strlen(s+1);
    last=rt=++n;for(int i=1;i<=m;++i) ins(s[i]-'A');
    while(owo--){
        char op[10];scanf("%s%s",op+1,s+1);m=strlen(s+1);decode(mask);
        if(op[1]=='A'){
            for(int i=1;i<=m;++i) ins(s[i]-'A');continue;
        }int res=query();mask^=res;printf("%d\n",res);
    }return 0;
}