[caioj1469][后缀数组]连续重复子串

问题描述

求两个字符串的最长公共子串。(长度不超过100000)

样例输入

yeshowmuchiloveyoumydearmotherreallyicannotbelieveit
Yeaphowmuchiloveyoumydearmother

样例输出

27

题解

咳咳首先掩饰一发蒟蒻妄图使用KMP&FFT来D本题的尴尬
然后马上发现kmp和FFT的定义都不是酱的。。
考虑后缀数组
将两个串接起来，跑一遍后缀数组。然后比较height值就ok了
这么简单什么鬼
当个模板用吧

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
int a[210000],Rsort[210000],Rank[210000],sa1[210000],sa2[210000];
int tt[210000];
char s[110000];
void get_sa(int len,int m)
{
    memcpy(Rank,a,sizeof(Rank));
    memset(Rsort,0,sizeof(Rsort));
    for(int i=1;i<=len;i++)Rsort[Rank[i]]++;
    for(int i=1;i<=m;i++)Rsort[i]+=Rsort[i-1];
    for(int i=len;i>=1;i--)sa1[Rsort[Rank[i]]--]=i;
    int ln=1,p=0;
    while(p<len)
    {
        int k=0;
        for(int i=len-ln+1;i<=len;i++)sa2[++k]=i;
        for(int i=1;i<=len;i++)if(sa1[i]-ln>0)sa2[++k]=sa1[i]-ln;

        memset(Rsort,0,sizeof(Rsort));
        for(int i=1;i<=len;i++)Rsort[Rank[i]]++;
        for(int i=1;i<=m;i++)Rsort[i]+=Rsort[i-1];
        for(int i=len;i>=1;i--)sa1[Rsort[Rank[sa2[i]]]--]=sa2[i];

        memcpy(tt,Rank,sizeof(tt));
        Rank[sa1[1]]=1;p=1;
        for(int i=2;i<=len;i++)
        {
            if(tt[sa1[i]]!=tt[sa1[i-1]] || tt[sa1[i]+ln]!=tt[sa1[i-1]+ln])p++;
            Rank[sa1[i]]=p;
        }
        ln*=2;m=p;
    }
}
int height[210000];
void get_he(int len)
{
    int j,k=0;
    for(int i=1;i<=len;i++)
    {
        j=sa1[Rank[i]-1];
        if(k!=0)k--;
        while(a[i+k]==a[j+k])k++;
        height[Rank[i]]=k;
    }
}
int main()
{
    int maxn=0;
    scanf("%s",s+1);
    int lenx=strlen(s+1);
    for(int i=1;i<=lenx;i++){a[i]=s[i]-'A'+1;maxn=max(maxn,a[i]);}
    a[lenx+1]=1;
    scanf("%s",s+1);
    int leny=strlen(s+1);
    for(int i=1;i<=leny;i++){a[i+lenx+1]=s[i]-'A'+1;maxn=max(maxn,a[i]);}
    int len=lenx+leny+1;
    get_sa(len,maxn);
    get_he(len);
    int ans=0;
    for(int i=2;i<=len;i++)
    {
        if((sa1[i-1]>lenx+1&&sa1[i]<=lenx) || (sa1[i-1]<=lenx&&sa1[i]>lenx+1))ans=max(ans,height[i]);
    }
    printf("%d\n",ans);
    return 0;
}