POJ 1743 (后缀数组)

思路：
首先构建后缀数组，这里要注意theme可以相差k，如abcdefgh中，d-a == e-b因此可以将原数组相邻元素做差，s[i] = s[i] - s[i+1]，这样就转化为大小为n-1的数组建立后缀树组。
然后利用二分法搜索可能的长度，判断该长度是不是无重叠重复子串的长度。对于每一个长度d，将height中大于等于d的保留，小于的去除，这样就将height分为几组，因为求最大重复子串就是求某个区间中最小的height，因此最长的无重叠子串一定在同一组内，判断是否重叠就是，利用sa数组看它们在实际数组中的位置间隔是否比d大，如果大则不重叠。
后缀数组是套用的吉林大学的模板（网上大多数的解题报告都是套的国家集训队的模板），真心给这些大牛跪了，感觉比国家集训队论文给的代码效率要高，并且简洁到令人发指的地步，只有望其项背的份了，不过简洁的代码太尼玛难懂了，看了整整一天才明白啊有木有！！

#include <iostream>using namespace std;const int N = 20050;int s[N];                         //store numbersint n, sa[N], height[N], rank[N], sum[N], tsa[N];inline int min(int a, int b){return a < b ? a : b;}inline int max(int a, int b){return a > b ? a : b;}void makesa(int na){                       // O(N * log N) int i, j, len; //n = strlen(s) + 1;//na = (n < 256 ? 256 : n); memset(sum, 0, na * sizeof(int)); for (i = 0; i < n ; i++) sum[ rank[i] = s[i]]++; for (i = 1; i < na; i++) sum[i] += sum[i - 1]; for (i = 0; i < n; i++) sa[ --sum[ rank[i] ] ] = i; //第一关键字排序 for (len = 1; len < n; len <<= 1) { //第二关键字排序：首先将sa看做第二关键字，j就是第一关键字，将第一关键字按照第二关键字//出现的顺序放入到tsa中，也就完成了双关键字的排序for (i = 0; i < n; i++) { j = sa[i] - len; if (j < 0) j += n; tsa[ sum[ rank[j] ]++ ] = j; } sa[ tsa[ sum[0] = 0 ] ] = j = 0; //j代表排名for (i = 1; i < n; i++) {//如果当前项的两个关键字中，第一个和第二个都和排在前面项的相同，则代表//该项和前一项排名相同，j不变if (rank[ tsa[i] ] != rank[ tsa[i-1] ] || rank[ tsa[i]+len ]!=rank[ tsa[i-1]+len ]) sum[++j] = i; //这里sum记录了相同项中，第一次出现的项的排名，也就是相同项中排名最高者 sa[ tsa[i] ] = j; //这里将rank值放入sa，sa这里只是作为一个缓冲} memcpy(rank, sa , n * sizeof(int)); memcpy(sa  , tsa, n * sizeof(int)); if (j >= n - 1) break; } }void lcp(){                          // O(4 * N) int i, j, k; j = rank[0];height[0] = 0;for (i = k = 0; i < n - 1; i++, k++) while (k >= 0 && s[i] != s[ sa[j-1] + k ]) height[j] = (k--), j = rank[ sa[j] + 1 ]; }/*利用二分法搜索可能的长度，然后判断是不是无重叠重复子串的长度。对于每一个长度d，将height中大于等于d的保留，小于的去除，这样就将height非为几组，因为求最大重复子串就是求某个区间中最小的height，因此最长的无重叠子串一定在同一组内，判断是否重叠就是，利用sa数组看它们在实际数组中的位置间隔是否比d大，如果da则不重叠。*/bool check(int mid){int low = sa[1], high = sa[1];for (int i = 2; i < n; ++i){if (height[i] >= mid){low = min(low, sa[i]);high = max(high, sa[i]);if (high - low >= mid) return true;}else{low = sa[i];high = sa[i];}}return false;}void solve(){int l = 1, r = n, mid, i, ans = 0;while (l < r){mid = (l + r) / 2;if (check(mid)) ans = mid, l = mid + 1;else r = mid;}if (ans >= 4) printf("%d\n", ans+1);else printf("0\n");}int main(){int i;while (scanf("%d", &n) && n){for (i = 0; i < n; ++i)scanf("%d", &s[i]);for (i = 0; i < n-1; ++i)s[i] = s[i] - s[i+1] + 100;s[i] = 0;makesa(200);lcp();solve();}return 0;}