[SCOI2012]喵星球上的点名

对于第一个询问，对于姓名串建广义后缀自动机，那么我们对于插入一只喵的姓名串，要将这两个（姓和名）串的所有字串位置打标记，怎么搞呢，其实可以直接插入暴力向上跳，复杂度的话：对于一个串假设长度为x，那么插入复杂度是min（x*x,x*n），而最多只有n/x个这样的串，算一下总的复杂度上界是n根号的。然后就是第二个询问，刚开始naive了，以为再对询问串建广义后缀自动机暴力跑姓名串，之后发现如果失配要保留最大前缀没法搞只能用ac自动机。去看了一下广义sam的题解，发现直接把询问串也建在姓名串的广义sam里就是对的了，因为这样的话暴力跑姓名串就不会出现失配的情况了。。

代码：

#include<bits/stdc++.h>
using namespace std;
const int N=2e5+100;

int n,m;
vector<int>fi[N],sc[N],q[N];

struct Sam{
    int las,tot,fa[N],sz[N],szz[N],mx[N],vis[N],res;
    map<int,int>nxt[N];
    Sam(){las=tot=1;}
    void ins(int x)
    {
        int p=las,np=nxt[las][x],t,nt;
        if(np)
        {
            if(mx[np]==mx[p]+1)las=np;
            else
            {
                nt=++tot,mx[nt]=mx[p]+1;
                fa[nt]=fa[np],fa[np]=nt;
                nxt[nt]=nxt[np],las=nt,sz[nt]=sz[np];
                for(;p&&nxt[p][x]==np;p=fa[p])nxt[p][x]=nt;
            }
            return;
        }
        np=++tot,mx[np]=mx[p]+1,las=np;
        for(;p&&!nxt[p][x];p=fa[p])nxt[p][x]=np;
        if(!p){fa[np]=1;return;}
        t=nxt[p][x];
        if(mx[t]==mx[p]+1)fa[np]=t;
        else
        {
            nt=++tot,mx[nt]=mx[p]+1,nxt[nt]=nxt[t];
            fa[nt]=fa[t],fa[t]=fa[np]=nt,sz[nt]=sz[t];
            for(;p&&nxt[p][x]==t;p=fa[p])nxt[p][x]=nt;
        }
    }
    void upd1(int x,int y)
    {
        if(vis[x]==y)return;
        vis[x]=y,++sz[x];
        if(fa[x]!=1)upd1(fa[x],y);
    }
    void upd2(int x,int y)
    {
        if(vis[x]==y)return;
        vis[x]=y,res+=szz[x];
        if(fa[x]!=1)upd2(fa[x],y);
    }
    void in1()
    {
        int l1,l2;
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&l1),las=1;
            fi[i].resize(l1+1);
            for(int j=1;j<=l1;j++)
                scanf("%d",&fi[i][j]),++fi[i][j],ins(fi[i][j]);
            scanf("%d",&l2),las=1;
            sc[i].resize(l2+1);
            for(int j=1;j<=l2;j++)
                scanf("%d",&sc[i][j]),++sc[i][j],ins(sc[i][j]);
            for(int j=1,nw=1;j<=l1;j++)upd1(nw=nxt[nw][fi[i][j]],i);
            for(int j=1,nw=1;j<=l2;j++)upd1(nw=nxt[nw][sc[i][j]],i);
        }
    }
    void in2()
    {
        int l;
        for(int i=1,nw;i<=m;i++)
        {
            scanf("%d",&l),nw=las=1;
            q[i].resize(l+1);
            for(int j=1;j<=l;j++)
                scanf("%d",&q[i][j]),++q[i][j];
        }
    }
    void sol1()
    {
        int l;bool can;
        for(int i=1,nw;i<=m;i++)
        {
            l=q[i].size()-1,nw=can=1;
            for(int j=1;j<=l;j++)
            {
                if(!nxt[nw][q[i][j]])
                {
                    can=0;
                    break;
                }
                nw=nxt[nw][q[i][j]];
            }
            if(!can)puts("0");
            else printf("%d\n",sz[nw]),szz[nw]++;
        }
    }
    void sol2()
    {
        int l;
        memset(vis,0,sizeof vis);
        for(int i=1,nw;i<=n;i++)
        {
            res=0;
            l=fi[i].size()-1,nw=1;
            for(int j=1;j<=l;j++)
                nw=nxt[nw][fi[i][j]],upd2(nw,i);
            l=sc[i].size()-1,nw=1;
            for(int j=1;j<=l;j++)
                nw=nxt[nw][sc[i][j]],upd2(nw,i);
            printf("%d ",res);
        }
        puts("");
    }
}sam;

int main()
{
    scanf("%d%d",&n,&m);
    sam.in1(),sam.in2(),sam.sol1(),sam.sol2();
}