【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 */