SP1812 LCS2

一、题目

点此看题

二、解法

可以先把这道题写了，只有两个串，是这道题的简化版。

我们先对于第一个串建出后缀自动机，然后把其他串丢进去匹配，大致方法如下（实际上第一种情况可以和第二种情况合并，减少代码量）：

• 如果当前点有转移，直接跳，长度++
• 如果没有转移，跳父亲，看它有没有转移，在第一个有的停下，操作即可
• 如果跳出了自动机，回到根节点，长度赋值$0$

自动机上的每个节点存一个最大长度$mx$，我们在匹配的时候把长度拿来更新跳到的点。由于子节点的最大长度是可以转移到父节点的，匹配完之后我们用子树的$mx$更新父亲的$mx$（具体就是把长度排序，注意$mx$不能超过节点的长度），我们对于每个节点再维护一个$mi$，表示历史所有$mx$的最小值，就是$lcp$选在这个点上的长度，最后我们找到$mi$的最大值即是答案。

时间复杂度$O(n)$，贴两个代码，第一个是弱化版的，第二个本题代码。第一个代码匹配是详细讨论了，第二个代码是写的缩略版$qwq$。

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int M = 2000005;
int read()
{
	int x=0,flag=1;char c;
	while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
	while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return x*flag;
}
int n,m;char s[M];
struct node
{
    int len,fa,ch[26];
    node() {memset(ch,0,sizeof ch);len=fa=0;}
};
struct automaton
{
    int cnt,last,ans;node a[M];
    automaton() {cnt=last=1;}
    void add(int c)
    {
        int p=last,np=last=++cnt;
        a[np].len=a[p].len+1;
        for(;p && !a[p].ch[c];p=a[p].fa) a[p].ch[c]=np;
        if(!p) a[np].fa=1;
        else
        {
            int q=a[p].ch[c];
            if(a[q].len==a[p].len+1) a[np].fa=q;
            else
            {
                int nq=++cnt;
                a[nq]=a[q];a[nq].len=a[p].len+1;
                a[q].fa=a[np].fa=nq;
                for(;p && a[p].ch[c]==q;p=a[p].fa) a[p].ch[c]=nq;
            }
        }
    }
    void solve()
    {
        scanf("%s",s);
        n=strlen(s);
        int p=1,len=0;
        for(int i=0;i<n;i++)
        {
            int c=s[i]-'a';
            if(a[p].ch[c])
                p=a[p].ch[c],len++;
            else
            {
                for(;p && !a[p].ch[c];p=a[p].fa);
                if(p) len=a[p].len+1,p=a[p].ch[c];
                else p=1,len=0;
            }
            ans=max(ans,len);
        }
        printf("%d\n",ans);
    }
}Sam;
int main()
{
    scanf("%s",s);
    n=strlen(s);
    for(int i=0;i<n;i++)
        Sam.add(s[i]-'a');
    Sam.solve();
}

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int M = 2000005;
int read()
{
 int x=0,flag=1;char c;
 while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
 while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
 return x*flag;
}
int n,ans;char s[M];
struct node
{
    int len,fa,ch[26];
    node() {memset(ch,0,sizeof ch);len=fa=0;}
};
struct automaton
{
    int cnt,last,c[M],A[M],mx[M],mi[M];node a[M];
    automaton() {cnt=last=1;}
    void add(int c)
    {
        int p=last,np=last=++cnt;
        a[np].len=a[p].len+1;
        for(;p && !a[p].ch[c];p=a[p].fa) a[p].ch[c]=np;
        if(!p) a[np].fa=1;
        else
        {
            int q=a[p].ch[c];
            if(a[q].len==a[p].len+1) a[np].fa=q;
            else
            {
                int nq=++cnt;
                a[nq]=a[q];a[nq].len=a[p].len+1;
                a[q].fa=a[np].fa=nq;
                for(;p && a[p].ch[c]==q;p=a[p].fa) a[p].ch[c]=nq;
            }
        }
    }
    void sort()
    {
        memset(mi,0x3f,sizeof mi);
        for(int i=1;i<=cnt;i++) c[a[i].len]++;
        for(int i=1;i<=n;i++) c[i]+=c[i-1];
        for(int i=1;i<=cnt;i++) A[c[a[i].len]--]=i;
    }
    void solve()
    {
        n=strlen(s);
        int p=1,len=0;
        for(int i=0;i<n;i++)
        {
            int c=s[i]-'a';
            for(;p && !a[p].ch[c];p=a[p].fa,len=a[p].len);
            if(p) p=a[p].ch[c],len++;
            else p=1,len=0;
            mx[p]=max(mx[p],len);
        }
        for(int i=cnt;i>=1;i--)
        {
            int x=A[i],y=a[x].fa;
            mx[y]=max(mx[y],min(a[y].len,mx[x]));
            mi[x]=min(mi[x],mx[x]);
            mx[x]=0;
        }
    }
}Sam;
int main()
{
    scanf("%s",s);
    n=strlen(s);
    for(int i=0;i<n;i++)
        Sam.add(s[i]-'a');
    Sam.sort();
    while(~scanf("%s",s))
        Sam.solve();
    for(int i=1;i<=Sam.cnt;i++)
        ans=max(ans,Sam.mi[i]);
    printf("%d\n",ans);
}