BZOJ 2434 [Noi2011]阿狸的打字机

Description

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

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

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

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

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

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

a

aa

ab

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

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

【题目分析】

AC自动机吼题啊。

利用trie图fail指针的含义（该点的后缀的下一个位置），然后将它反向，就成为了统计子树大小的题目了，用自动机DFS序来加入、删除点，然后用反向图的DFS序来维护。

【代码】

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>

#include <set>
#include <map>
#include <string>
#include <algorithm>
#include <vector>
#include <iostream>
#include <queue>

using namespace std;

#define maxn 100005
#define inf (0x3f3f3f3f)

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 trie[maxn][26],fail[maxn],cnt=0,len,sta[maxn],top=0,tot=1;
char s[maxn];
int id[maxn],in[maxn],out[maxn],idx=0;
int t[maxn],pos[maxn];
int h[maxn<<1],to[maxn<<1],ne[maxn<<1],en=0;
vector <int> que[maxn];
vector <int> bel[maxn];
int ans[maxn];

void addedge(int a,int b)
{
//  printf("add an edge %d to %d\n",a,b);
to[en]=b;
ne[en]=h[a];
h[a]=en++;
}

void add(int x,int f)
{
//  printf("%d add %d\n",x,f);
for (;x<=maxn;x+=x&(-x)) t[x]+=f;
}

int gs(int x)
{
int ret=0;
for (;x;x-=x&(-x)) ret+=t[x];
return ret;
}

void getfail()
{
queue <int> q;
while (!q.empty()) q.pop();
q.push(1);
while (!q.empty())
{
int x=q.front(); q.pop();
//      printf("%d\n",x);
for (int i=0;i<26;++i)
if (trie[x][i])
{
int j=fail[x];
while (j&&!trie[j][i]) j=fail[j];
fail[trie[x][i]]=j?trie[j][i]:1;
if(j==0) fail[trie[x][i]]=1;
q.push(trie[x][i]);
}
}
}

void dfs(int k)
{
in[k]=++idx;
for (int i=h[k];i>=0;i=ne[i])
dfs(to[i]);
out[k]=idx;
}

void solve(int k)
{
add(in[k],1);
for (int i=0;i<que[id[k]].size();++i)
{
ans[bel[id[k]][i]]=gs(out[pos[que[id[k]][i]]])-gs(in[pos[que[id[k]][i]]]-1);
//      printf("ans %d from %d to %d\n",bel[id[k]][i],out[pos[que[id[k]][i]]],in[pos[que[id[k]][i]]]);
}
for (int i=0;i<26;++i) if (trie[k][i]) solve(trie[k][i]);
add(in[k],-1);
}

int main()
{
memset(h,-1,sizeof h);
scanf("%s",s+1); len=strlen(s+1);
int now=1;sta[0]=1;
for (int i=1;i<=len;++i)
{
if (s[i]=='P') id[now]=++cnt,pos[cnt]=now;
else if (s[i]=='B') now=sta[--top];
else
{
if (!trie[now][s[i]-'a']) trie[now][s[i]-'a']=++tot;
now=trie[now][s[i]-'a'];
sta[++top]=now;
}
//      printf("now is %d\n",now);
//      printf("sta is :");
//      for (int i=1;i<=top;++i) printf("%d ",sta[i]);
//      printf("\n");
}
getfail();
for (int i=2;i<=tot;++i) addedge(fail[i],i);
//  printf("pos is :");
//  for (int i=1;i<=tot;++i) printf("%d ",pos[i]); printf("\n");
//  printf("\n");
dfs(1);
int m=read();
for (int i=1;i<=m;++i)
{
int x,y;
scanf("%d%d",&x,&y);
que[y].push_back(x);
bel[y].push_back(i);
}
solve(1);
for (int i=1;i<=m;++i) printf("%d\n",ans[i]);
}