[Longest Common Substring II] SPOJ - LCS2

题意:

求n个字符串的最长公共子串的长度

题解:

这个题可以类比只有两个字符串的情况来做。
我们仍然选择对第一个字符串建立SAM,然后对于其他串在SAM上匹配,过程中记录一下每个字符串到达SAM中某个节点时的长度,对于每个状态取 m i n min min,然后从最小值中选取最大值即为答案。

A C AC AC代码:

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 2e5+50;
const int INF = 0x3f3f3f3f;
char s[MAXN],ss[MAXN];
struct Suffix{
    int nxt[MAXN][26],fa[MAXN],len[MAXN];
    int c[MAXN],a[MAXN],l[MAXN],res[MAXN];
    int last=1,tot=1;
    inline void Insert(int x){
        int p = last,np = ++tot;
        last = np,len[np] = len[p] + 1;
        for(;p && !nxt[p][x];p=fa[p]) nxt[p][x] = np;
        if(!p) fa[np] = 1;
        else {
            int q = nxt[p][x];
            if(len[p]+1 == len[q]) fa[np] = q;
            else {
                int nq = ++tot;
                len[nq] = len[p] + 1;
                memcpy(nxt[nq],nxt[q],sizeof(nxt[q]));
                fa[nq] = fa[q];
                fa[q] = fa[np] = nq;
                for(;nxt[p][x]==q;p=fa[p])
                    nxt[p][x] = nq;
            }
        }
    }
    inline void Sort(int n){
        for(int i=1;i<=tot;i++) c[len[i]]++;
        for(int i=1;i<=n;i++) c[i]+=c[i-1];
        for(int i=1;i<=tot;i++) a[c[len[i]]--]=i;
        memset(res,INF,sizeof(res));
    }
    inline void Solve(int n){
        int p = 1,ans = 0;
        memset(l,0,sizeof(l));
        for(int i=1;i<=n;i++){
            int x = ss[i]-'a';
            if(nxt[p][x]) p = nxt[p][x],ans ++;
            else {
                for(;p && !nxt[p][x];p=fa[p]);
                if(!p) p = 1,ans = 0;
                else ans = len[p]+1, p = nxt[p][x];
            }
            l[p] = max(l[p],ans);
        }
        for(int i=tot;i;i--){
            int x = a[i];
            res[x] = min(res[x],l[x]);
            if(l[x] && fa[x])
                l[fa[x]] = len[fa[x]];//如果它能够匹配到,那么它的fa也应该能匹配到
        }
    }
    inline void Work(){
        int ans = 0;
        for(int i=1;i<=tot;i++)
            ans = max(ans,res[i]);
        printf("%d\n",ans);
    }
}SAM;
int main(){
    scanf("%s",s+1);
    int m = strlen(s+1);
    for(int i=1;i<=m;i++) SAM.Insert(s[i]-'a');
    SAM.Sort(m);
    while(~scanf("%s",ss+1)){
        int n = strlen(ss+1);
        SAM.Solve(n);
    }
    SAM.Work();
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
### 回答1: 最长公共子序列(Longest Common Subsequence)指的是在两个序列中找到最长的公共子序列,这个公共子序列可以不连续,但是需要保持相对顺序不变。例如,对于序列ABCD和ACDFG,它们的最长公共子序列是ACD。 ### 回答2: 最长公共子序列(Longest Common Subsequence,简称LCS)是指在给定多个序列中,找到最长的一个子序列,该子序列同时出现在这些序列中,并且其他元素的相对顺序保持一致。 举个例子,假设有两个序列A和B,A为[1, 2, 3, 4, 5],B为[2, 4, 5, 6]。它们的一个最长公共子序列是[2, 4, 5],该子序列同时存在于A和B中。 求解LCS的问题可以用动态规划的方法来解决。我们可以构建一个二维数组dp,其中dp[i][j]表示序列A的前i个元素和序列B的前j个元素的LCS长度。那么dp[i][j]可以通过以下方式得到: 1. 如果A[i]等于B[j],则dp[i][j]等于dp[i-1][j-1] + 1; 2. 如果A[i]不等于B[j],则dp[i][j]等于max(dp[i-1][j], dp[i][j-1])。 通过填充整个dp数组,最终可以得到序列A和序列B的LCS长度。要找到具体的LCS序列,则可以通过反向遍历dp数组进行构建。 LCS问题在字符串处理、DNA序列匹配、版本控制等领域都有广泛的应用。其时间复杂度为O(m*n),其中m和n分别为序列A和序列B的长度。 ### 回答3: 最长公共子序列(Longest Common Subsequence)是一个经典的计算机科学问题。给定两个序列S和T,我们要找出它们之间最长的公共子序列。 子序列是从给定序列中按顺序选择几个元素而组成的序列。而公共子序列指的是同时是序列S和T的子序列的序列。 为了解决这个问题,可以使用动态规划的方法。我们可以定义一个二维数组dp,其中dp[i][j]表示序列S的前i个元素和序列T的前j个元素之间的最长公共子序列的长度。 接下来,我们可以使用以下递推关系来填充dp数组: 如果S[i]等于T[j],则dp[i][j] = dp[i-1][j-1] + 1; 如果S[i]不等于T[j],则dp[i][j] = max(dp[i-1][j], dp[i][j-1])。 最后,我们可以通过查看dp[S.length()][T.length()]来得到最长公共子序列的长度。 此外,我们也可以用回溯法来还原最长公共子序列本身。我们可以从dp[S.length()][T.length()]开始,如果S[i]等于T[j],则将S[i]添加到结果序列中,并向左上方移动,即i = i-1,j = j-1。如果S[i]不等于T[j],则根据dp数组的值选择向上(i = i-1)或向左(j = j-1)移动。 总之,最长公共子序列问题是一个经典的计算机科学问题,可以使用动态规划的方法解决。我们可以通过构建二维dp数组来计算最长公共子序列的长度,并可以使用回溯法来还原它本身。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值