[SPOJ-LCS2]Longest Common Substring II

vjudge

题意

求多个串的最长公共子串

sol

首先拿第一个串建SAM。
然后对于每一个串，在SAM上跑匹配。
那么最终就是要把多个串跑出匹配的信息合并起来。
具体来说，对每个串记录一下在每个状态上的最长长度，然后同一个状态上对每个串的最长长度取个\(min\)，最后再把每个状态取个\(max\)。

code

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 2e5+5;
int n,last=1,tot=1,tr[N][26],fa[N],len[N],t[N],a[N],sum[N],ans[N],Ans;
char s[N];
void extend(int c)
{
    int v=last,u=++tot;last=u;
    len[u]=len[v]+1;
    while (v&&!tr[v][c]) tr[v][c]=u,v=fa[v];
    if (!v) fa[u]=1;
    else
    {
        int x=tr[v][c];
        if (len[x]==len[v]+1) fa[u]=x;
        else
        {
            int y=++tot;
            memcpy(tr[y],tr[x],sizeof(tr[y]));
            fa[y]=fa[x];fa[x]=fa[u]=y;len[y]=len[v]+1;
            while (v&&tr[v][c]==x) tr[v][c]=y,v=fa[v];
        }
    }
}
int main()
{
    scanf("%s",s+1);n=strlen(s+1);
    for (int i=1;i<=n;++i) extend(s[i]-'a');
    for (int i=1;i<=tot;++i) ++t[len[i]];
    for (int i=1;i<=tot;++i) t[i]+=t[i-1];
    for (int i=1;i<=tot;++i) a[t[len[i]]--]=i;
    for (int i=1;i<=tot;++i) ans[i]=len[i];
    while (scanf("%s",s+1)!=EOF)
    {
        n=strlen(s+1);
        memset(sum,0,sizeof(sum));
        for (int i=1,now=1,cnt=0;i<=n;++i)
        {
            int c=s[i]-'a';
            if (tr[now][c]) ++cnt,now=tr[now][c];
            else
            {
                while (now&&!tr[now][c]) now=fa[now];
                if (!now) cnt=0,now=1;
                else cnt=len[now]+1,now=tr[now][c];
            }
            sum[now]=max(sum[now],cnt);
        }
        for (int i=tot;i;--i) sum[fa[a[i]]]=max(sum[fa[a[i]]],sum[a[i]]);
        for (int i=1;i<=tot;++i) ans[i]=min(ans[i],sum[i]);
    }
    for (int i=1;i<=tot;++i) Ans=max(Ans,ans[i]);
    printf("%d\n",Ans);return 0;
}

转载于:https://www.cnblogs.com/zhoushuyu/p/8660359.html