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

时间:2022-08-31 21:57:26

[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

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

 #include<map>
#include<set>
#include<cmath>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
using namespace std;
int read()
{
int x=;char ch=getchar();
while(ch<''||ch>'')ch=getchar();
while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
return x;
}
int n,m,cnt;
int T,l[],r[],t[];
int last[],lastq[],pos[],ans[];
char ch[];
struct fail_tree{
int to,next,v;
}e[],que[]; void add(int x,int val)
{
for(int i=x;i<=T;i+=i&-i)t[i]+=val;
}
int query(int x)
{
int sum=;
for(int i=x;i;i-=i&-i)sum+=t[i];
return sum;
}
void insert(int u,int v)
{
e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
}
struct acm{
int cnt;
int fa[],fail[],next[][];
int q[];
acm(){
cnt=;
for(int i=;i<;i++)
next[][i]=;
}
void build(){
int now=,id=;
for(int i=;i<n;i++)
if(ch[i]=='P')pos[++id]=now;
else if(ch[i]=='B')now=fa[now];
else
{
if(!next[now][ch[i]-'a'])
{
next[now][ch[i]-'a']=++cnt;
fa[cnt]=now;
}
now=next[now][ch[i]-'a'];
}
}
void buildfail(){
int head=,tail=;
fail[]=;q[]=;
while(head!=tail)
{
int now=q[head];head++;
for(int i=;i<;i++)
if(next[now][i])
{
int v=next[now][i];
int k=fail[now];
while(!next[k][i])k=fail[k];
fail[v]=next[k][i];
q[tail++]=v;
}
}
}
void solve(){
int now=,id=;
add(l[],);
for(int i=;i<n;i++)
if(ch[i]=='P')
{
id++;
for(int x=lastq[id];x;x=que[x].next)
{
int t=pos[que[x].to];
ans[x]=query(r[t])-query(l[t]-);
}
}
else if(ch[i]=='B')add(l[now],-),now=fa[now];
else now=next[now][ch[i]-'a'],add(l[now],);;
}
}acm;
void dfs(int x)
{
l[x]=++T;
for(int i=last[x];i;i=e[i].next)
dfs(e[i].to);
r[x]=++T;
}
int main()
{
scanf("%s",ch);
n=strlen(ch);
acm.build();
acm.buildfail();
for(int i=;i<=acm.cnt;i++)
insert(acm.fail[i],i);
m=read();
for(int i=;i<=m;i++)
{
int x=read(),y=read();
que[i].next=lastq[y];
lastq[y]=i;
que[i].to=x;
}
dfs();
acm.solve();
for(int i=;i<=m;i++)printf("%d\n",ans[i]);
}