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

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

传送门:POJ 3415: Common Substrings

题意:

公共子串:统计A和B长度不小于K的公共子串个数。

思路: 一开始想当然的求LCP数组,把LCP遍历一遍,大于K的,ans不断累加ans += lcp[i] - k + 1,但这只是一部分解,因为所求的只是相邻后缀的最长公共前缀。

一种直观的做法时,两个for循环,对于当前位置属于S的LCP时,向下找属于T的后缀,把每个后缀的最长对应公共前缀加一遍,但TLE了,此处需要利用单调性来简化时间复杂度。采用记忆化手段,用栈。

自己写了一种,但各别样例过不了,并不清楚原因。此处代码参考hankcs的代码,具体思路如下:对应于S的lcp,先不管三七二十一,贡献值都以最大的加,但我们知道,中间如果lcp变小,意味着后缀的公共前缀在不断减小,此时需要按照最小的算,只需要在贡献值上,删除原先多加的即可。针对T的lcp,同样处理。

代码如下:

import java.io.BufferedReader;
import java.io.File;
import java.io.FileInputStream;
import java.io.IOException;
import java.io.InputStream;
import java.io.InputStreamReader;
import java.io.PrintWriter;
import java.util.Arrays;
import java.util.Comparator;
import java.util.StringTokenizer;

public class Main{

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

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

    class SuffixArray{
        int k = 1;
        int n;
        Integer[] sa;
        int[] rank, tmp, 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
            int h = 0;
            for (int i = 0; i <= n; ++i) rank[sa[i]] = i;
            for (int i = 0; i < n; ++i) {
                int j = sa[rank[i] - 1];
                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;
            }
        }

        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 read() {
        while (true) {
            int k = ni();
            if (k == 0) break;
            String s = ns();
            String t = ns();
            SuffixArray sa = new SuffixArray((s + '$' + t).toCharArray());
            long ans = 0;
            ans += solve(true, s.length(), k, sa);
            ans += solve(false, s.length(), k, sa);
            out.println(ans);
        }
    }

    int[][] stack;
    int top;
    long sum;
    long solve(boolean isS, int sLen, int k, SuffixArray sa) {
        long ans = 0;
        sum = 0;
        top = 0;
        stack = new int[sa.n + 16][2];
        for (int i = 0; i < sa.n; ++i) {
            if (sa.lcp[i] < k) {
                top = 0;
                sum = 0;
            }
            else {
                int size = 0;
                if ((isS && sa.sa[i] < sLen) || (!isS && sa.sa[i] > sLen)) {
                    size ++;
                    sum += sa.lcp[i] - k + 1;
                }
                while (top > 0 && sa.lcp[i] < stack[top - 1][0]) {
                    --top;
                    sum -= stack[top][1] * (stack[top][0] - sa.lcp[i]);
                    size += stack[top][1];
                }
                if (size > 0) {
                    stack[top][0] = sa.lcp[i];
                    stack[top][1] = size;
                    top ++;
                }

                if ((isS && sa.sa[i + 1] > sLen) || (!isS && sa.sa[i + 1] < sLen)) {
                    ans += sum;
                }
            }
        }
        return ans;
    }

    FastScanner in;
    PrintWriter out;

    void run() throws IOException {
        boolean oj;
        try {
            oj = ! System.getProperty("user.dir").equals("F:\\java_workspace\\leetcode");
        } catch (Exception e) {
            oj = System.getProperty("ONLINE_JUDGE") != null;
        }

        InputStream is = oj ? System.in : new FileInputStream(new File(INPUT));
        in = new FastScanner(is);
        out = new PrintWriter(System.out);
        long s = System.currentTimeMillis();
        read();
        out.flush();
        if (!oj){
            System.out.println("[" + (System.currentTimeMillis() - s) + "ms]");
        }
    }

    public boolean more(){
        return in.hasNext();
    }

    public int ni(){
        return in.nextInt();
    }

    public long nl(){
        return in.nextLong();
    }

