poj 2774 求两个字符串的最长公共子串

利用后缀数组求两个字符串的最长公共子串。把两个串通过一个特殊字符隔开，合成一个字符串。

然后求height数组，利用二分查找最长公共前缀 且开头下标分别在原来的两个串中。

代码如下：

#include<iostream>
#include<string.h>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define MAXD 200010
#define MAXL 300
using namespace std;
char r[MAXD];
int Rank[MAXD], height[MAXD];
int sa[MAXD], wa[MAXD], wb[MAXD], WS[MAXL], wv[MAXD]; //WS 的大小是最大字符的大小。其他的是串长大小
//r数组存放字符的值
int cmp(int *p, int x, int y, int l)
{
    return p[x] == p[y] && p[x + l] == p[y + l];
}
void da(int n, int m)      //n为字符串长度，m为字符的最大值
{
    int i, j, p, *x = wa, *y = wb, *t;
    for(i = 0; i < m; i ++)
        WS[i] = 0;
    for(i = 0; i < n; i ++)
        ++ WS[x[i] = r[i]];
    for(i = 1; i < m; i ++)
        WS[i] += WS[i - 1];
    for(i = n - 1; i >= 0; i --)
        sa[-- WS[x[i]]] = i;
    for(p = 1, j = 1; p < n; j *= 2, m = p)
    {
        for(p = 0, i = n - j; i < n; i ++)
            y[p ++] = i;
        for(i = 0; i < n; i ++)
            if(sa[i] >= j)
                y[p ++] = sa[i] - j;
        for(i = 0; i < n; i ++)
            wv[i] = x[y[i]];
        for(i = 0; i < m; i ++)
            WS[i] = 0;
        for(i = 0; i < n; i ++)
            ++ WS[wv[i]];
        for(i = 1; i < m; i ++)
            WS[i] += WS[i - 1];
        for(i = n - 1; i >= 0; i --)
            sa[-- WS[wv[i]]] = y[i];
        for(t = x, x = y, y = t, x[sa[0]] = 0, p = 1, i = 1; i < n; i ++)
            x[sa[i]] = cmp(y, sa[i - 1], sa[i], j) ? p - 1: p ++;
    }
}
void calheight(int n)  //n为串的长度减一
{
    int i, j, k = 0;
    for(i = 1; i <= n; i ++)
        Rank[sa[i]] = i;
    for(i = 0; i < n; height[Rank[i ++]] = k)
        for(k ? -- k : 0, j = sa[Rank[i] - 1]; r[i + k] == r[j + k]; k ++);
}
/*****************以上都是后缀数组模板，下面的是具体问题的处理。******************/
char b[MAXD];
int len1,len;
bool judge(int k,int n)
{
    int i;
    int MM=-1,mm=n;
    for(i=1;i<=n;i++)
    {
        if(height[i]>=k)
        {
            MM=max(MM,sa[i]);
            mm=min(mm,sa[i]);
            if(MM>len1&&mm<len1) return 1;
        }
        else {MM=mm=sa[i];}   // 一定要是sa[i].
    }
    return 0;
}
int main()
{
    int i,j;
    gets(r);
    gets(b);
    len1=strlen(r);
    r[len1]='$';
    r[len1+1]=0;
    strcat(r,b);
    len=strlen(r);
    da(len+1,300);
    calheight(len);
    int MIN=min(len1,len-len1-1);
    int L=0,R=MIN;
    while(L<=R)
    {
        int mid=(L+R)/2;
        if(judge(mid,len))
            L=mid+1;
        else R=mid-1;
    }
    printf("%d\n",R);
    return 0;
}