bzoj 2434 [Noi2011]阿狸的打字机 AC自动机

[Noi2011]阿狸的打字机

Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 4001 Solved: 2198
[Submit][Status][Discuss]

Description

阿狸喜欢收藏各种稀奇古怪的东西，最近他淘到一台老式的打字机。打字机上只有28个按键，分别印有26个小写英文字母和'B'、'P'两个字母。

经阿狸研究发现，这个打字机是这样工作的：

l 输入小写字母，打字机的一个凹槽中会加入这个字母(这个字母加在凹槽的最后)。

l 按一下印有'B'的按键，打字机凹槽中最后一个字母会消失。

l 按一下印有'P'的按键，打字机会在纸上打印出凹槽中现有的所有字母并换行，但凹槽中的字母不会消失。

例如，阿狸输入aPaPBbP，纸上被打印的字符如下：

a

aa

ab

我们把纸上打印出来的字符串从1开始顺序编号，一直到n。打字机有一个非常有趣的功能，在打字机中暗藏一个带数字的小键盘，在小键盘上输入两个数(x,y)（其中1≤x,y≤n），打字机会显示第x个打印的字符串在第y个打印的字符串中出现了多少次。

阿狸发现了这个功能以后很兴奋，他想写个程序完成同样的功能，你能帮助他么？

Input

输入的第一行包含一个字符串，按阿狸的输入顺序给出所有阿狸输入的字符。

第二行包含一个整数m，表示询问个数。

接下来m行描述所有由小键盘输入的询问。其中第i行包含两个整数x, y，表示第i个询问为(x, y)。

Output

输出m行，其中第i行包含一个整数，表示第i个询问的答案。

Sample Input

aPaPBbP

3

1 2

1 3

2 3

Sample Output

2

1

0

HINT

1<=N<=10^5

1<=M<=10^5

输入总长<=10^5

首先构建ac自动机的fail树

得出dfs序，得出每个结点进出时间l[x]，r[x]，考虑这样一种暴力

对于一个询问x，y，查询自动机上root-y的每一个结点，沿着fail指针是否会走到x的结尾点

如果可以即答案+1

而在fail树中，变为查询自动机上root-y的所有结点中，有多少个在x的子树中？

只要在自动机上再重新走一遍，走到一个结点y，则将1-l[y]都+1，解决询问x，y（把y相同的链表串起来），即查询l[x]到r[x]的和。。。当遇到一个B时1-l[y]都-1

树状数组实现加减和区间求和

  1 #include<map>
  2 #include<set>
  3 #include<cmath>
  4 #include<queue>
  5 #include<cstdio>
  6 #include<vector>
  7 #include<cstring>
  8 #include<cstdlib>
  9 #include<iostream>
 10 using namespace std;
 11 int read()
 12 {
 13     int x=0;char ch=getchar();
 14     while(ch<'0'||ch>'9')ch=getchar();
 15     while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
 16     return x;
 17 }
 18 int n,m,cnt;
 19 int T,l[100005],r[100005],t[200005];
 20 int last[100005],lastq[100005],pos[100005],ans[100005];
 21 char ch[100005];
 22 struct fail_tree{
 23     int to,next,v;
 24 }e[100005],que[100005];
 25 
 26 void add(int x,int val)
 27 {
 28     for(int i=x;i<=T;i+=i&-i)t[i]+=val;
 29 }
 30 int query(int x)
 31 {
 32     int sum=0;
 33     for(int i=x;i;i-=i&-i)sum+=t[i];
 34     return sum;
 35 }
 36 void insert(int u,int v)
 37 {
 38     e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
 39 }
 40 struct acm{
 41     int cnt;
 42     int fa[100005],fail[100005],next[100005][26];
 43     int q[100005];
 44     acm(){
 45         cnt=1;
 46         for(int i=0;i<26;i++)
 47             next[0][i]=1;
 48     }
 49     void build(){
 50         int now=1,id=0;
 51         for(int i=0;i<n;i++)
 52             if(ch[i]=='P')pos[++id]=now;
 53             else if(ch[i]=='B')now=fa[now];
 54             else 
 55             {
 56                 if(!next[now][ch[i]-'a'])
 57                 {
 58                     next[now][ch[i]-'a']=++cnt;
 59                     fa[cnt]=now;
 60                 }
 61                 now=next[now][ch[i]-'a'];
 62             }
 63     }
 64     void buildfail(){
 65         int head=0,tail=1;
 66         fail[1]=0;q[0]=1;
 67         while(head!=tail)
 68         {
 69             int now=q[head];head++;
 70             for(int i=0;i<26;i++)
 71                 if(next[now][i])
 72                 {
 73                     int v=next[now][i];
 74                     int k=fail[now];
 75                     while(!next[k][i])k=fail[k];
 76                     fail[v]=next[k][i];
 77                     q[tail++]=v;
 78                 }
 79         }
 80     }
 81     void solve(){
 82         int now=1,id=0;
 83         add(l[1],1);
 84         for(int i=0;i<n;i++)
 85             if(ch[i]=='P')
 86             {
 87                 id++;
 88                 for(int x=lastq[id];x;x=que[x].next)
 89                 {
 90                     int t=pos[que[x].to];
 91                     ans[x]=query(r[t])-query(l[t]-1);
 92                 }
 93             }
 94             else if(ch[i]=='B')add(l[now],-1),now=fa[now];
 95             else now=next[now][ch[i]-'a'],add(l[now],1);;
 96     }    
 97 }acm;
 98 void dfs(int x)
 99 {
100     l[x]=++T;
101     for(int i=last[x];i;i=e[i].next)
102         dfs(e[i].to);
103     r[x]=++T;
104 }
105 int main()
106 {
107     scanf("%s",ch);
108     n=strlen(ch);
109     acm.build();
110     acm.buildfail();
111     for(int i=1;i<=acm.cnt;i++)
112         insert(acm.fail[i],i);
113     m=read();
114     for(int i=1;i<=m;i++)
115     {
116         int x=read(),y=read();
117         que[i].next=lastq[y];
118         lastq[y]=i;
119         que[i].to=x;
120     }
121     dfs(0);
122     acm.solve();
123     for(int i=1;i<=m;i++)printf("%d
",ans[i]);
124 }