【WC2016模拟】打击目标

Description

Input

Output

Sample Input

3 0
xyz xy x
1 2
2
1 3 xyz
2 3 zzzzxy

Sample Output

3
2

Data Constraint

题解

这一题在比赛的时候我已经想到了和题解基本一样的做法，然而因为不会树剖（雾，所以只打了50point，结果数据超级水，要不是开小了数组就AC了QvQ
下午赶紧学习了一下树剖，然后发现是一个超级简单的东西，就是维护出每一个点对应的重链的上面的顶点是什么，然后看哪边深就跳哪边就好了

那么这一题怎么做呢？

收先看到多串匹配条件反射的就想到了SA和AC自动机
看了一下发现就是一个AC自动机
某一个串经过的自动机中某一节点沿着fail能跳到的就是它的子串

那么我们不妨对自动机中的每一个点挂一个可持久化线段树，然后沿着fail边不断的更新就好了

询问的话因为树剖了一下所以只有log段，复杂度应该是两个log的

贴代码

打的不是很优美。。。

#include<iostream>#include<algorithm>#include<cstdio>#include<cstring>#include<cmath>#define fo(i,a,b) for(i=a;i<=b;i++)#define fo1(i,b,a) for(i=b;i>=a;i--)#define max(x,y) ((x)>(y)?(x):(y))using namespace std;const int maxn=1e5+5;int tree[maxn*10][27],fail[maxn*10],h[maxn*10];int tr[maxn*50][3];int ff[maxn][18];int cc[maxn*10],de[maxn],dep[maxn*10];int fi[maxn],ne[maxn*2],dui[maxn*2],qc[maxn],siz[maxn];int fi1[maxn*3],ne1[maxn*5],dui1[maxn*5],qc1[maxn*3];int fa[maxn],dfn[maxn],top[maxn],go[maxn],root[maxn*10];bool bq;int i,j,k,l,m,n,x,y,z,o,ans,now,q,xz,zo,zz,cz,mx,xx,yy,p,xc;char s[maxn];void add(int x,int y){    if (fi[x]==0) fi[x]=++mx; else ne[qc[x]]=++mx;    dui[mx]=y; qc[x]=mx;}void add1(int x,int y){    if (fi1[x]==0) fi1[x]=++mx; else ne1[qc1[x]]=++mx;    dui1[mx]=y; qc1[x]=mx;}void ge_a(){    k=0;    fo(i,1,26) if (tree[0][i]) h[++k]=tree[0][i];    i=k; j=0;    while (i>j){        j++; x=h[j];        fo(k,1,26) if (tree[x][k]){            now=fail[x];            while (tree[now][k]==0 && now>0) now=fail[now];            if (tree[now][k]) now=tree[now][k];            fail[tree[x][k]]=now;            h[++i]=tree[x][k];        }    }}void ge_ff(){    fo(j,1,16)        fo(i,1,n) ff[i][j]=ff[ff[i][j-1]][j-1];}int lca(int x,int y){    if (de[x]<de[y]){        int z=x; x=y; y=z;    }    fo1(i,16,0) if (de[x]-(1<<i)>=de[y]) x=ff[x][i];    fo1(i,16,0) if (ff[x][i]!=ff[y][i]){        x=ff[x][i]; y=ff[y][i];    }    if (x!=y) x=ff[x][0];    return x;}void dfs(int x){    int i=fi[x]; siz[x]=1;     while (i){        de[dui[i]]=de[x]+1;        dfs(dui[i]);        siz[x]+=siz[dui[i]];        if (siz[dui[i]]>siz[go[x]]) go[x]=dui[i];        i=ne[i];    }}void dfss(int x){    dfn[x]=++now;    top[go[x]]=top[x];    if (go[x]) dfss(go[x]);    int i=fi[x];    while (i){        if (dui[i]==go[x]){            i=ne[i]; continue;        }        top[dui[i]]=dui[i];        dfss(dui[i]);        i=ne[i];    }}void maketree(int v,int s1,int s2,int l,int r){    tr[++p][0]=tr[v][0]; tr[p][1]=tr[v][1]; tr[p][2]=tr[v][2];    int q=p;    if (l==r){        tr[p][2]=s2; return;    }    int mid=(l+r)/2;    if (s1<=mid){        tr[p][0]=p+1;        maketree(tr[v][0],s1,s2,l,mid);    } else{        tr[p][1]=p+1;        maketree(tr[v][1],s1,s2,mid+1,r);    }    tr[q][2]=max(tr[tr[q][0]][2],tr[tr[q][1]][2]);}void wor(int x){    if (! root[fail[x]] && fail[x]) wor(fail[x]);    if (fi1[x]){        int pp=fi1[x];        root[x]=p+1;        maketree(root[fail[x]],dfn[dui1[pp]],dep[x],1,n); pp=ne1[pp];        while (pp){            int qq=p+1;            maketree(root[x],dfn[dui1[pp]],dep[x],1,n);            root[x]=qq;            pp=ne1[pp];        }    }        else root[x]=root[fail[x]];}void find(int v,int l,int r,int x,int y){    if (l==x && r==y) ans=max(ans,tr[v][2]); else{        int mid=(l+r)/2;        if (y<=mid) find(tr[v][0],l,mid,x,y); else        if (x>mid) find(tr[v][1],mid+1,r,x,y); else{            find(tr[v][0],l,mid,x,mid);            find(tr[v][1],mid+1,r,mid+1,y);        }    }}void tur(int q){    xx=x; yy=y;    while (top[xx]!=top[yy]){        if (de[top[xx]]<de[top[yy]]){            find(root[q],1,n,dfn[top[yy]],dfn[yy]);            yy=fa[top[yy]];        } else{            find(root[q],1,n,dfn[top[xx]],dfn[xx]);            xx=fa[top[xx]];        }    }    if (dfn[xx]>dfn[yy]){        z=xx; xx=yy; yy=z;    }    find(root[q],1,n,dfn[xx],dfn[yy]);}int main(){//  freopen("t2.in","r",stdin);//  freopen("t2.out","w",stdout);    scanf("%d%d",&n,&o);    fo(i,1,n){        scanf("%s",s+1); l=strlen(s+1);        now=0;        fo(j,1,l){            x=s[j]-96;            if (tree[now][x]) now=tree[now][x]; else{                tree[now][x]=zo+1; dep[++zo]=dep[now]+1;                 now=tree[now][x];            }        }        add1(now,i);    }    ge_a(); mx=0;    fo(i,2,n){        scanf("%d",&ff[i][0]); fa[i]=ff[i][0];        add(ff[i][0],i);    }    ge_ff();    now=0;    dfs(1); top[1]=1;     dfss(1); p=0;    fo(i,1,zo)        if (! root[i]) wor(i);    scanf("%d",&q); ans=0;    fo(xz,1,q){        scanf("%d%d",&x,&y);        if (o){            x^=ans; y^=ans;        } xx=x; yy=y;        z=lca(x,y);        now=0; scanf("%s",s+1); l=strlen(s+1); ans=0;        fo(i,1,l){            xc=s[i]-96;            while (now>0 && tree[now][xc]==0) now=fail[now];            if (tree[now][xc]) now=tree[now][xc];            tur(now);        }        printf("%d\n",ans);    }    return 0;}