最长公共子串问题

求最长公共子串，基本就是求两个串或多个串~
方法也就是把一个串建SAM，剩下的进行匹配就好啦

1.spoj1811 LCS

题意：求两个串的最长公共子串

方法：用第一个串建立SAM，第二个串在第一个串上找就可以。有点类似ac自动机跑fail指针吧，匹配失败了就往前跳

不过好像不能用c++交诶..就很气

#include <stdio.h>
#include <string.h>
#define N 250010
int len=0,ans=0,n=0,last,root,son[N<<1][26],fa[N<<1],mx[N<<1];  
char str1[N],str2[N];
void ins(int ch){
    int p=last,np=++n;last=np;mx[np]=mx[p]+1;
    while(p && !son[p][ch]) son[p][ch]=np,p=fa[p];
    if(!p) fa[np]=root;
    else{
        int q=son[p][ch];
        if(mx[q]==mx[p]+1) fa[np]=q;
        else{
            int nq=++n;mx[nq]=mx[p]+1;
            memcpy(son[nq],son[q],sizeof(son[q]));
            fa[nq]=fa[q];fa[q]=fa[np]=nq;
            while(son[p][ch]==q) son[p][ch]=nq,p=fa[p];
        }
    }
} 
int main(){
    scanf("%s%s",str1+1,str2+1);
    int len1=strlen(str1+1),len2=strlen(str2+1);root=last=++n;
    for(int i=1;i<=len1;i++) ins(str1[i]-'a');
    int p=root;
    for(int i=1;i<=len2;i++){
        if(son[p][str2[i]-'a']) len++,p=son[p][str2[i]-'a'];
        else{
            while(p && !son[p][str2[i]-'a']) p=fa[p];
            if(!p) len=0,p=root;
            else len=mx[p]+1,p=son[p][str2[i]-'a'];
        }if(len>ans) ans=len;
    }printf("%d\n",ans);
    return 0;
} 

2.spoj1812 LCS2

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

方法：先根据第一个字符串建立SAM，然后后面的串在SAM上匹配。匹配时记录每个串在每个点能匹配最大值mx，然后根据拓扑序传给前面的点。再记录一个最小值表示到这个点被所有串匹配的值（因为要求所有串的lcs啊，不能只有一个能匹配）。结果就是这个最小值中的最大值

#include <stdio.h>
#include <string.h>
#define N 100010
#define inf 0x7fffffff
char str[N];
int n,last,root,cnt=0,ans=0,len[N<<1],son[N<<1][26],fa[N<<1],mx[N<<1],mn[N<<1],a[N<<1],c[N<<1];
inline int max(int x,int y){
    if(x>y) return x;
    return y;
}
inline int min(int x,int y){
    if(x<y) return x;
    return y;
}
inline void ins(int ch){
    int p=last,np=++cnt;last=np;len[np]=len[p]+1;
    while(p && !son[p][ch]) son[p][ch]=np,p=fa[p];
    if(!p) fa[np]=root;
    else{
        int q=son[p][ch];
        if(len[q]==len[p]+1) fa[np]=q;
        else{
            int nq=++cnt;len[nq]=len[p]+1;
            memcpy(son[nq],son[q],sizeof(son[nq]));
            fa[nq]=fa[q];fa[q]=fa[np]=nq;
            while(son[p][ch]==q) son[p][ch]=nq,p=fa[p];
        }
    }
}
int main(){
    scanf("%s",str+1);last=root=++cnt;n=strlen(str+1);
    for(int i=1;i<=n;i++) ins(str[i]-'a');
    for(int i=1;i<=cnt;i++) c[len[i]]++;
    for(int i=1;i<=cnt;i++) c[i]+=c[i-1];
    for(int i=cnt;i>=1;i--) a[c[len[i]]--]=i;
    for(int i=1;i<=cnt;i++) mn[i]=inf;
    memset(mx,0,sizeof(mx));
    while(scanf("%s",str+1)>0){
        n=strlen(str+1);int p=root,sz=0;
        for(int i=1;i<=n;i++){
            if(son[p][str[i]-'a']) sz++,p=son[p][str[i]-'a'];
            else{
                while(p&&!son[p][str[i]-'a']) p=fa[p];
                if(!p) sz=0,p=root;
                else sz=len[p]+1,p=son[p][str[i]-'a'];
            }mx[p]=max(mx[p],sz);
        }for(int i=1;i<=cnt;i++){
            int p=a[i];
            mn[p]=min(mn[p],mx[p]);
            mx[fa[p]]=max(mx[fa[p]],mx[p]);
            mx[p]=0;
        }
    }for(int i=1;i<=cnt;i++) if(mn[i]!=inf) ans=max(mn[i],ans);
    printf("%d\n",ans);
    return 0;
}