SPOJ LCS2 1812

题意：多个字符串的最长公共子串

解法：SAM

• 将第一个串建立自动机
• 后面的每个字符串就行匹配，维护每个位置的值，然后所有串中此位置取最小值
• 拓扑序，对父节点进行更新

/************************************************************************* 
    > File Name: lcs2.cpp 
    > Author: cy 
    > Mail: [e-mail]  
    > Created Time: 14/11/25 20:48:22 
 ************************************************************************/ 
 
#include<iostream> 
#include<cstring> 
#include <algorithm> 
#include<cstdlib> 
#include<vector> 
#include<cmath> 
#include<stdlib.h> 
#include<iomanip> 
#include<list> 
#include<deque> 
#include<map> 
#include <stdio.h> 
#include <queue> 
 
 
#define inf 0x3f3f3f3f 
  #define INF 0x3FFFFFFFFFFFFFFFLL 
#define rep(i,n) for(i=0;i<n;i++) 
 #define reP(i,n) for(i=1;i<=n;i++) 
 
#define ull unsigned long long 
 #define ll long long 
 
#define cle(a) memset(a,0,sizeof(a)) 
 
#define MAXN 500010+5 
#define MAXC 26 
using namespace std; 
char ch[MAXN]; 
int rk[MAXN],sz[MAXN],ans[MAXN]; 
 
struct Suffix_Automaton{ 
    int child[MAXN][MAXC]; 
 
    int f[MAXN];//指向最近的接受态 
    int l[MAXN];//此节点代表的最长串的长度 
    int head,tail,tot;//头节点,最后的接受态,总节点 
    void clear(){ 
        memset(child,0,sizeof(child)); 
        memset(f,0,sizeof(f)); 
        memset(l,0,sizeof(l)); 
        head=tail=tot=1; 
    } 
    void add(int c){ 
        int p=tail,np=++tot;l[np]=l[p]+1;tail=np; 
        for(;p&&!child[p][c];p=f[p]) child[p][c]=np; 
        if(!p) f[np]=head; 
        else if(l[child[p][c]]==l[p]+1) f[np]=child[p][c]; 
        else{ 
            int q=child[p][c],r=++tot; 
            memcpy(child[r],child[q],MAXC*sizeof(int)); 
            f[r]=f[q];l[r]=l[p]+1;f[q]=f[np]=r; 
            for(;p&&child[p][c]==q;p=f[p]) child[p][c]=r; 
        } 
    } 
    void insert(char*ch){ 
        int l=strlen(ch); 
        for(int i=0;i<l;i++) add(ch[i]-'a'); 
    } 
    int ws[MAXN]; 
    void tpu(){ 
        int len=strlen(ch); 
        for(int i=0;i<=len;i++)ws[i]=0; 
        for(int i=1;i<=tot;i++)ws[l[i]]++; 
        for(int i=1;i<=len;i++)ws[i]+=ws[i-1]; 
        for(int i=tot;i>0;i--)rk[ws[l[i]]--]=i; 
    } 
    int ans[MAXN]; 
    //为什么要沿着父节点更新? 
    int now[MAXN]; 
    void solve(){ 
        for(int i=1;i<=tot;i++)//最长可能性 
        { 
            now[i]=l[i]; 
        } 
        while(scanf("%s",&ch)!=EOF){ 
            cle(ans); 
            int len=strlen(ch); 
            int p=head; 
            int temp=0; 
            for(int i=0;i<len;i++) 
            { 
                int tt=ch[i]-'a'; 
                if(child[p][tt])temp++,p=child[p][tt]; 
                else{ 
                    while(true){ 
                        if(!p)break; 
                        if(child[p][tt])break; 
                        p=f[p]; 
                    } 
                    if(!p){ 
                        temp=0; 
                        p=head; 
                    } 
                    else{ 
                        temp=l[p]+1; 
                        p=child[p][tt]; 
                    } 
                } 
                ans[p]=max(temp,ans[p]); 
            } 
            for(int i=tot;i>0;i--) 
            { 
                ans[f[rk[i]]]=max(ans[f[rk[i]]],ans[rk[i]]); 
            } 
            for(int i=1;i<=tot;i++) 
            { 
                now[i]=min(now[i],ans[i]); 
            } 
        } 
        int ret=0; 
        for(int i=1;i<=tot;i++){ 
            ret=max(now[i],ret); 
        } 
        printf("%d\\n",ret); 
    } 
}SAM; 
int main() 
{ 
#ifndef ONLINE_JUDGE 
     freopen("in.txt","r",stdin); 
     //freopen("out.txt","w",stdout); 
#endif 
    scanf("%s",&ch); 
    SAM.clear(); 
    SAM.insert(ch); 
    SAM.tpu(); 
    SAM.solve(); 
    return 0; 
}