[USACO5.1]音楽をテーマにミュージカルのテーマ - 接尾辞配列

説明

繰り返されない部分文字列を重ね最長の文字列を探しています。

解決

二つの同一のサブ文字列がある場合、それは2つの接尾辞です\（LCP \）と同じ。

だから、二分法の解答考える\（K \を）の連続期間がある場合には、\（高さ\）未満ではありません（K \）\、その後、2つのサフィックス間の間隔のこの時期、\（LCPが\）未満ではありませんK（\ \）、レコード何間隔\（\最大\ {sa_i \ } \） と\（\分\ {sa_i \} \） 、もし\（\最大\ {sa_i \ } - \分\ {sa_i \} \ルK \）、2つの部分文字列が重複しないこと。

長い間の差として、それを解決することができるようまた、この問題は、問題の「再デプロイ」を持っています。しかしなお、2つの時間微分決意条件がなければならない後の\（\最大\ {sa_i \} - \分\ {sa_i \} <K \） 。

コード

#include <bits/stdc++.h>
using namespace std;

inline int ty() {
    char ch = getchar(); int x = 0, f = 1;
    while (ch < '0' || ch > '9') { if (ch == '-') f = -1; ch = getchar(); }
    while (ch >= '0' && ch <= '9') { x = x * 10 + ch - '0'; ch = getchar(); }
    return x * f;
}

const int _ = 5000 + 10;
const int INF = 0x3f3f3f3f;
int N, s[_], rnk[_], sa[_], height[_];

void SA() {
    static int t[_], a[_], buc[_], fir[_], sec[_], tmp[_];
    copy(s + 1, s + N + 1, t + 1);
    sort(t + 1, t + N + 1);
    int *end = unique(t + 1, t + N + 1);
    for (int i = 1; i <= N; ++i) a[i] = lower_bound(t + 1, end, s[i]) - t;
    for (int i = 1; i <= N; ++i) ++buc[a[i]];
    for (int i = 1; i <= N; ++i) buc[i] += buc[i - 1];
    for (int i = 1; i <= N; ++i) rnk[i] = buc[a[i] - 1] + 1;
    for (int len = 1; len <= N; len <<= 1) {
        for (int i = 1; i <= N; ++i) {
            fir[i] = rnk[i];
            sec[i] = i + len > N ? 0 : rnk[i + len];
        }
        fill(buc + 1, buc + N + 1, 0);
        for (int i = 1; i <= N; ++i) ++buc[sec[i]];
        for (int i = 1; i <= N; ++i) buc[i] += buc[i - 1];
        for (int i = 1; i <= N; ++i) tmp[N - --buc[sec[i]]] = i;
        fill(buc + 1, buc + N + 1, 0);
        for (int i = 1; i <= N; ++i) ++buc[fir[i]];
        for (int i = 1; i <= N; ++i) buc[i] += buc[i - 1];
        for (int i, j = 1; j <= N; ++j) {
            i = tmp[j];
            sa[buc[fir[i]]--] = i;
        }
        bool same = false;
        for (int i, j = 1, last = 0; j <= N; ++j) {
            i = sa[j];
            if (!last) rnk[i] = 1;
            else if (fir[i] == fir[last] && sec[i] == sec[last]) rnk[i] = rnk[last], same = true;
            else rnk[i] = rnk[last] + 1;
            last = i;
        }
        if (!same) break;
    }
    for (int i = 1, k = 0; i <= N; ++i) {
        if (rnk[i] == 1) k = 0;
        else {
            if (k > 0) --k;
            int j = sa[rnk[i] - 1];
            while (i + k <= N && j + k <= N && a[i + k] == a[j + k]) ++k;
        }
        height[rnk[i]] = k;
    }
}

bool check(int mid) {
    int minsa = INF, maxsa = -1;
    for (int i = 1; i <= N; ++i) {
        if (height[i] < mid) minsa = maxsa = sa[i];
        else {
            minsa = min(minsa, sa[i]);
            maxsa = max(maxsa, sa[i]);
            if (maxsa - minsa > mid) return true;
        }
    }
 return false;
}

int main() {
#ifndef ONLINE_JUDGE
    freopen("theme.in", "r", stdin);
    freopen("theme.out", "w", stdout);
#endif
    N = ty();
    for (int i = 1; i <= N; ++i) s[i] = ty();
    for (int i = 1; i <= N; ++i) s[i] = s[i + 1] - s[i] + 100;
    SA();
    // for (int i = 1; i <= N; ++i) printf("%d ", height[i]);
    int l = 0, r = N;
    while (l < r) {
        int mid = (l + r + 1) >> 1;
        if (check(mid)) l = mid;
        else r = mid - 1;
    }
    printf("%d\n", l >= 4 ? l + 1 : 0);
    return 0;
}