POJ2774

题意：给定两个字符串，求两个串的最长公共子路。两个串的长度不超过100000。

算法：把两个串用一个连接符连接起来，然后求height[]和sa[]，扫描height[]，如果sa[i]和sa[i-1]分别在两个串中的最大hgight[]的值。

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=200010;
char c[100005],cc[100005];
int s[maxn],sa[maxn],rank[maxn],height[maxn],saf[maxn],cs[maxn];
int n,m,mid;
void init()
{
    scanf("%s",c);scanf("%s",cc);
    int p=0;
    for(int i=0;c[i];++i)s[++p]=c[i];++p;
    mid=p;
    for(int i=0;cc[i];++i)s[++p]=cc[i];
    n=p;
    s[0]=s[n+1]=-1;
}
void rsort()
{
    for(int i=0;i<=m;++i)cs[i]=0;
    for(int i=1;i<=n;++i)cs[rank[saf[i]]]++;
    for(int i=1;i<=m;++i)cs[i]+=cs[i-1];
    for(int i=n;i>=1;--i)sa[cs[rank[saf[i]]]--]=saf[i];
}
int cmp(int *f,int x,int y,int w)
{
    return f[x]==f[y] && f[x+w]==f[y+w];
}
void suffix()
{
    for(int i=1;i<=n;++i)rank[i]=s[i],saf[i]=i;
    m=127;rsort();
    for(int p=1,w=1,i;p<n;w<<=1,m=p)
    {
        for(i=n-w+1,p=0;i<=n;++i)saf[++p]=i;
        for(i=1;i<=n;++i)if(sa[i]>w)saf[++p]=sa[i]-w;
        rsort();swap(rank,saf);rank[sa[1]]=p=1;
        for(i=2;i<=n;++i)rank[sa[i]]=cmp(saf,sa[i],sa[i-1],w)?p:++p;
    }
    int j,k=0;
    for(int i=1;i<=n;height[rank[i++]]=k)
        for(k=k?k-1:k,j=sa[rank[i]-1];s[i+k]==s[j+k];++k);
}
bool ind1(int x)
{
    return x>=1 &&x<mid;
}
bool ind2(int x)
{
    return x>mid && x<=n;
}
void work()
{
    int ans=0;
    for(int i=2;i<=n;++i)
    if(height[i]>ans && ((ind1(sa[i]) && ind2(sa[i-1]))||(ind2(sa[i]) && ind1(sa[i-1]))))
    ans=height[i];
    cout<<ans;
}
int main()
{
    init();
    suffix();
    work();
    return 0;
}