URAL - 1517 Freedom of Choice （最长公共子串）

小松萘

于 2020-02-22 14:30:43 发布

阅读量98

点赞数

分类专栏：字符串---后缀树

本文链接：https://blog.csdn.net/weixin_44116061/article/details/104443382

版权

字符串---后缀树专栏收录该内容

3 篇文章 0 订阅

订阅专栏

🎑 🎑 🎑

求两个字符串的最长公共子串。
利用height数组的时候，为了防止求出的最长子串是连接两个字符串形成的，我们在两个字符串中之间增加一个没有出现过的字符串；为了防止求得的子串是出自同一个字符串的，我们对sa[ i i]与sa[ i -1 ]的范围做出规定

还发现加了关流就会WA，不知道为什么

int n,m;
char s[MAXN];
int sa[MAXN],height[MAXN],x[MAXN],y[MAXN],c[MAXN];
void get_nxt()
{
    //
    rep(i,m+1) c[i]=0;
    rpp(i,n) ++c[x[i]=s[i]],y[i]=i;
    rpp(i,m) c[i]+=c[i-1];
    for(int i=n;i>=1;--i) sa[c[x[y[i]]]--] = y[i];
    //
    for(int k=1;k<=n;k<<=1)
    {
        int num=0;
        for(int i=n-k+1;i<=n;++i) y[++num]=i;
        rpp(i,n) if(sa[i]>k) y[++num]=sa[i]-k;
        rep(i,m+1) c[i]=0;
        rpp(i,n) ++c[x[i]];
        rpp(i,m) c[i]+=c[i-1];
        for(int i=n;i>=1;--i) sa[c[x[y[i]]]--] = y[i];
        swap(x,y);
        num=1,x[sa[1]]=1;
        for(int i=2;i<=n;++i)
            x[sa[i]]=((y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k])?num:++num);
        if(num==n) break;
        m=num;
    }
    int k=0;
    for (int i=1;i<=n;++i)  
    {
        if (x[i]==1) continue;//第一名height为0 
        if (k) --k;//h[i]>=h[i-1]+1;
        int j=sa[x[i]-1];
        while (j+k<=n && i+k<=n && s[i+k]==s[j+k]) ++k;
        height[x[i]]=k;//h[i]=height[rk[i]];
    }
}
void solve()
{
    int tmp=n;
    scanf("%s",s+1);
    s[n+1]=1;
    scanf("%s",s+n+2);
    n=strlen(s+1),m=122;
    get_nxt();
    int ans=0,res=0;
    for(int i=2;i<=n;++i)
    {
        if(ans<height[i])
            if((sa[i]<=tmp&&sa[i-1]>tmp+1)||(sa[i]>tmp+1&&sa[i-1]<=tmp))//
                ans=height[i],res=sa[i];
    }
    s[res+ans]='\0';
    printf("%s\n",s+res);
}
signed main()
{
    while(cin>>n)  solve();
    return 0;
}