题目地址:http://poj.org/problem?id=2774

题目:


给出两个字符串,求他们的最长公共子串的长度

 

解题思路:


后缀数组经典题。

两个字符串通过一个特殊字符连接在一起,求这个新的字符串的sa数组和height数组

求最长公共子串转化为求新字符串的两个排名相邻的后缀的LCP问题(因为是LCP,所以这两个后缀的排名一定是紧挨着的,不需要st表)

并且需要保证一个后缀的开始下标在给出的第一个字符串内,另一个后缀的开始下标在给出的第二个字符串内。

多组输入!

 

ac代码:


#include <iostream>
#include <stdio.h>
#include <string.h>
using namespace std;
const int maxn = 200020;
typedef long long ll;
//不能声明ws,保留字
int sa[maxn], wv[maxn], wss[maxn], wa[maxn], wb[maxn], r[maxn];
char s1[maxn/2], s2[maxn/2];
bool cmp(int *r, int a, int b, int l)
{
    return r[a] == r[b] && r[a + l] == r[b + l];
}
//O(nlogn)读入下标从0开始
void get_sa(int *r, int *sa, int n, int m)
{
    int *x=wa, *y=wb;
    int p =0, i, j;
    for(i = 0; i < m; i++) wss[i] = 0;
    for(i = 0; i < n; i++) wss[ x[i]=r[i] ]++;
    for(i = 1; i <= m; i++) wss[i] += wss[i - 1];
    for(i = n - 1; i >= 0; i--) sa[--wss[x[i]]] = i;
    for(j = 1, p = 1; p < n; j *= 2, m = p)
    {
        //对第二关键字排序
        for(p = 0, i = n - j; i < n; i++) // [n-j,n)没有内容
            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]]; //排名为i的第二关键字对应的第一关键字的排名,x此时相当于rankk,y相当于第二关键字的sa
        for(i = 0; i < m; i++) wss[i] = 0;
        for(i = 0; i < n; i++) wss[wv[i]]++;
        for(i = 1; i <= m; i++) wss[i] += wss[i - 1];
        for(i = n - 1; i >= 0; i--) sa[--wss[wv[i]]] = y[i];
        //相同的字符串排名相同
        swap(x,y);
        for(i = 1, p = 1, x[sa[0]] = 0; i < n; i++)
            x[sa[i]] = cmp(y, sa[i-1], sa[i], j) ? p - 1 : p++;
    }
}
//O(n)
int rankk[maxn], height[maxn];
void get_height(int n)
{
    int k = 0;
    for(int i = 1; i <= n; i++) rankk[sa[i]] = i;
    for(int i = 0; i < n; i++)
    {
        k ? k-- : 0;//根据性质height[rank[i]] ≥ (height[rank[i-1]] -1)
        int j = sa[rankk[i] - 1];//上一名的开始下标
        while(r[i + k] == r[j + k]) k++;
        height[rankk[i]] = k;
    }
}
bool judge(int i, int len1)
{
    if(0 <= sa[i-1] && sa[i-1] < len1 && len1 < sa[i])
        return  true;
    else if(0 <= sa[i] && sa[i] < len1 && len1 < sa[i-1])
        return true;
    else return  false;
}
int main() {
    //freopen("/Users/zhangkanqi/Desktop/11.txt","r",stdin);
    while (~scanf("%s", s1))
    {
    int n = 0;
    int len1 = strlen(s1);
    for (int i = 0; i < len1; i++)
        r[n++] = (int) s1[i];
    r[n++] = '*';
    scanf("%s", s2);
    int len2 = strlen(s2);
    for (int i = 0; i < len2; i++)
        r[n++] = (int) s2[i];
    r[n] = 0;
    get_sa(r, sa, n + 1, 255);
    get_height(n);
    int ans = 0;
    for (int i = 1; i <= n; i++)
        if (height[i] > ans && judge(i, len1))
            ans = height[i];
    printf("%d\n", ans);
    }
    return 0;
}