HDU 1403

题目描述：给出两个字符串，求最长公共子串的长度。

题解：后缀数组。将两个字符串合并成一个，中间加一个字符'0'，然后求出height数组。因为以最长公共子串为前缀的子串，他们的rank肯定相差为一，然后遍历sa数组，判断sa[i]和sa[i-1]的位置是否位于字符'0'的两侧，符合条件的最大的height值即为答案。

代码：

#include<bits/stdc++.h>
#define N 200010
using namespace std;

char s[N];
int sa[N],t[N],t2[N],c[N],n,rak[N],height[N];

void build_sa(int m,char *s)
{
    int i,*x=t,*y=t2;
    for (i=0;i<m;i++)c[i]=0;
    for (i=0;i<n;i++)c[x[i]=s[i]]++;
    for (i=1;i<m;i++)c[i]+=c[i-1];
    for (i=n-1;i>=0;i--)sa[--c[x[i]]]=i;
    for (int k=1;k<=n;k<<=1)
    {
        int p=0;
        for (i=n-k;i<n;i++)y[p++]=i;
        for (i=0;i<n;i++)if (sa[i]>=k)y[p++]=sa[i]-k;
        for (i=0;i<m;i++)c[i]=0;
        for (i=0;i<n;i++)c[x[y[i]]]++;
        for (i=0;i<m;i++)c[i]+=c[i-1];
        for (i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
        swap(x,y);
        p=1; x[sa[0]]=0;
        for (i=1;i<n;i++)
            x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
        if (p>=n) break;
        m=p;
    }
}

void getheight()
{
    int i,j,k=0;
    for (i=0;i<n;i++)rak[sa[i]]=i;
    for (i=0;i<n;i++)
    {
        if (k)k--;
        if (!rak[i])continue;
        j=sa[rak[i]-1];
        while (s[i+k]==s[j+k])k++;
        height[rak[i]]=k;
    }
}

int main()
{
    while (~scanf("%s",s))
    {
        n=strlen(s);
        s[n]='0';   int nn=n,ans=0;
        scanf("%s",s+n+1);
        n=strlen(s);
        build_sa(200,s);
        getheight();
        for (int i=1;i<n;i++)
        {
            int x=sa[i-1],y=sa[i];
            if (x>y) swap(x,y);
            if (x<nn && nn<y) ans=max(ans,height[i]);
        }
        printf("%d\n",ans);
    }
    return 0;
}