【BZOJ2555】SubString（后缀自动机，LCT）

题意：给你一个字符串init，要求你支持两个操作

(1):在当前字符串的后面插入一个字符串

(2):询问字符串s在当前字符串中出现了几次？(作为连续子串)

你必须在线支持这些操作。

长度 <= 600000，询问次数<= 10000,询问总长度<= 3000000

思路：因为有加边，删边，加点操作，需要动态维护SAM中每个right集合的大小，所以使用LCT维护

需要维护根节点=1号点到每个结点路径上的和，因为固定了1号点为根节点所以不需要makeroot和rev标记

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long ll;
 typedef unsigned int uint;
 typedef unsigned long long ull;
 typedef pair<int,int> PII;
 typedef pair<ll,ll> Pll;
 typedef vector<int> VI;
 typedef vector<PII> VII;
 typedef pair<ll,int>P;
 #define N  1200010
 #define M  210000
 #define fi first
 #define se second
 #define MP make_pair
 #define pi acos(-1)
 #define mem(a,b) memset(a,b,sizeof(a))
 #define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
 #define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
 #define lowbit(x) x&(-x)
 #define Rand (rand()*(1<<16)+rand())
 #define id(x) ((x)<=B?(x):m-n/(x)+1)
 #define ls p<<1
 #define rs p<<1|1

 const int MOD=1e9+,inv2=(MOD+)/;
       double eps=1e-;
       int INF=<<;
       ll inf=5e13;
       int dx[]={-,,,};
       int dy[]={,,-,};

 int mask;
 char s[];
 string chars;

 int read()
 {
    int v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 void gets(int mask)
 {
     scanf("%s",s);
     chars=s;
     for(int j=;j<chars.length();j++)
     {
         mask=(mask*+j)%chars.length();
         char t=chars[j];
         chars[j]=chars[mask];
         chars[mask]=t;
     }
 }

 struct lct
 {
     int t[N][],w[N],fa[N],q[N],rev[N],tag[N],top;

     void add(int x,int y)
     {
         if(x)
         {
             w[x]+=y;
             tag[x]+=y;
         }
     }

     int isroot(int x)
     {
         return ((t[fa[x]][]!=x)&&(t[fa[x]][]!=x));
     }

     void pushdown(int x)
     {
         int l=t[x][],r=t[x][];
         if(tag[x])
         {
             add(l,tag[x]);
             add(r,tag[x]);
             tag[x]=;
         }
     }

     void rotate(int x)
     {
         int y=fa[x],z=fa[y];
         int l=(t[y][]==x),r=l^;
         if(!isroot(y)) t[z][t[z][]==y]=x;
         fa[t[x][r]]=y,fa[y]=x,fa[x]=z;
         t[y][l]=t[x][r],t[x][r]=y;
     }

     void splay(int x)
     {
         top=;
         q[++top]=x;
         for(int i=x;!isroot(i);i=fa[i]) q[++top]=fa[i];
         while(top) pushdown(q[top--]);
         while(!isroot(x))
         {
             int y=fa[x],z=fa[y];
             if(!isroot(y))
             {
                 if((t[y][]==x)^(t[z][]==y)) rotate(x);
                  else rotate(y);
             }
             rotate(x);
         }
     }

     void access(int x)
     {
         for(int k=;x;k=x,x=fa[x])
         {
             splay(x);
             t[x][]=k;
         }
     }

     void link(int x,int y)
     {
         fa[x]=y;
         access(y);
         splay(y);
         add(y,w[x]);
     }

     void cut(int x)
     {
         access(x);
         splay(x);
         add(t[x][],-w[x]);
         fa[t[x][]]=;
         t[x][]=;
     }

 }lct;

 struct sam
 {
     int cnt;
     int F[N],ch[N][];
     int st[N],b[N],bl[N],c[N];
     int p,np,q,nq;

     sam()
     {
         cnt=np=;
     }

     void extend(int x)
     {
         p=np; st[np=++cnt]=st[p]+;
         while(p&&!ch[p][x])
         {
             ch[p][x]=np;
             p=F[p];
         }
         lct.w[np]=;
         if(!p) F[np]=,lct.link(np,);
          else if(st[p]+==st[q=ch[p][x]]) F[np]=q,lct.link(np,q);
           else
           {
               st[nq=++cnt]=st[p]+;
               memcpy(ch[nq],ch[q],sizeof ch[q]);
               F[nq]=F[q];
               lct.link(nq,F[q]);
               F[np]=F[q]=nq;
               lct.cut(q);
               lct.link(q,nq);
               lct.link(np,nq);
               while(ch[p][x]==q)
               {
                     ch[p][x]=nq;
                     p=F[p];
               }
           }
     }

     void build()
     {
         scanf("%s",s);
         int n=strlen(s);
         rep(i,,n-) extend(s[i]-'A');
     }

     void add()
     {
         gets(mask);
         int n=chars.length();
         rep(i,,n-) extend(chars[i]-'A');
     }

     int query()
     {
         gets(mask);
         int p=,n=chars.length();
         rep(i,,n-)
         {
             p=ch[p][chars[i]-'A'];
             if(!p) return ;
         }
         lct.splay(p);
         return lct.w[p];
     }

 }sam;

 int main()
 {
     int Q;
     scanf("%d",&Q);
     sam.build();
     while(Q--)
     {
         scanf("%s",s);
         if(s[]=='A') sam.add();
          else
          {
             int ans=sam.query();
             printf("%d\n",ans);
             mask^=ans;
          }
     }
     return ;
 }