hdu 4416 Good Article Good sentence(后缀自动机)

题目链接：hdu 4416 Good Article Good sentence

题意：

给你一个串A和n个串B，问你A有多少个子串不是这n个B的子串。

题解：

将A串建立后缀自动机，对于每个B串都拿去匹配一下，并记录后缀自动机中每个节点的最大匹配长度。

然后拓扑排序，更新每个节点的fail节点。最后对于每个节点的贡献就是ml[i]-max(is[i],mx[f[i]]) (is[i]是该节点的最大匹配长度)

#include<bits/stdc++.h>
#define F(i,a,b) for(int i=a;i<=b;++i)
#define mst(a,b) memset(a,b,sizeof(a))
using namespace std;

const int N=1e5+7,tyn=26,M=N*2;
struct SAM{
    int tr[M][tyn],f[M],ml[M],ed,last,p,x,r,q;
    int cnt[M],b[M],d[M],is[M];
    inline int gid(char x){return x-'a';}
    inline void nc(int s,int &p){
        ml[p=++ed]=s,f[ed]=cnt[ed]=0,mst(tr[ed],0);
    }    
    void clear(){ed=0,nc(0,last);}
    void add(int w){
        nc(ml[x=last]+1,p),last=p,cnt[p]=1;
        while(x&&!tr[x][w])tr[x][w]=p,x=f[x];
        if(!x)f[p]=1;
        else if(ml[x]+1==ml[q=tr[x][w]])f[p]=q;
        else{
            nc(ml[x]+1,r),f[r]=f[q],f[p]=f[q]=r;
            memcpy(tr[r],tr[q],sizeof(tr[r]));
            while(x&&tr[x][w]==q)tr[x][w]=r,x=f[x];
        }
    }
    void build(char *s){for(int i=1;s[i];i++)add(gid(s[i]));}
    void upright(int mx=0){
        F(i,0,ed)d[i]=0;
        F(i,1,ed)d[mx<ml[i]?mx=ml[i]:ml[i]]++;
        F(i,1,mx)d[i]+=d[i-1];
        F(i,1,ed)b[d[ml[i]]--]=i;
        for(int i=ed;i;--i)if(f[b[i]])is[f[b[i]]]=max(is[f[b[i]]],is[b[i]]);
    }
    void go(char *s,int x=1)
    {
        int cl=0;//与s的匹配长度
        for(int i=1;s[i];i++)
        {
            int w=gid(s[i]);
            if(tr[x][w])x=tr[x][w],cl++;
            else
            {
                while(x&&!tr[x][w])x=f[x];
                if(!x)cl=0,x=1;else cl=ml[x]+1,x=tr[x][w];
            }
            is[x]=max(is[x],cl);
        }
    }
}sam;

char s[N];
int t,n,cas;
long long ans;

int main()
{
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%s",&n,s+1),ans=0;
        sam.clear(),sam.build(s);
        F(i,1,sam.ed)sam.is[i]=0;
        F(i,1,n)scanf("%s",s+1),sam.go(s);
        sam.upright();
        F(i,1,sam.ed)if(sam.is[i]<sam.ml[i])
        ans+=sam.ml[i]-max(sam.is[i],sam.f[i]!=0?sam.ml[sam.f[i]]:0);
        printf("Case %d: %lld
",++cas,ans);
    }
    return 0;
}