【SPOJ

题意

  求出多个串的最长公共子串。

分析

  　刚学SAM想做这个题的话最好先去做一下那道codevs3160。求两个串的LCS应该怎么求？把一个串s1建自动机，然后跑另一个串s2，然后找出s2每个前缀的最长公共后缀。那么多个的时候，我们也用这种类似的方法，但是我们求最长公共后缀的时候要求第一个串的。我们把其中一个串建SAM，然后把其他的串都在上面跑，维护两个值，Max[u]和Min[u]。自动机中每个状态u的Right存的是结尾集合。那么对于一个字符串，我们可以求出他和自动机中每个状态的最长公共后缀。然后，我们通过Max[fa[u]]=max(Max[fa[u]],Max[u])来确定左右状态的最长公共后缀，然后更新Min[o]。

下面的代码没有AC，但是我找了好久BUG没找到··如果有人看完并且找到了bug麻烦跟我说一下万分感谢！

　　

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

using namespace std;
const int maxn=300000+100;
char s[maxn];
struct state{
    int len,link;
    int next[26];
}st[2*maxn];
int n,last,cur,sz;
int Min[2*maxn],Max[2*maxn],c[2*maxn];

void init(){
    sz=1;
    cur=last=0;
    st[0].link=-1;
    st[0].len=0;
}

void build_sam(int c){
    cur=sz++;
    st[cur].len=st[last].len+1;
    Min[cur]=st[cur].len;
    int p;
    for(p=last;p!=-1&&st[p].next[c]==0;p=st[p].link)
        st[p].next[c]=cur;
    if(p==-1)
        st[cur].link=0;
    else{
        int q=st[p].next[c];
        if(st[q].len==st[p].len+1)
            st[cur].link=q;
        else{
            int clone=sz++;
            st[clone].len=st[p].len+1;
            Min[clone]=st[clone].len;
            st[clone].link=st[q].link;
            for(int i=0;i<26;i++)
                st[clone].next[i]=st[q].next[i];
            for(;p!=-1&&st[p].next[c]==q;p=st[p].link)
                st[p].next[c]=clone;
            st[cur].link=st[q].link=clone;
        }
    }
    last=cur;
}

int cmp(int a,int b){
    return st[a].len<st[b].len;
}

int main(){
    scanf("%s",s);
    n=strlen(s);
    init();

    for(int i=0;i<n;i++)
        build_sam(s[i]-'a');
    for(int i=0;i<sz;i++)
        c[i]=i;
    sort(c,c+sz,cmp);

    while(scanf("%s",s)!=EOF){
        n=strlen(s);
        int u=0,len=0;
        for(int i=0;i<n;i++){
            int c=s[i]-'a';
            if(u!=-1&&st[u].next[c]==0)
                u=st[u].link,len=st[u].len;
            if(u==-1)
                u=0,len=0;
            else{
                u=st[u].next[c];
                len++;
                Max[u]=max(Max[u],len);
            }
        }

        for(int i=sz-1;i>=0;i--){
            int o=c[i];
            Min[o]=min(Min[o],Max[o]);
            if(st[o].link!=-1){
                Max[st[o].link]=max(Max[st[o].link],Max[o]);
            }
            Max[o]=0;
        }
    }
    int ans=0;
    for(int i=0;i<sz;i++){
        ans=max(ans,Min[i]);
    }
    printf("%d
",ans);

return 0;
}