【bzoj2434】[Noi2011]阿狸的打字机 AC自动机+fail树+dfs序+树状数组

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

Source

Trie

竟然过样例就一遍A了…

学了fail树发现这题不难啊…细节处理稍麻烦

建好AC自动机，建好fail树，要求x串在y串中出现次数，可以离线一下，把y串上每个点权值+1，然后查询x串的子树权值和…

因为权值是动态的，还要维护子树，可以考虑dfs序，由于每次修改一个字符，所以树状数组维护即可。

#include<cstring>
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;

const int SZ = 1000010;

int ch[SZ][30],sz = 0,val[SZ],fa[SZ];
int num[SZ],id = 0;

void insert(char s[])
{
int p = 0;
int l = strlen(s);
for(int i = 0;i < l;i ++)
{
if(s[i] == 'B') p = fa[p];
else if(s[i] == 'P')
num[++ id] = p;
else
{
int c = s[i] - 'a' + 1;
if(!ch[p][c]) ch[p][c] = ++ sz;
fa[sz] = p;
p = ch[p][c];
}
}
}

int head[SZ],nxt[SZ],to[SZ],tot = 0;

void build(int f,int t)
{
to[++ tot] = t;
nxt[tot] = head[f];
head[f] = tot;
}

int fail[SZ];

queue<int> q;

void build_ac()
{
fail[0] = 0;
for(int i = 1;i <= 26;i ++)
{
int u = ch[0][i];
if(u) { q.push(u); fail[u] = 0; }
}
while(q.size())
{
int f = q.front(); q.pop();
for(int i = 1;i <= 26;i ++)
{
int u = ch[f][i];
if(!u) continue;
q.push(u);
int v = fail[f];
while(v && !ch[v][i]) v = fail[v];
fail[u] = ch[v][i];
}
}
}

int dfs_clock = 0,pre[SZ],suf[SZ];
int dfs_num[SZ];

void dfs(int u)
{
dfs_num[++ dfs_clock] = u;
pre[u] = dfs_clock;
for(int i = head[u];i;i = nxt[i])
{
dfs(to[i]);
}
dfs_num[++ dfs_clock] = u;
suf[u] = dfs_clock;
}

int bits[SZ];

void add(int i,int d)
{
while(i < dfs_clock)
{
bits[i] += d;
i += i & -i;
}
}

int sum(int i)
{
int ans = 0;
while(i)
{
ans += bits[i];
i -= i & -i;
}
return ans;
}

struct ques{
int x,y,id,ans;
}ask[SZ];

bool cmp1(ques a,ques b) { return a.y < b.y; }
bool cmp2(ques a,ques b) { return a.id < b.id; }

char s[SZ];

void solve()
{
int l = strlen(s);
int now = 0;
int p = 0;
for(int i = 0,j = 1;i < l;i ++)
{
if(s[i] == 'P')
{
now ++;
while(ask[j].y == now)
{
int pos = num[ask[j].x];
ask[j].ans = sum(suf[pos]) - sum(pre[pos] - 1);
j ++;
}
}
else if(s[i] == 'B')
{
add(pre[p],-1);
p = fa[p];
}
else
{
//  printf("%d %d\n",p,pre[p]);
p = ch[p][s[i] - 'a' + 1];
add(pre[p],1);
}
}
}

int main()
{

scanf("%s",s);
int n;
scanf("%d",&n);
for(int i = 1;i <= n;i ++)
{
scanf("%d%d",&ask[i].x,&ask[i].y);
ask[i].id = i;
}
sort(ask + 1,ask + 1 + n,cmp1);

insert(s);
build_ac();
for(int i = 1;i <= sz;i ++)
build(fail[i],i);
dfs(0);
//  for(int i = 1;i <= dfs_clock;i ++)
//      printf("%d ",dfs_num[i]); puts("");
solve();

sort(ask + 1,ask + 1 + n,cmp2);

for(int i = 1;i <= n;i ++)
printf("%d\n",ask[i].ans);

return 0;
}

/*
aPaPBbP
3
1 2
1 3
2 3

*/