SPOJ1812 LCS2 - Longest Common Substring II（SAM）

题目链接

分析：
SAM求多串最长公共子串
详解戳这里

#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

const int N=100002;
int dis[N<<1],ch[N<<1][26],fa[N<<1],len,root=1,last=1,sz=1;
int l[N<<1],c[N<<1],pos[N<<1],mn[N];
char s[N];

void insert(int x)
{
    int now=++sz,pre=last;
    last=now;
    dis[now]=dis[pre]+1;
    for (;pre&&!ch[pre][x];pre=fa[pre]) ch[pre][x]=now;

    if (!pre) fa[now]=root;
    else
    {
        int q=ch[pre][x];
        if (dis[q]==dis[pre]+1) fa[now]=q;
        else
        {
            int nows=++sz;
            dis[nows]=dis[pre]+1;
            memcpy(ch[nows],ch[q],sizeof(ch[nows]));
            fa[nows]=fa[q]; fa[q]=fa[now]=nows;
            for (;pre&&ch[pre][x]==q;pre=fa[pre]) ch[pre][x]=nows;
        }
    }
}

void solve()
{
    memset(l,0,sizeof(l));
    int now=root,tmp=0;
    for (int i=1;i<=len;i++)
    {
        int x=s[i]-'a';
        if (ch[now][x]) now=ch[now][x],tmp++;
        else
        {
            while (now&&!ch[now][x]) now=fa[now];
            if (!now) now=root,tmp=0;
            else tmp=dis[now]+1,now=ch[now][x];
        }
        l[now]=max(l[now],tmp);
    }

    for (int i=sz;i>=1;i--)
    {
        int t=pos[i];        //排名第i的结点 
        mn[t]=min(mn[t],l[t]);
        if (fa[t]&&l[t]) l[fa[t]]=dis[fa[t]];
    }
}

int main()
{
    freopen("a.in","r",stdin);
    scanf("%s",s+1);
    len=strlen(s+1);
    for (int i=1;i<=len;i++) insert(s[i]-'a');

    for (int i=1;i<=sz;i++) ++c[dis[i]];
    for (int i=1;i<=len;i++) c[i]+=c[i-1];
    for (int i=1;i<=sz;i++) pos[c[dis[i]]--]=i;   //dis从大到小 

    memset(mn,127,sizeof(mn));
    int ans=0;
    while (scanf("%s",s+1)!=EOF)
    {
        len=strlen(s+1);
        solve();
    }
    for (int i=1;i<=sz;i++) ans=max(ans,mn[i]);
    printf("%d",ans);
}