[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:

aPaPBbP 3 1 2 1 3 2 3 输出样例#1: 2 1 0 说明

数据范围:

对于100%的数据,n<=100000,m<=100000,第一行总长度<=100000。

题解:

这道题看着题解搞了一上午,但也学到了一些有关AC自动机的神奇的东西。

对于AC自动机中的每一个节点,如果节点A的fail指向节点B,则B对应的字符串一定在A对应的字符串中出现

所以这道题就变成了在B串中找谁的fail指针指向了A串的最后一位

这样,我们可以将这些fail指针化成边

由now的fail指针指向的点向now连边

这样构建出的东西好像叫 Fail树

这样本题就成了查询以A串的最后一位为根的子树中有多少个属于b串的节点,即为答案

但是这道题的询问太多辣,如果在线查询的话肯定会T掉

于是我们可以考虑离线处理

先读入所有的询问,然后按照y排个序

处理出所有节点的dfs序

那么这样这棵fail树就被我们压成了一个序列

我们都知道一个点的子树的dfs序肯定是连续的一段

所以我们可以考虑用树状数组维护

从0开始dfs一遍

并把这个点的值+1

如果这个点是一段字符串的结尾

就处理以这个点为y的所有询问:直接查询这个询问的x的子树(id[x]~id[x]+size[x]-1/dfn[x]~low[x])的值

回溯的时候再减掉就行辣


#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
const int M = 1000005 ;
const int N = 100005 ;
using namespace std;
inline int read(){
	char c=getchar(); int x=0,w=1;
	while(c>'9'||c<'0'){ if(c=='-') w=-1 ; c=getchar();}
	while(c>='0'&&c<='9'){ x=x*10+c-'0'; c=getchar(); }
	return x*w;
}
struct Trie{
	int son[26];
	int ch[26];
	int Fail;
	int father;
	int End ;
}t[M];
struct E{
	int nex,to;
}edge[M<<1];
struct Ques{
	int x,y;
	int id,Ans;
	friend inline bool operator < (Ques a , Ques b){
		return a.y<b.y;
	}
}p[M];
int num=0,hea[M];
inline void add_edge(int from,int to){
	edge[++num].nex=hea[from];
	edge[num].to=to;
	hea[from]=num;
}
int Ans[M];
int len,tot,n,m,tail[M];
char s[N];
int dfn[M],low[M],tim;
int l[M],r[M];
inline void Get_Fail(){
	queue<int>q;
	for(int i=0;i<=25;i++)
	  if(t[0].son[i]){
	  	t[t[0].son[i]].Fail=0;
	  	q.push(t[0].son[i]);
	  }
	while(!q.empty()){
		int u=q.front(); q.pop();
		for(int i=0;i<=25;i++){
			if(t[u].son[i]){
				t[t[u].son[i]].Fail=t[t[u].Fail].son[i];
				q.push(t[u].son[i]);
			}
			else t[u].son[i]=t[t[u].Fail].son[i];
		}
	}
}
struct BIT{
	int sum[M];
	inline int lowbit(int x){
		return x&(-x);
	}
	inline void change(int k,int x){
		while(k<=tim){
			sum[k]+=x;
			k+=lowbit(k);
		}
	}
	inline int query(int k){
		int temp=0;
		while(k){
			temp+=sum[k];
			k-=lowbit(k);
		}
		return temp;
	}
}T;
inline void dfs1(int u){
	dfn[u]=++tim;
	for(int i=hea[u];i;i=edge[i].nex){
		int v=edge[i].to;
		dfs1(v);
	}
	low[u]=tim;
}
void Dfs(int u){
	T.change(dfn[u],1);
	if(t[u].End)
	  for(int i=l[t[u].End];i<=r[t[u].End];i++)
		p[i].Ans=T.query(low[tail[p[i].x]])-T.query(dfn[tail[p[i].x]]-1);			
    for(int i=0;i<=25;i++)
      if(t[u].ch[i])
        Dfs(t[u].ch[i]);
    T.change(dfn[u],-1);
}
int main(){
	scanf("%s",s+1);
	int now=0;
	len=strlen(s+1);
	for(int i=1;i<=len;i++){
		if(s[i]>='a'&&s[i]<='z'){
			int v=s[i]-'a';
			if(!t[now].son[v]){
				t[now].son[v]=++tot;
				t[tot].father=now;
				t[now].ch[v]=tot;
			}
			now=t[now].son[v];
		}
		else if(s[i]=='B') now=t[now].father;
		else if(s[i]=='P') {
			tail[++n]=now;
			t[now].End=n;
		}
	}
	m=read();
	for(int i=1;i<=m;i++){
		p[i].x=read();
		p[i].y=read();
		p[i].id=i;
	}
	sort(p+1,p+m+1);
	t[0].Fail=0;
	Get_Fail();
	for(int i=1;i<=tot;i++)
		add_edge(t[i].Fail,i);
	dfs1(0);
	int tmp=0;
	for(int i=1;i<=m;i=tmp){
		l[p[i].y]=i;
		while(p[i].y==p[tmp].y) ++tmp;
	    r[p[i].y]=tmp-1;
	}
	Dfs(0);
	for(int i=1;i<=m;i++)
	  Ans[p[i].id]=p[i].Ans;
	for(int i=1;i<=m;i++)
	  printf("%d\n",Ans[i]);
	return 0;
}