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]);}