最长公共子串(后缀数组+LCP)

求两个串的最长公共子串,我们可以将两个串中间插入一个两个串都没有的字符串起来,然后计算出他的后缀数组(sa)和高度数组(lcp),不在同一串中的sa中相邻后缀的最长公共前缀(即高度数组的最大值)就是两个串的最长公共子串。

#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;


const int maxn = 200000+5;
const int oo = 0x3f3f3f3f;
typedef long long ll;
#define FOR(i,n) for(int i=0;i<n;i++)
#define clr(a,x) memset(a,x,sizeof a)

char S[maxn];
char T[maxn];
int ranks[maxn];
int sa[maxn];
int lcp[maxn];
int tmp[maxn];
int k = 0;
int n;
bool compare(int i, int j) {
    if (ranks[i] != ranks[j]) return ranks[i] < ranks[j];
    else {
        int ri = i + k <= n ? ranks[i + k] : -1;
        int rj = j + k <= n ? ranks[j + k] : -1;
        return ri < rj;
    }
}
void calc_sa(char *s,int *sa) {
    for (int i = 0; i <= n; i++) {
        sa[i] = i;
        ranks[i] = i < n ? s[i] : -1;
    }
    for (k = 1; k <= n; k *= 2) {
        sort(sa, sa + n + 1, compare);
        tmp[sa[0]] = 0;
        for (int i = 1; i <= n; i++) {
            tmp[sa[i]] = tmp[sa[i - 1]] + (compare(sa[i - 1], sa[i]) ? 1 : 0);
        }
        for (int i = 0; i <= n; i++) {
            ranks[i] = tmp[i];
        }
    }
}

void calc_lcp(char *s, int *sa, int *lcp) {
    int h = 0;
    lcp[0] = 0;
    for (int i = 0; i <= n; i++) {
        int j = sa[ranks[i] + 1];
        if (h > 0) h--;
        for (; j + h < n&&i + h < n; h++) {
            if (s[j + h] != s[i + h]) break;
        }
        lcp[ranks[i]] = h;
    }
}

void solve() {
    int len = strlen(S);
    int m = strlen(T);
    n = len + m + 1;
    S[len] = '\0';
    for (int i = 0; i < m; i++) {
        S[len + i + 1] = T[i];
    }
    S[m + len + 1] = 0;
    int ans = 0;
    calc_sa(S, sa);
    calc_lcp(S, sa, lcp);
    for (int i = 0; i < n; i++) {
        if ((sa[i] < len) != (sa[i + 1] < len)) {
            ans = max(ans, lcp[i]);
        }
    }
    printf("%d\n", ans);
}
int main() {
    while (scanf("%s%s", S, T) != EOF) {
        solve();
    }
    return 0;
}
  • 2
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值