hdu4416[多串后缀自动机]

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4416
题意：给一个A串和多个B串，求A串中有多少个子串不是任何一个B串的子串。
分析：统计子串的个数，最直接的就想到后缀自动机。首先把第一个串建成后缀自动机，当把其他的串加入的时候，我们首先让last指针指向根节点，然后与单串类似，逐个加入字符，不同的是现在的当前后缀可能是原本的某个后缀，即$ch[last][c]!=null$，这时候需要比较$ml[ch[last][c]]$是否和$ml[p]+1$是否相等（与构造算法类似，不明白的话可以把abab ba 画出来就懂了）

#include<bits/stdc++.h>
using namespace std;
const int Maxn=4023456;
typedef long long Int;
char s[100020];
int n,cs=1;

struct suffix_automaton
{
    int ml[Maxn],f[Maxn],ch[Maxn][26];
    int sz,last,maxlen;
    int col[Maxn];
    int id[Maxn],cnt[100020];
    void init()
    {
        memset(ch[1],0,sizeof(ch[1]));
        sz=2;
        last=1;
    }
    int getidx(char c){return c-'a';}
    void add(char x,int w)
    {
        int c=getidx(x);
        int p=last;
        if(ch[p][c])
        {
            if(ml[ch[p][c]]==ml[p]+1){last=ch[p][c];col[last]=w;}
            else
            {
                int nq=sz++,q=ch[p][c];
                last=nq;
                col[nq]=w;
                memcpy(ch[nq],ch[q],sizeof(ch[q]));
                ml[nq]=ml[p]+1;
                f[nq]=f[q];
                f[q]=nq;
                for(;p&&ch[p][c]==q;p=f[p])ch[p][c]=nq;
            }
            return;
        }
        int np=sz++;ml[np]=ml[p]+1;last=np;
        col[np]=w;
        memset(ch[np],0,sizeof(ch[np]));
        for(;p&&!ch[p][c];p=f[p])ch[p][c]=np;
        if(!p){f[np]=1;return;}
        int q=ch[p][c];
        if(ml[q]==ml[p]+1){f[np]=q;return;}
        int nq=sz++;
        col[nq]=w;
        memcpy(ch[nq],ch[q],sizeof(ch[q]));
        ml[nq]=ml[p]+1;
        f[nq]=f[q];
        f[q]=nq;
        for(;p&&ch[p][c]==q;p=f[p])ch[p][c]=nq;
        f[np]=nq;
    }
    void solve()
    {
        scanf("%s",s);
        int i,j;
        for(i=0;s[i];i++)add(s[i],1);maxlen=i;
        for(i=0;i<n;i++)
        {
            last=1;
            scanf("%s",s);
            for(j=0;s[j];j++)add(s[j],0);maxlen=max(maxlen,j);
        }
        for(i=1;i<=maxlen;i++)cnt[i]=0;
        for(i=1;i<sz;i++)cnt[ml[i]]++;
        for(i=1;i<=maxlen;i++)cnt[i]+=cnt[i-1];
        for(i=1;i<sz;i++)id[cnt[ml[i]]--]=i;
        for(i=sz-1;i>=1;i--)col[f[id[i]]]&=col[id[i]];
        Int ans=0;
        for(i=1;i<sz;i++)
            if(col[i])
                ans+=ml[i]-ml[f[i]];
        printf("Case %d: %lld\n",cs++,ans);
    }
}solver;
int main()
{
    int _;scanf("%d",&_);
    while(_--)
    {
        scanf("%d",&n);
        solver.init();
        solver.solve();
    }
}