NOI2011阿狸的打字机

题目背景
阿狸喜欢收藏各种稀奇古怪的东西，最近他淘到一台老式的打字机。

题目描述
打字机上只有28个按键，分别印有26个小写英文字母和’B’、’P’两个字母。经阿狸研究发现，这个打字机是这样工作的：
·输入小写字母，打字机的一个凹槽中会加入这个字母(这个字母加在凹槽的最后)。
·按一下印有’B’的按键，打字机凹槽中最后一个字母会消失。
·按一下印有’P’的按键，打字机会在纸上打印出凹槽中现有的所有字母并换行，但凹槽中的字母不会消失。
例如，阿狸输入aPaPBbP，纸上被打印的字符如下：
a aa ab 我们把纸上打印出来的字符串从1开始顺序编号，一直到n。打字机有一个非常有趣的功能，在打字机中暗藏一个带数字的小键盘，在小键盘上输入两个数(x,y)（其中1≤x,y≤n），打字机会显示第x个打印的字符串在第y个打印的字符串中出现了多少次。
阿狸发现了这个功能以后很兴奋，他想写个程序完成同样的功能，你能帮助他么？

输入输出格式

输入格式：
输入的第一行包含一个字符串，按阿狸的输入顺序给出所有阿狸输入的字符。
第二行包含一个整数m，表示询问个数。
接下来m行描述所有由小键盘输入的询问。其中第i行包含两个整数x, y，表示第i个询问为(x, y)。
输出格式：
输出m行，其中第i行包含一个整数，表示第i个询问的答案。

分析：
1.我们发现如果y串上的某一节点在fail边上可以走到x的结束节点，那么答案就会+1
2.考虑用AC自动机构造fail树，那么问题就转化为求以x为根的子树中有多少个是root到y上的节点
3.利用dfs序，每遍历到trie树上的一个点时，就把它的in[x]+1，遇到B时，它的in[x]就减1。
4.用离线做法，把y相同的询问一起考虑，每次到达一个y，就用sum(out[x])-sum(in[x]-1)更新答案。发现树状数组可以很好的解决。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
using namespace std;
const int maxn=100010;
char s[maxn];
vector<int>g[maxn],w[maxn];
int c[maxn];
int lowbit(int x){
    return x&(-x);
}
void Add(int x,int a){
    while(x<maxn){
        c[x]+=a;
        x+=lowbit(x);
    }
}
int sum(int x){
    int ret=0;
    while(x){
        ret+=c[x];
        x-=lowbit(x);
    }
    return ret;
}
int to[maxn],Next[maxn],Begin[maxn],e;
int len,cnt;
int in[maxn],out[maxn],fa[maxn],ans[maxn];
struct acm{
    int ch[maxn][26],pos[maxn],fail[maxn];
    void build(){
        int u=0,id=0;
        for(int i=1;i<=len;i++){
            if(s[i]=='B'){
                u=fa[u];
            }else if(s[i]=='P'){
                pos[++id]=u;
            }else{
                if(!ch[u][s[i]-'a']){
                    ch[u][s[i]-'a']=++cnt;
                    fa[cnt]=u;
                }
                u=ch[u][s[i]-'a'];
            }
        }
    }
    void build_fail(){
        queue<int>q;
        for(int i=0;i<26;i++)if(ch[0][i])q.push(ch[0][i]);
        while(!q.empty()){
            int u=q.front();q.pop();
            for(int i=0;i<26;i++){
                int& c=ch[u][i],o=ch[fail[u]][i];
                if(c){
                    q.push(c);
                    fail[c]=o;
                }else c=o;
            }
        }
    }
    void solve(){
        int id=0,u=0;
        for(int i=1;i<=len;i++){
            if(s[i]=='P'){
                id++;
                for(int j=0;j<g[id].size();j++){
                    int v=pos[g[id][j]];
                    ans[w[id][j]]=sum(out[v])-sum(in[v]-1);
                }
            }else if(s[i]=='B') Add(in[u],-1),u=fa[u];
            else u=ch[u][s[i]-'a'],Add(in[u],1);
        }
    }
}AC;
int m;
void add(int x,int y){
    to[++e]=y;
    Next[e]=Begin[x];
    Begin[x]=e;
}
int dfn;
void dfs(int x){
    in[x]=++dfn;
    for(int i=Begin[x];i;i=Next[i]) dfs(to[i]);
    out[x]=dfn;
}
int main(){
    scanf("%s",s+1);
    len=strlen(s+1);
    AC.build();
    AC.build_fail();
    for(int i=1;i<=cnt;i++) add(AC.fail[i],i);
    scanf("%d",&m);
    for(int i=1;i<=m;i++){
        int x,y;
        scanf("%d%d",&x,&y);
        g[y].push_back(x);
        w[y].push_back(i);
    }
    dfs(0);
    AC.solve();
    for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
    return 0;
}

^_^