挑战程序竞赛系列(71):4.7高度数组(1)

挑战程序竞赛系列(71):4.7高度数组(1)

传送门:POJ 2217: Secretary

题意:

给定两个字符串S和T。请计算两个字符串最长的公共字符串子串的长度。

LCP概念:在后缀数组中相邻两个后缀的最长公共前缀(Longest Common Prefix)

的确很取巧,相邻位置的最长公共前缀可以由前一个位置的最长公共前缀得出,至少在h-1个位置上不需要再比较,而直接比较后续的字符即可。

代码如下:

import java.io.IOException;
import java.util.Arrays;
import java.util.Comparator;
import java.util.Scanner;

public class Main{

    String INPUT = "./data/judge/201709/P2217.txt";

    static Scanner in;
    public static void main(String[] args) throws IOException {
        in = new Scanner(System.in);
        new Main().run();
    }

    void read() {
        int t = in.nextInt();
        in.nextLine();
        while (t --> 0) {
            String S = in.nextLine();
            String T = in.nextLine();
            StringBuilder sb = new StringBuilder();
            sb.append(S + '$' + T);
            SuffixArray sa = new SuffixArray(sb.toString().toCharArray());
            int max = 0;
            for (int i = 0; i < sb.length(); ++i) {
                if (sa.sa[i] < S.length() != sa.sa[i + 1] < S.length())
                    max = Math.max(max, sa.lcp[i]);
            }
            System.out.println("Nejdelsi spolecny retezec ma delku " + max + ".");
        }
    }

    class SuffixArray{
        int k = 1;
        int n;
        Integer[] sa;
        int[] rank, tmp;
        int[] lcp;

        SuffixArray(char[] cs){
            n    = cs.length;
            sa   = new Integer[n + 1];
            rank = new int[n + 1];
            tmp  = new int[n + 1];
            lcp  = new int[n];
            for (int i = 0; i <= n; ++i) {
                sa[i] = i;
                rank[i] = i < n ? cs[i] : -1;
            }

            for (k = 1; k <= n; k <<= 1) {
                Arrays.sort(sa, cmp);
                tmp[sa[0]] = 0;
                for (int i = 1; i <= n; ++i) {
                    tmp[sa[i]] = tmp[sa[i - 1]] + (cmp.compare(sa[i - 1], sa[i]) < 0 ? 1 : 0);  
                }

                for (int i = 0; i <= n; ++i) {
                    rank[i] = tmp[i];
                }
            }

            // lcp
            for (int i = 0; i <= n; ++i) rank[sa[i]] = i;  //求出每个后缀的当前排名
            int h = 0;
            for (int i = 0; i < n; ++i) {
                int j = sa[rank[i] - 1]; // 后缀i的前一个后缀
                if (h > 0) h--;
                for (;j + h < n && i + h < n; ++h) {
                    if (cs[j + h] != cs[i + h]) break;
                }
                lcp[rank[i] -1] = h;
            }

        }



        private Comparator<Integer> cmp = new Comparator<Integer>() {
            @Override
            public int compare(Integer o1, Integer o2) {
                int i = o1;
                int j = o2;
                if (rank[i] != rank[j]) return rank[i] - rank[j];
                else {
                    int ri = i + k <= n ? rank[i + k] : -1;
                    int rj = j + k <= n ? rank[j + k] : -1;
                    return ri - rj;
                }
            }
        };
    }

    void run() throws IOException {
        read();
    }
}

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏吴伟祥

泛型方法 原

问题:编写一个泛型方法add(),当传入不同数字类型的值时,能够进行加法运算,(如可以传入int、long、float、double类型,但要对传入的值做一定的...

713
来自专栏云霄雨霁

查找----基于有序数组

1330
来自专栏owent

数论模板(个人模板)

若a与n互质(即GCD(a,n) = 1),则a^Ψ(n) = 1 (mod n)a^{\varphi(n)} \equiv 1 \pmod n

934
来自专栏Python疯子

数据分析之numpy

ndarray概述 创建n维数组 接收的是列表类型,所有元素类型必须相同 shape表示各维度大小的元组 dtype表示数组数据类型对象

571
来自专栏Fish

蓝桥杯 幸运数

先说题意。题意有点长,我还是复制粘贴吧。。。 问题描述 幸运数是波兰数学家乌拉姆命名的。它采用与生成素数类似的“筛法”生成。 首先从1开始写出自然数1,...

1976
来自专栏闵开慧

java概念2

1 文件操作 public class FileTest { public static void main(String[] args) throw...

3178
来自专栏zaking's

用js来实现那些数据结构08(链表02-双向链表)

  其实无论在任何语言中,一种数据结构往往会有很多的延伸和变种以应对不同场景的需要。其实前面我们所学过的栈和队列也是可以用链表来实现的。有兴趣的小伙伴可以自己尝...

2726
来自专栏linjinhe的专栏

密码学工具箱

对称加密算法(Symmetric-key_algorithm)是指在加密和解密时使用同一密钥的方式,如AES。

1403
来自专栏机器学习从入门到成神

几种有关排序的常见面试问题

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/sinat_35512245/articl...

672
来自专栏开发与安全

散列表(二):冲突处理的方法之链地址法的实现(哈希查找)

首先需要澄清的一点是,这里讲的是hash table/hash map ,即数据项所存储的表要用数组来实现。 一、链地址法 这种基本思想:将所有哈希地址为i 的...

1890

扫码关注云+社区