bzoj 3796 Mushroom追妹纸 后缀数组 kmp

一开始写了个后缀自动机，还调了半天，然后发现gg了。。。
果然后缀自动机不支持删除所有包含一个子串的节点（好像什么都不支持。。。。）

先把s1,s2连到一起，中间用非法字符隔开。
先kmp求有哪些位置有s3，然后可以求每个点开始往后不包含s3的最长前缀。
求一个后缀数组。从前往后扫一遍后缀数组，维护一下s2到当前位置的最长不包含s3的前缀，如果遇到一个s1的后缀那么更新一下答案。
再从后往前扫一遍后缀数组，更新一遍答案。

#include <bits/stdc++.h>
using namespace std;
#define A 128
#define N 110000
int n,m,ans;
int nex[N],vis[N],f[N],bel[N];
char s[N],s1[N];
int sa[N],h[N],rank[N],tr[N],has[N];
void kmp()
{
    for(int i=1,j=0;i<=m;)
    {
        nex[i]=j;
        while(j&&s1[j]!=s1[i])j=nex[j];
        i++;j++;
    }
    for(int i=1,j=1;i<=n;)
    {
        while(j&&s1[j]!=s[i])j=nex[j];
        if(j==m)
        {
            vis[i-m+1]=1,j=nex[j];
            while(j&&s1[j]!=s[i])j=nex[j];
        }
        i++;j++;
    }
}
int cmp(int x,int y,int k)
{
    if(x+k>n||y+k>n)return 0;
    return rank[x]==rank[y]&&rank[x+k]==rank[y+k];
}
void getsa()
{
    int i,cnt;
    for(i=1;i<=n;i++)has[s[i]]++;
    for(i=1,cnt=0;i<=A;i++)if(has[i])tr[i]=++cnt;
    for(i=1;i<=A;i++)has[i]+=has[i-1];
    for(i=1;i<=n;i++)rank[i]=tr[s[i]],sa[has[s[i]]--]=i;
    for(int k=1;cnt!=n;k<<=1)
    {
        for(i=1;i<=n;i++)has[i]=0;
        for(i=1;i<=n;i++)has[rank[i]]++;
        for(i=1;i<=n;i++)has[i]+=has[i-1];
        for(i=n;i>=1;i--)if(sa[i]>k)tr[sa[i]-k]=has[rank[sa[i]-k]]--;
        for(i=1;i<=k;i++)tr[n-i+1]=has[rank[n-i+1]]--;
        for(i=1;i<=n;i++)sa[tr[i]]=i;
        for(i=1,cnt=0;i<=n;i++)tr[sa[i]]=cmp(sa[i],sa[i-1],k) ? cnt:++cnt;
        for(i=1;i<=n;i++)rank[i]=tr[i];
    }
    for(i=1;i<=n;i++)
    {
        if(rank[i]==1)continue;
        for(int j=max(h[rank[i-1]]-1,1);;j++)
        {
            if(s[i+j-1]==s[sa[rank[i]-1]+j-1])h[rank[i]]=j;
            else break;
        }
    }
}
int main()
{
    //freopen("tt.in","r",stdin);
    scanf("%s",s+1);n=strlen(s+1);s[++n]='#';
    for(int i=1;i<n;i++)bel[i]=1;
    scanf("%s",s+n+1);n+=strlen(s+n+1);
    scanf("%s",s1+1);m=strlen(s1+1);
    kmp();
    for(int i=n;i>=1;i--)
    {
        if(vis[i])f[i]=0;
        else if(vis[i+1])f[i]=m;
        else f[i]=f[i+1]+1;
    }
    getsa();
    for(int i=1,tmp=0;i<=n;i++)
    {
        tmp=min(tmp,h[i]);
        if(bel[sa[i]])ans=max(ans,min(tmp,f[sa[i]]));
        else tmp=max(tmp,f[sa[i]]);
    }
    for(int i=n,tmp=0;i>=1;i--)
    {
        if(bel[sa[i]])ans=max(ans,min(tmp,f[sa[i]]));
        else tmp=max(tmp,f[sa[i]]);
        tmp=min(tmp,h[i]);
    }
    printf("%d\n",ans);
    return 0;
}