URAL - 1517 Freedom of Choice (最长公共子串)

🎑 🎑 🎑

求两个字符串的最长公共子串。
利用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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值