diff options
| -rw-r--r-- | geometry/geometry.tex | 8 | ||||
| -rw-r--r-- | math/math.tex | 1 | ||||
| -rw-r--r-- | string/suffixArray.cpp | 36 |
3 files changed, 27 insertions, 18 deletions
diff --git a/geometry/geometry.tex b/geometry/geometry.tex index b0d8454..d3e1671 100644 --- a/geometry/geometry.tex +++ b/geometry/geometry.tex @@ -7,14 +7,14 @@ \sourcecode{geometry/closestPair.cpp} \end{algorithm} -\begin{algorithm}{Konvexe Hülle} +\begin{algorithm}{Konvexehülle} \begin{methods} \method{convexHull}{berechnet Konvexehülle}{n\*\log(n)} \end{methods} \begin{itemize} - \item Konvexehülle gegen den Uhrzeigersinn Sortiert + \item Konvexehülle gegen den Uhrzeigersinn sortiert \item nur Eckpunkte enthalten(für alle Punkte = im CCW Test entfernen) - \item Erster und Letzter Punkt sind identisch + \item erster und letzter Punkt sind identisch \end{itemize} \sourcecode{geometry/convexHull.cpp} \end{algorithm} @@ -24,7 +24,7 @@ \begin{methods} \method{antipodalPoints}{berechnet antipodale Punkte}{n} \end{methods} - \textbf{WICHTIG:} Punkte müssen gegen den Uhrzeigersinn Sortiert sein und konvexes Polygon bilden! + \textbf{WICHTIG:} Punkte müssen gegen den Uhrzeigersinn sortiert sein und konvexes Polygon bilden! \sourcecode{geometry/antipodalPoints.cpp} \end{algorithm} diff --git a/math/math.tex b/math/math.tex index 3ffbb41..8ccc55e 100644 --- a/math/math.tex +++ b/math/math.tex @@ -7,6 +7,7 @@ \end{itemize}
\sourcecode{math/longestIncreasingSubsequence.cpp}
\end{algorithm}
+\columnbreak
\begin{algorithm}{Zykel Erkennung}
\begin{methods}
diff --git a/string/suffixArray.cpp b/string/suffixArray.cpp index 9e6960e..54deb84 100644 --- a/string/suffixArray.cpp +++ b/string/suffixArray.cpp @@ -1,30 +1,38 @@ +constexpr int MAX_CHAR = 256; struct SuffixArray { int n; vector<int> SA, LCP; vector<vector<int>> P; - SuffixArray(const string& s) : n(sz(s)), SA(n), LCP(n) { - vector<ll> L(n); - P.assign(__lg(n - 1) + 2, vector<int>(n, 1)); - P[0] = vector<int>(all(s)); + SuffixArray(const string& s) : n(sz(s)), SA(n), LCP(n), + P(__lg(n - 1) + 2, vector<int>(n)) { + P[0].assign(all(s)); iota(all(SA), 0); + sort(all(SA), [&](int a, int b) {return s[a] < s[b];}); + vector<int> x(n); for (int k = 1, c = 1; c < n; k++, c *= 2) { - for (int i = 0; i < n; i++) L[i] = (ll)P[k - 1][i] << 32; - for (int i = 0; i+c < n; i++) L[i] |= P[k - 1][i + c]; - sort(all(SA), [&](int a, int b) {return L[a] < L[b];}); + iota(all(x), n - c); + for (int ptr = c; int i : SA) if (i >= c) x[ptr++] = i - c; + + vector<int> cnt(k == 1 ? MAX_CHAR : n); + for (int i : P[k-1]) cnt[i]++; + partial_sum(all(cnt), begin(cnt)); + for (int i : x | views::reverse) SA[--cnt[P[k-1][i]]] = i; + + auto p = [&](int i) {return i < n ? P[k-1][i] : -1;}; for (int i = 1; i < n; i++) { - P[k][SA[i]] = P[k][SA[i-1]] + (L[SA[i]] != L[SA[i-1]]); + int a = SA[i-1], b = SA[i]; + P[k][b] = P[k][a] + (p(a) != p(b) || p(a+c) != p(b+c)); }} for (int i = 1; i < n; i++) LCP[i] = lcp(SA[i-1], SA[i]); } - // x and y are text-indices, not SA-indices. - int lcp(int x, int y) { - int ret = 0; + int lcp(int x, int y) {//x & y are text-indices, not SA-indices if (x == y) return n - x; - for (int i = sz(P) - 1; i >= 0 && max(x, y) + ret < n; i--) { - if (P[i][x + ret] == P[i][y + ret]) ret += 1 << i; + int res = 0; + for (int i = sz(P) - 1; i >= 0 && max(x, y) + res < n; i--) { + if (P[i][x + res] == P[i][y + res]) res |= 1 << i; } - return ret; + return res; } }; |