[Noi2011]阿狸的打字机
Time Limit: 10 Sec Memory Limit: 256 MB[Submit][Status][Discuss]
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
3
1 2
1 3
2 3
Sample Output
2
1
0
1
0
HINT
1<=N<=10^5
1<=M<=10^5
输入总长<=10^5
题目相当于问在fail树中，字符串x最后一个节点的子树内 有多少属于 字符串y的节点
暴力法：对于字符串y的每个节点，沿fail指针一直往上走，若能走到字符串x的最后一个节点，ans++
据说能得70分
AC做法：离线
1、边输入字符串边建立AC自动机，如果是‘B’，now就跳到trie中父节点
2、抽离fail树构图
3、dfs 遍历fail树，得到每个点的入栈时间和出栈时间（即子树的dfs序范围）
4、如果询问x，y，则在字符串y的最后一个字母 上 挂上 x的最后一个字母
这样遍历到y的时候，就可以处理与 某字符串 在y中出现多少次的所有询问
5、在AC自动机上再跑一遍字符串，
如果是小写字母，当前dfs序位置+1，到往下走；
‘B’，当前dfs序位置-1，到父节点
‘P’，统计答案，做法：枚举这个y的每一个x，对于这一对x、y，答案就是x的入栈时间——出栈时间有多少个1
统计一段区间1的个数，支持单点修改操作：树状数组
#include<queue> #include<cstdio> #include<cstring> #define N 100101 using namespace std; int tot,trie[N][26],f[N]; int fa[N]; int fail_front[N],fail_nxt[N],fail_to[N],cnt; int tree_front[N],tree_nxt[N],tree_to[N]; char s[N]; int len,now,sum_p; int p[N]; int c[N]; int ans[N],in[N],out[N]; queue<int>q; void add_tree(int u,int v,int k) { tree_to[++tot]=v; tree_nxt[tot]=tree_front[u]; tree_front[u]=tot; ans[tot]=k; } void fail_add(int u,int v) { fail_to[++tot]=v; fail_nxt[tot]=fail_front[u]; fail_front[u]=tot; } int lowbit(int x) { return x&(-x); } void add(int x,int w) { for(int i=x;i<=cnt;i+=lowbit(i)) c[i]+=w; } int query(int x) { int ans=0; for(int i=x;i;i-=lowbit(i)) ans+=c[i]; return ans; } void dfs(int x) { in[x]=++cnt; for(int i=fail_front[x];i;i=fail_nxt[i]) dfs(fail_to[i]); out[x]=cnt; } void getfail() { for(int i=0;i<26;i++) trie[0][i]=1; q.push(1); int j; while(!q.empty()) { now=q.front(); q.pop(); for(int i=0;i<26;i++) { if(!trie[now][i]) { trie[now][i]=trie[f[now]][i]; continue; } q.push(trie[now][i]); j=f[now]; f[trie[now][i]]=trie[j][i]; fail_add(trie[j][i],trie[now][i]); } } } int main() { scanf("%s",s); len=strlen(s); now=1; tot=1; for(int i=0;i<len;i++) { if(s[i]>='a') { trie[now][s[i]-'a']=++tot; fa[tot]=now; now=trie[now][s[i]-'a']; } else if(s[i]=='B') now=fa[now]; else p[++sum_p]=now; } tot=0; getfail(); dfs(1); tot=0; int m,x,y; scanf("%d",&m); for(int i=1;i<=m;i++) { scanf("%d%d",&x,&y); add_tree(p[y],p[x],i); } now=1; for(int i=0;i<len;i++) { if(s[i]=='P') { for(int i=tree_front[now];i;i=tree_nxt[i]) ans[i]=query(out[tree_to[i]])-query(in[tree_to[i]]-1); } else if(s[i]=='B') { add(in[now],-1); now=fa[now]; } else { add(in[trie[now][s[i]-'a']],1); now=trie[now][s[i]-'a']; } } for(int i=1;i<=m;i++) printf("%d ",ans[i]); }