    public double nd(){
        return in.nextDouble();
    }

    public String ns(){
        return in.nextString();
    }

    public char nc(){
        return in.nextChar();
    }

    class FastScanner {
        BufferedReader br;
        StringTokenizer st;
        boolean hasNext;

        public FastScanner(InputStream is) throws IOException {
            br = new BufferedReader(new InputStreamReader(is));
            hasNext = true;
        }

        public String nextToken() {
            while (st == null || !st.hasMoreTokens()) {
                try {
                    st = new StringTokenizer(br.readLine());
                } catch (Exception e) {
                    hasNext = false;
                    return "##";
                }
            }
            return st.nextToken();
        }

        String next = null;
        public boolean hasNext(){
            next = nextToken();
            return hasNext;
        }

        public int nextInt() {
            if (next == null){
                hasNext();
            }
            String more = next;
            next = null;
            return Integer.parseInt(more);
        }

        public long nextLong() {
            if (next == null){
                hasNext();
            }
            String more = next;
            next = null;
            return Long.parseLong(more);
        }

        public double nextDouble() {
            if (next == null){
                hasNext();
            }
            String more = next;
            next = null;
            return Double.parseDouble(more);
        }

        public String nextString(){
            if (next == null){
                hasNext();
            }
            String more = next;
            next = null;
            return more;
        }

        public char nextChar(){
            if (next == null){
                hasNext();
            }
            String more = next;
            next = null;
            return more.charAt(0);
        }
    }
}

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

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏拂晓风起

java 将数组排序 指这种数组:[] ArrayList之类的请用Collection.sort

991
来自专栏机器学习入门

POJ 刷题系列:1573. Robot Motion

题意: 一张地图包含N,S,W,E的指令,从指定起点[1, p]出发,是否有一条路径能够走到地图边缘,有则输出路径数,无说明走入了死循环,输出走入循环前和循环...

20710
来自专栏闻道于事

Java面向对象之抽象类,接口

抽象类: 含有抽象方法的类被声明为抽象类 抽象方法由子类去实现 含有抽象方法的类必须被声明为抽象类 抽象类被子类继承,子类(如果不是抽象类)必须重写抽象类中...

2895
来自专栏java学习

面试题53(考察求职者对String声明变量在jvm中的存储方法)

(单选题) 1、有如下一段代码,请选择其运行结果() public class StringDemo{ private static final Stri...

2873
来自专栏大数据钻研

44 个 JavaScript 变态题解析

当初笔者做这套题的时候不仅怀疑智商, 连人生都开始怀疑了…. 不过, 对于基础知识的理解是深入编程的前提. 让我们一起来看看这些变态题到底变态不变态吧! 第1题...

3366
来自专栏禁心尽力

Java Collection开发技巧

Java Collection(集合) ? 集合中的一些技巧: 通过Collections类的静态方法,可以对集合进行一些操作 1 java....

21510
来自专栏机器学习入门

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

挑战程序竞赛系列(73):4.7高度数组(3) 传送门:POJ 3729: Facer’s String 题意: 公共子串: 给出两个字符串A,B,求...

2048
来自专栏salesforce零基础学习

salesforce 零基础开发入门学习(二)变量基础知识,集合,表达式,流程控制语句

salesforce如果简单的说可以大概分成两个部分:Apex,VisualForce Page. 其中Apex语言和java很多的语法类似,今天总结的是一些简...

2867
来自专栏微信公众号:Java团长

Java集合类详解

上述类图中,实线边框的是实现类,比如ArrayList,LinkedList,HashMap等,折线边框的是抽象类,比如AbstractCollectio...

981
来自专栏我是攻城师

关于Java里面的字符串拼接,你了解多少?

字符串拼接是我们日常开发中很常见的操作,虽然常见,但要是使用不当的的话,很有可能让你的程序处理效率降低一大半,所以我们有必要来重新了解一下Java里面的字符串操...

1453

扫码关注云+社区