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 
#include 
#include 
#include 

#include 
#include 
#include 
#include 
#include 
#include 
#include 

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