/****************************************************************************** * Compilation: javac SuffixArray.java * Execution: java SuffixArray < input.txt * Dependencies: StdIn.java StdOut.java * Data files: https://algs4.cs.princeton.edu/63suffix/abra.txt * * A data type that computes the suffix array of a string. * * % java SuffixArray < abra.txt * i ind lcp rnk select * --------------------------- * 0 11 - 0 "!" * 1 10 0 1 "A!" * 2 7 1 2 "ABRA!" * 3 0 4 3 "ABRACADABRA!" * 4 3 1 4 "ACADABRA!" * 5 5 1 5 "ADABRA!" * 6 8 0 6 "BRA!" * 7 1 3 7 "BRACADABRA!" * 8 4 0 8 "CADABRA!" * 9 6 0 9 "DABRA!" * 10 9 0 10 "RA!" * 11 2 2 11 "RACADABRA!" * * See SuffixArrayX.java for an optimized version that uses 3-way * radix quicksort and does not use the nested class Suffix. * ******************************************************************************/ package edu.princeton.cs.algs4; import java.util.Arrays; /** * The {@code SuffixArray} class represents a suffix array of a string of * length n. * It supports the selecting the ith smallest suffix, * getting the index of the ith smallest suffix, * computing the length of the longest common prefix between the * ith smallest suffix and the i-1st smallest suffix, * and determining the rank of a query string (which is the number * of suffixes strictly less than the query string). *

* This implementation uses a nested class {@code Suffix} to represent * a suffix of a string (using constant time and space) and * {@code Arrays.sort()} to sort the array of suffixes. * The index and length operations takes constant time * in the worst case. The lcp operation takes time proportional to the * length of the longest common prefix. * The select operation takes time proportional * to the length of the suffix and should be used primarily for debugging. *

* For alternate implementations of the same API, see * {@link SuffixArrayX}, which is faster in practice (uses 3-way radix quicksort) * and uses less memory (does not create {@code Suffix} objects) * and SuffixArrayJava6.java, * which relies on the constant-time substring extraction method that existed * in Java 6. *

* For additional documentation, see Section 6.3 of * Algorithms, 4th Edition by Robert Sedgewick and Kevin Wayne. */ public class SuffixArray { private Suffix[] suffixes; /** * Initializes a suffix array for the given {@code text} string. * @param text the input string */ public SuffixArray(String text) { int n = text.length(); this.suffixes = new Suffix[n]; for (int i = 0; i < n; i++) suffixes[i] = new Suffix(text, i); Arrays.sort(suffixes); } private static class Suffix implements Comparable { private final String text; private final int index; private Suffix(String text, int index) { this.text = text; this.index = index; } private int length() { return text.length() - index; } private char charAt(int i) { return text.charAt(index + i); } public int compareTo(Suffix that) { if (this == that) return 0; // optimization int n = Math.min(this.length(), that.length()); for (int i = 0; i < n; i++) { if (this.charAt(i) < that.charAt(i)) return -1; if (this.charAt(i) > that.charAt(i)) return +1; } return this.length() - that.length(); } public String toString() { return text.substring(index); } } /** * Returns the length of the input string. * @return the length of the input string */ public int length() { return suffixes.length; } /** * Returns the index into the original string of the ith smallest suffix. * That is, {@code text.substring(sa.index(i))} is the ith smallest suffix. * @param i an integer between 0 and n-1 * @return the index into the original string of the ith smallest suffix * @throws java.lang.IllegalArgumentException unless {@code 0 <= i < n} */ public int index(int i) { if (i < 0 || i >= suffixes.length) throw new IllegalArgumentException(); return suffixes[i].index; } /** * Returns the length of the longest common prefix of the ith * smallest suffix and the i-1st smallest suffix. * @param i an integer between 1 and n-1 * @return the length of the longest common prefix of the ith * smallest suffix and the i-1st smallest suffix. * @throws java.lang.IllegalArgumentException unless {@code 1 <= i < n} */ public int lcp(int i) { if (i < 1 || i >= suffixes.length) throw new IllegalArgumentException(); return lcpSuffix(suffixes[i], suffixes[i-1]); } // longest common prefix of s and t private static int lcpSuffix(Suffix s, Suffix t) { int n = Math.min(s.length(), t.length()); for (int i = 0; i < n; i++) { if (s.charAt(i) != t.charAt(i)) return i; } return n; } /** * Returns the ith smallest suffix as a string. * @param i the index * @return the i smallest suffix as a string * @throws java.lang.IllegalArgumentException unless {@code 0 <= i < n} */ public String select(int i) { if (i < 0 || i >= suffixes.length) throw new IllegalArgumentException(); return suffixes[i].toString(); } /** * Returns the number of suffixes strictly less than the {@code query} string. * We note that {@code rank(select(i))} equals {@code i} for each {@code i} * between 0 and n-1. * @param query the query string * @return the number of suffixes strictly less than {@code query} */ public int rank(String query) { int lo = 0, hi = suffixes.length - 1; while (lo <= hi) { int mid = lo + (hi - lo) / 2; int cmp = compare(query, suffixes[mid]); if (cmp < 0) hi = mid - 1; else if (cmp > 0) lo = mid + 1; else return mid; } return lo; } // compare query string to suffix private static int compare(String query, Suffix suffix) { int n = Math.min(query.length(), suffix.length()); for (int i = 0; i < n; i++) { if (query.charAt(i) < suffix.charAt(i)) return -1; if (query.charAt(i) > suffix.charAt(i)) return +1; } return query.length() - suffix.length(); } /** * Unit tests the {@code SuffixArray} data type. * * @param args the command-line arguments */ public static void main(String[] args) { String s = StdIn.readAll().replaceAll("\\s+", " ").trim(); SuffixArray suffix = new SuffixArray(s); // StdOut.println("rank(" + args[0] + ") = " + suffix.rank(args[0])); StdOut.println(" i ind lcp rnk select"); StdOut.println("---------------------------"); for (int i = 0; i < s.length(); i++) { int index = suffix.index(i); String ith = "\"" + s.substring(index, Math.min(index + 50, s.length())) + "\""; assert s.substring(index).equals(suffix.select(i)); int rank = suffix.rank(s.substring(index)); if (i == 0) { StdOut.printf("%3d %3d %3s %3d %s\n", i, index, "-", rank, ith); } else { int lcp = suffix.lcp(i); StdOut.printf("%3d %3d %3d %3d %s\n", i, index, lcp, rank, ith); } } } } /****************************************************************************** * Copyright 2002-2022, Robert Sedgewick and Kevin Wayne. * * This file is part of algs4.jar, which accompanies the textbook * * Algorithms, 4th edition by Robert Sedgewick and Kevin Wayne, * Addison-Wesley Professional, 2011, ISBN 0-321-57351-X. * http://algs4.cs.princeton.edu * * * algs4.jar is free software: you can redistribute it and/or modify * it under the terms of the GNU General Public License as published by * the Free Software Foundation, either version 3 of the License, or * (at your option) any later version. * * algs4.jar is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU General Public License for more details. * * You should have received a copy of the GNU General Public License * along with algs4.jar. If not, see http://www.gnu.org/licenses. ******************************************************************************/