洛谷P2414 jzoj 2784.【NOI2011】阿狸的打字机(AC自动机)

传送门

Description:

阿狸喜欢收藏各种稀奇古怪的东西,最近他淘到一台老式的打字机。打字机上只有28个按键,分别印有 26个小写英文字母和’B’、’P’两个字母。 经阿狸研究发现,这个打字机是这样工作的:
输入小写字母,打字机的一个凹槽中会加入这个字母(按 P 前凹槽中至少有一个字母)。
按一下印有’B’的按键,打字机凹槽中最后一个字母会消失。
按一下印有’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


分析

不会AC自动机的可以看看我的blog
首先,我们构建一个AC自动机。
朴素做法是找到y在tire中的路径,沿着这条路径上每一个节点的fail指针走,找到一个x答案就加1。
最坏时间复杂度 O ( n m )
估计只能拿30分。
换一种思路,从x沿着fail指针倒着走,统计走到y上多少次。
这样我们就需要引进fail树。
fail树就是以AC自动机中的fail边构建的一棵树。
fail树有一个性质:对于一个节点x,它的子树中的所有节点所表示的字符串都包含x所表示的字符串(联想一下fail指针的构建过程)。
构建fail树后用dfs序标号,那么就可以用树状数组维护子树。
我们沿着tire走,走到一个点就把它按dfs序标号加入树状数组。
如何统计答案?假设我们有一个询问x,y,在AC自动机中我们走到了y,那么我们只需要查询x的子树就可以了。
回溯时将节点对应树状数组的值减去。


code:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#define low(x) (x&-x)
using namespace std;
const int N=1e5+10;
char S[N];
int n,len,fail[N],a[N][26],tot,fa[N],dfn[N],tr[N],b[N],l[N],r[N],tt,ans[N];
bool bz[N][26];
struct node{
    int top,ne[N],to[N],last[N],num[N];
    void add(int x,int y,int z) {
        to[++top]=y;ne[top]=last[x];last[x]=top;
        num[top]=z;
    }
}ac,e;
queue<int > d;

void change(int x,int y) {
    while (x<=tot) 
        tr[x]+=y,x+=low(x);
}
int sum(int x) {
    int y=0;
    while (x) 
        y+=tr[x],x-=low(x);
    return y;
}

int makenew() {
    memset(a[++tot],-1,sizeof(a[tot]));
    return tot;
}

void build() {
    int i;
    for (i=0;i<26;i++) {
        int &x=a[1][i];
        if (x==-1) x=1,bz[1][i]=true; else fail[x]=1,d.push(x);
    }
    while (!d.empty()) {
        int x=d.front();
        for (i=0;i<26;i++) {
            int y=a[fail[x]][i];
            if (a[x][i]==-1) a[x][i]=y,bz[x][i]=true; 
            else fail[a[x][i]]=y,d.push(a[x][i]);
        }
        d.pop();
    }
    for (i=1;i<=tot;i++) ac.add(fail[i],i,0);
}

void dfs(int x) {
    dfn[x]=l[x]=++tt;
    for (int i=ac.last[x];i;i=ac.ne[i]) 
        dfs(ac.to[i]);
    r[x]=tt;
}

void dfs1(int x) {
    change(dfn[x],1);
    for (int i=0;i<26;i++)
        if (!bz[x][i]) {
            dfs1(a[x][i]);
        }
    for (int i=e.last[x];i;i=e.ne[i]){
        int y=e.to[i],z=e.num[i];
        if (z==1)
            printf("");
        ans[z]=sum(r[y])-sum(l[y]-1);
    }
    change(dfn[x],-1);
}

int main() {
    int i,j,k;
    scanf("%s",S);
    len=strlen(S);
    int x=1;
    makenew();
    for (i=0;i<len;i++) {
        char c=S[i];
        if (c=='B') {
            x=fa[x];
        }
        else {
            if (c!='P') {
                if (a[x][c-'a']==-1) a[x][c-'a']=makenew(),fa[tot]=x;
                x=a[x][c-'a'];
            }
            else b[++b[0]]=x;
        }
    }
    build();
    dfs(1);
    scanf("%d",&n);
    for (i=1;i<=n;i++) {
        int x,y;scanf("%d%d",&x,&y);
        e.add(b[y],b[x],i);
    }
    dfs1(1);
    for (i=1;i<=n;i++) printf("%d\n",ans[i]);
}

猜你喜欢

转载自blog.csdn.net/kerGH/article/details/80882038