Hiho 122 后缀数组三·重复旋律3(多个串的最长公共重复子串)

using System;

namespace Hiho
{
    class _122
    {
        static readonly int MAXN = 1000050;
        static char[] ch = new char[MAXN], All = new char[MAXN];
        static int[] SA = new int[MAXN], rank = new int[MAXN], Height = new int[MAXN], tax = new int[MAXN], tp = new int[MAXN], original = new int[MAXN];
        static int n, m;
        static void RSort()
        {
            for (int i = 0; i <= m; i++) tax[i] = 0;
            for (int i = 1; i <= n; i++) tax[rank[tp[i]]]++;
            for (int i = 1; i <= m; i++) tax[i] += tax[i - 1];
            for (int i = n; i >= 1; i--) SA[tax[rank[tp[i]]]--] = tp[i];
        }
        static bool cmp(int[] f, int x, int y, int w) { return f[x] == f[y] && f[x + w] == f[y + w]; }
        static void Suffix()
        {
            //SA
            for (int i = 1; i <= n; i++)
            {
                rank[i] = original[i];
                tp[i] = i;
            }
            m = 2000; RSort();
            for (int w = 1, p = 1, i; p < n; w += w, m = p)
            {
                for (p = 0, i = n - w + 1; i <= n; i++) tp[++p] = i;
                for (i = 1; i <= n; i++) if (SA[i] > w) tp[++p] = SA[i] - w;
                RSort();
                int[] temp = rank;
                rank = tp;
                tp = temp;
                rank[SA[1]] = p = 1;
                for (i = 2; i <= n; i++) rank[SA[i]] = cmp(tp, SA[i], SA[i - 1], w) ? p : ++p;
            }
            int j, k = 0;
            for (int i = 1; i <= n; Height[rank[i++]] = k)
                for (k = k != 0 ? k - 1 : k, j = SA[rank[i] - 1]; original[i + k] == original[j + k]; ++k) ;
        }
        static void Main(string[] args)
        {
            string str = Console.ReadLine();
            int index = str.Length;
            str += "$" + Console.ReadLine();
            n = str.Length;
            for (int i = 0; i < n; i++) original[i + 1] = str[i];
            Suffix();
            int ans = 0;
            for (int i = 1; i < n; ++i)
                if (Height[i] > ans && ((SA[i] > index && SA[i - 1] < index) || (SA[i] < index && SA[i - 1] > index)))
                    ans = Height[i];
            Console.WriteLine(ans);
        }
    }
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值