[bzoj2434][AC自动机][树状数组]阿狸的打字机

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树
对于每个节点，我们把他和他的fail连起来，那么这样就会形成一棵fail树
如果对于一个串y的在Tire上的某个节点，让他去跳fail树的祖先，能跳到x这个串的末尾，那么可以说明x在y的这个位置是出现了的
那么我们可以先构建出fail树，然后每次询问x的末尾这个节点的子树中，有多少来自y的节点
直接做不好做，考虑离线
我们可以先构出dfs序，因为子树问题一般均需要dfs序解决
然后深搜Tire，进入一个节点时把这个点在Tire中的权加1，离开时减1消去影响。当询问到了一个串的末尾时，我们可以通过一个链表找到关于他的询问，然后树状数组统计答案回答即可

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
struct node{int x,y,next;}a[210000];int len,last[110000];
void ins(int x,int y){len++;a[len].x=x;a[len].y=y;a[len].next=last[x];last[x]=len;}
struct ACmachine
{
    int son[30],bel,fail,f;
    ACmachine(){memset(son,-1,sizeof(son));bel=-1;}
}tr[210000];int tot,root;
struct ask{int x,y,op;}w[110000];
bool cmp(ask n1,ask n2){return n1.y<n2.y;}
int las[110000],pos[110000];
char ch[210000];
int pre[210000],cnt;
int belong[210000];//第x个打印的字符串的结尾在belong[x] 
void add()
{
    int len=strlen(ch+1);int p=root,cnt=0;
    memset(pre,-1,sizeof(pre));
    for(int i=1;i<=len;i++)
    {
        if(ch[i]=='B')p=tr[p].f;
        else if(ch[i]=='P')
        {
            cnt++;pre[cnt]=tr[p].bel;
            tr[p].bel=cnt;belong[cnt]=p;
        }
        else
        {
            int y=ch[i]-'a'+1;
            if(tr[p].son[y]==-1)tr[p].son[y]=++tot;
            tr[tr[p].son[y]].f=p;
            p=tr[p].son[y];
        }
    }
}
queue<int> q;
void buildfail()
{
    q.push(0);
    while(!q.empty())
    {
        int x=q.front();
        for(int i=1;i<=26;i++)
        {
            int p=tr[x].son[i];if(p==-1)continue;
            if(x==0)tr[p].fail=0;
            else
            {
                int j=tr[x].fail;
                while(j && tr[j].son[i]==-1)j=tr[j].fail;
                tr[p].fail=max(tr[j].son[i],0);
            }
            ins(tr[p].fail,p);
            q.push(p);
        }
        q.pop();
    }
}
int in[110000],ot[110000],dfn;
void pre_tree_node(int x)
{
    in[x]=++dfn;
    for(int k=last[x];k;k=a[k].next)
        pre_tree_node(a[k].y);
    ot[x]=dfn;
}
int s[110000];
int lowbit(int x){return x&-x;}
void change(int x,int c){while(x<=dfn){s[x]+=c;x+=lowbit(x);}}
int findsum(int x){int ret=0;while(x>=1)ret+=s[x],x-=lowbit(x);return ret;}
int m,answer[110000];
void sol(int x)
{
    change(in[x],1);
    if(tr[x].bel!=-1)
    {
        int p=tr[x].bel;
        while(p!=-1)
        {
            int u=pos[p];
            while(u!=-1)
            {
                answer[u]=findsum(ot[belong[w[u].x]])-findsum(in[belong[w[u].x]]-1);
                u=las[u];
            }
            p=pre[p];
        }
    }
    for(int i=1;i<=26;i++)
        if(tr[x].son[i]!=-1)sol(tr[x].son[i]);
    change(in[x],-1);
}
int main()
{
    scanf("%s",ch+1);
    add();
    buildfail();
    scanf("%d",&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&w[i].x,&w[i].y);
        w[i].op=i;
    }
    memset(las,-1,sizeof(las));memset(pos,-1,sizeof(pos));
    for(int i=1;i<=m;i++)las[i]=pos[w[i].y],pos[w[i].y]=i;//第i个询问前一个询问是什么 
    pre_tree_node(0);
    sol(root);
    for(int i=1;i<=m;i++)printf("%d\n",answer[i]);
    return 0;
}