bzoj 2555(linkcuttree+后缀自动机)

2555: SubString

Time Limit: 30 Sec   Memory Limit: 512 MB
Submit: 1702   Solved: 491
[ Submit][ Status][ Discuss]

Description

  
    懒得写背景了，给你一个字符串init，要求你支持两个操作
    
    (1):在当前字符串的后面插入一个字符串
    
    (2):询问字符串s在当前字符串中出现了几次？(作为连续子串)
    
    你必须在线支持这些操作。
    

Input

    第一行一个数Q表示操作个数
    
    第二行一个字符串表示初始字符串init
    
    接下来Q行，每行2个字符串Type,Str 
    
    Type是ADD的话表示在后面插入字符串。
    
    Type是QUERY的话表示询问某字符串在当前字符串中出现了几次。
    
    为了体现在线操作，你需要维护一个变量mask，初始值为0
   
    
    读入串Str之后，使用这个过程将之解码成真正询问的串TrueStr。
    询问的时候，对TrueStr询问后输出一行答案Result
    然后mask = mask xor Result  
    插入的时候，将TrueStr插到当前字符串后面即可。

HINT:ADD和QUERY操作的字符串都需要解压
   

Output

Sample Input

2

A

QUERY B

ADD BBABBBBAAB

Sample Output

0

HINT

 40 % 的数据字符串最终长度 <= 20000，询问次数<= 1000，询问总长度<= 10000

    

100 % 的数据字符串最终长度 <= 600000，询问次数<= 10000,询问总长度<= 3000000

新加数据一组--2015.05.20

解题思路：后缀自动机，但是如果暴力更新right会超时，然后就想到用linkcuttree更新

各种细节，坑死。。。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
int qg,mask,cnt,last;
char c[600011],s[600011];
string chars;
int pre[3100000],ch[3100000][27],mx[3100000];
int fa[3100000],l[3100000],r[3100000],zhi[3100000],add[3100000];
int q[3100000];


inline int read()
{
char y; int x=0,f=1; y=getchar();
while (y<'0' || y>'9') {if (y=='-') f=-1; y=getchar();}
while (y>='0' && y<='9') {x=x*10+int(y)-48; y=getchar();}
return x*f;
}


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


bool isroot(int now)
 {
  if (l[fa[now]]==now || r[fa[now]]==now) return false;else return true;
 }


void pushdown(int now)
 {
  if (add[now])
  {
  if (l[now]!=0) add[l[now]]+=add[now],zhi[l[now]]+=add[now];
  if (r[now]!=0) add[r[now]]+=add[now],zhi[r[now]]+=add[now];
  add[now]=0;
 }
 }


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


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


void access(int now)
 {
  splay(now);
  while (fa[now]!=0)
  {
   splay(fa[now]); r[now]=0;
   r[fa[now]]=now; 
   splay(now);
}
r[now]=0; 
 }


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


void del(int x,int y)
 {
  access(y); splay(y); 
zhi[l[y]]-=zhi[y];add[l[y]]-=zhi[y];
    fa[l[y]]=0; l[y]=0;
 }


void extend(int x)
 {
  int p=last; last=++cnt; int np=cnt; zhi[cnt]=1; mx[np]=mx[p]+1;
  while (!ch[p][x] && p!=0) ch[p][x]=np,p=pre[p];
  if (!p) pre[np]=1,link(1,np);else
  {
  int q=ch[p][x];
  if (mx[q]==mx[p]+1)
  {
    pre[np]=q; link(q,np);
}else
   {
    int nq=++cnt; mx[nq]=mx[p]+1; del(pre[q],q); link(nq,q); link(nq,np); link(pre[q],nq); 
    pre[nq]=pre[q]; 
       pre[q]=pre[np]=nq; 
       memcpy(ch[nq],ch[q],sizeof(ch[q]));
       while (p!=0 && ch[p][x]==q) ch[p][x]=nq,p=pre[p];
}
 }
 }


int main()
{
mask=0;
qg=read();
scanf("%s",c);
cnt=1; last=cnt;
    for (int i=1;i<=strlen(c);++i)
     {
      extend(c[i-1]-'A');
}
for (int i=1;i<=qg;++i)
{
scanf("%s",c); gets(mask);
if (c[0]=='A')
{
for (int j=1;j<=chars.length();++j)
{
extend(chars[j-1]-'A');
}
 }else
  {
  int now=1,g; int j;
  for (j=1;j<=chars.length();++j) 
   {
    g=chars[j-1]-'A';
    if (!ch[now][g]) break;
    now=ch[now][g];
 }
if (!ch[now][g] && j<=chars.length()) 
{
mask^=0; printf("0\n");
}else
 {
  splay(now); printf("%d\n",zhi[now]);
  mask=mask^zhi[now];
 }
  }
}
}