hdu 4416 后缀自动机 问在S中有多少个不同子串满足它不是s1~sn中任意一个字符串的子串

题意：给一个字符串S和一系列字符串s1~sn，问在S中有多少个不同子串满足它不是s1~sn中任意一个字符串的子串

后缀自动机 对sa建立串，然后找每个串在母串上每个点的最大匹配长度，然后最后从后往前拓扑更新就好，因为叶子的max是能影响父亲的max，最后统计的时候 看Maxi是否有值，一旦有值，必定是这个点上的合法的值，那么就是统计上 Max[p]-Maxi[p] 而且往父节点更新，否则就是Max[p]-Max[fail[p]]

#include <bits/stdc++.h>
using namespace std;

const int maxn = 200000+5;
int last,tail,in[maxn],Min[maxn];
int Max[maxn],cnt[maxn],vis[maxn];
int nxt[maxn][26],fail[maxn];

char sa[maxn],sb[maxn];
int Maxi[maxn];


inline void build(char *s)
{
    while(*s)
    {
        int p=last,t=++tail,c=*s++-'a';
        Max[t]=Max[p]+1;
        Maxi[t]=0;
        while(p&&!nxt[p][c])
            nxt[p][c]=t,p=fail[p];
        if(p)
        {
            int q=nxt[p][c];
            if(Max[q]==Max[p]+1)
                fail[t]=q,Min[t]=Max[q]+1;
            else
            {
                int k=++tail;
                fail[k]=fail[q];
                fail[t]=fail[q]=k;
                Max[k]=Max[p]+1;
                Maxi[k]=0;
                memcpy(nxt[k],nxt[q],sizeof(nxt[q]));
                while(p&&nxt[p][c]==q)
                    nxt[p][c]=k,p=fail[p];
            }
        }
        else 
            fail[t]=Min[t]=1;
        last=t;
    }
}


int b[maxn];

int main()
{
    int t;
    scanf("%d",&t);
    int cas=1;
    while(t--)
    {   
        int n;
        memset(nxt,0,sizeof(nxt));
        memset(b,0,sizeof(b));
        memset(cnt,0,sizeof(cnt));
        scanf("%d",&n);
        scanf("%s",sa);
        last=1,tail=1;
        build(sa);
        int hh=strlen(sa);

        for(int i=1;i<=tail;i++) cnt[Max[i]]++;
        for(int i=1;i<=hh;i++) cnt[i]+=cnt[i-1];
        for(int i=1;i<=tail;i++) b[cnt[Max[i]]--]=i;

        for(int i=1;i<=n;i++)
        {
            scanf("%s",sb);
            int h=strlen(sb);
            int p=1;
            int len=0;
            for(int j=0;j<h;j++)
            {
                int c=sb[j]-'a';
                if(nxt[p][c])
                    len++,p=nxt[p][c];
                else{
                    while(p&&!nxt[p][c])
                        p=fail[p];
                    if(!p) p=1,len=0;
                    else
                    {
                        len=Max[p]+1;
                        p=nxt[p][c];
                    }
                }
                Maxi[p]=max(Maxi[p],len);
            }
        }
        long long ans=0;
        for(int i=tail;i>=1;i--)
        {
            int p=b[i];
            if(Maxi[p])
            {
                if(fail[p]) 
                    Maxi[fail[p]]=max(Maxi[fail[p]],Maxi[p]);
                if(Maxi[p]<Max[p]) ans+=Max[p]-Maxi[p];
            }
            else ans+=Max[p]-Max[fail[p]];
        }
        printf("Case %d: %lld\n",cas++,ans );
    }
}