并不对劲的spoj1812

题意是求多个串的lcs。

这也是道后缀自动机的模板题。对于任意一个字符串建后缀自动机，用其他串查询就行。对于后缀自动机的每个状态要额外记匹配到当前状态的最大长度。

和spoj1811的区别在于这道题不方便后缀数组做。当然，如果不嫌很多个串用奇怪的字符连起来麻烦、判断时常数极大的话，也可以试试。 字符串总数不多，好像还真行。

#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define maxn 100010
using namespace std;
int ans,len,p;
int read(){
    int f=,x=;char ch=getchar();
    while(isdigit(ch)== && ch!='-')ch=getchar();
    if(ch=='-')f=-,ch=getchar();
    while(isdigit(ch))x=x*+ch-'',ch=getchar();
    return x*f;
}
void write(int x){
    int ff=;char ch[];
    while(x)ch[++ff]=(x%)+'',x/=;
    if(ff==)putchar('');
    while(ff)putchar(ch[ff--]);
    putchar(' ');
}
typedef struct node{
    int to[],dis,fa;
}spot;
struct SAM{
    spot x[maxn*];
    int ls,ls2,q,cnt,rt,lst,c[maxn*],ord[maxn*],ans[maxn*],maxl[maxn*];
    char s[maxn];
    void start(){
        lst=rt=++cnt;
        scanf("%s",s+);
        ls=strlen(s+);
        for(int i=;i<=ls;i++)
            extend(i);
        for(int i=;i<=cnt;i++)
            ans[i]=x[i].dis;
    }
    void extend(int pos){
        int val=s[pos]-'a',p=lst,np=++cnt;
        lst=np,x[np].dis=pos;
        for(;p&&x[p].to[val]==;p=x[p].fa)x[p].to[val]=np;
        if(p==)x[np].fa=rt;
        else{
            int q=x[p].to[val];
            if(x[q].dis==x[p].dis+)x[np].fa=q;
            else{
                int nq=++cnt;
                x[nq].dis=x[p].dis+;
                memcpy(x[nq].to,x[q].to,sizeof(x[q].to));
                x[nq].fa=x[q].fa,x[np].fa=x[q].fa=nq;
                for(;x[p].to[val]==q;p=x[p].fa)x[p].to[val]=nq;
            }
        }
    }
    void qsort(){
        for(int i=;i<=cnt;i++)
            c[x[i].dis]++;
        for(int i=;i<=ls;i++)
            c[i]+=c[i-];
        for(int i=;i<=cnt;i++)
            ord[c[x[i].dis]--]=i;
    }
    int work(){
        while(scanf("%s",s+)!=EOF){
            len=;
            memset(maxl,,sizeof(maxl));
            ls2=strlen(s+);
            for(int i=;i<=ls2;i++){
                int val=s[i]-'a';
                if(x[p].to[val])len++,p=x[p].to[val];
                else{
                    for(;p&&x[p].to[val]==;p=x[p].fa);
                    if(p==)p=rt,len=;
                    else len=x[p].dis+,p=x[p].to[val];
                }
                maxl[p]=max(maxl[p],len);
            }
            for(int i=cnt;i>=;i--){
                int u=ord[i];
                ans[u]=min(ans[u],maxl[u]);
                if(x[u].fa && maxl[u])maxl[x[u].fa]=x[x[u].fa].dis;
            }
        }
        int res=;
        for(int i=;i<=cnt;i++){
            res=max(res,ans[i]);
        }
        return res;
    }
}t;
int main(){
    t.start();
    t.qsort();
    int res=t.work();
    write(res);
    return ;
}

并不对劲的lcs