[CODEVS3160]最长公共子串|后缀数组|后缀自动机

  第一种做法是把两个字符串接起来，中间放一个奇怪的字符，然后建这个串的后缀数组，求出h数组，对于h[i]表示lcp（sa[i],sa[i-1]），如果sa[i]和sa[i-1]分布在奇怪的字符的两边就用h[i]更新答案。。
  然后补了一下后缀自动机，就是裸题了，对一个串建SAM，另一个串在上面匹配即可。。
SA:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<memory.h>
#define N 400005
#define clr(a) memset(a,0,sizeof(a))
using namespace std;
char s[N],b[N];
int n,i,d,sa[N],rank[N],h[N],buc[N],sec[N],rk[N];
bool cmp(int a,int b)
{
    return rk[a]==rk[b]&&rk[a+d]==rk[b+d];
}
void getsa()
{
    int i,j,p=1;
    d=1;clr(buc);clr(rank);
    for (i=1;i<=n;i++) buc[rank[i]=s[i]-'a'+1]++;
    for (i=1;i<=256;i++) buc[i]+=buc[i-1];
    for (i=n;i;i--) sa[buc[rank[i]]--]=i;
    while (p<n)
    {
        p=0;clr(buc);
        for (i=n-d+1;i<=n;i++) sec[++p]=i;
        for (i=1;i<=n;i++) if (sa[i]>d) sec[++p]=sa[i]-d;
        for (i=1;i<=n;i++) buc[rank[sec[i]]]++,rk[i]=rank[i];
        for (i=1;i<=n;i++) buc[i]+=buc[i-1];
        for (i=n;i;i--) sa[buc[rank[sec[i]]]--]=sec[i];
        for (i=2,p=1,rank[sa[1]]=1;i<=n;i++)
            rank[sa[i]]=cmp(sa[i-1],sa[i])? p:++p;
        d<<=1;
    }
}
void geth()
{
    int i,j,k=0;
    h[1]=0;
    for (i=1;i<=n;i++)
        if (rank[i]>1)
        {
            j=sa[rank[i]-1];
            while (s[i+k]==s[j+k]) k++;
            h[rank[i]]=k;if (k) k--;
        }
}
int main()
{
    scanf("%s%s",s+1,b+1);
    n=strlen(s+1);
    int n1=n+1;
    s[n+1]='z'+1;
    for (i=n+2;i<=n+1+strlen(b+1);i++) s[i]=b[i-n-1];
    n+=strlen(b+1)+1;
    getsa();geth();
    int ans=0;  
    for (i=2;i<=n;i++)
        if ((sa[i]<n1&&sa[i-1]>n1)||(sa[i]>n1&&sa[i-1]<n1)) ans=max(ans,h[i]);
    printf("%d\n",ans);
}

SAM:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<memory.h>
#define N 200005
using namespace std;
char s1[N],s2[N];
int n,ans,i,now,nd,a[N][26],fa[N],l[N];
void extend(int c)
{
    int p=now,np=now=++nd;
    l[np]=l[p]+1;
    while (!a[p][c]&&p) a[p][c]=np,p=fa[p];
    if (!p) fa[np]=1;
    else
    {
        int q=a[p][c];
        if (l[q]==l[p]+1) fa[np]=q;
        else
        {
            int nq=++nd;l[nq]=l[p]+1;
            fa[nq]=fa[q];
            fa[q]=fa[np]=nq;
            for (int i=0;i<26;i++) a[nq][i]=a[q][i];
            while (p&&a[p][c]==q) a[p][c]=nq,p=fa[p];
        }
    }
}
void build()
{
    for (int i=1;i<=strlen(s1+1);i++) extend(s1[i]-'a');
}
void solve()
{
    ans=0;now=1;
    int tmp=0;
    for (i=1;i<=strlen(s2+1);i++)
    {
        if (a[now][s2[i]-'a']) tmp++,now=a[now][s2[i]-'a'];
        else
        {
            while (!a[now][s2[i]-'a']) now=fa[now];
            if (!now) tmp=0,now=1;else tmp=l[now]+1,now=a[now][s2[i]-'a'];  
        }
        ans=max(ans,tmp);
    }
}
int main()
{
    scanf("%s",s1+1);
    memset(a,0,sizeof(a));
    l[1]=0;now=nd=1;fa[1]=0;
    for (i=0;i<26;i++) a[0][i]=1;
    build();
    scanf("%s",s2+1);
    solve();
    printf("%d\n",ans);
}