POJ 1509 后缀数组

链接：

http://poj.org/problem?id=1509

题意：

给你一个环形字符串，问从哪个地方截断使得得到的字符串字典序最小

题解：

计算s+s的后缀数组，然后遍历sa数组，找到第一个小于n的地方 输出即可

但是直接这样做会wa掉，因为有多个结果时，题目要求输出最小的，看了别人的博客知道在后面再加一个char('z'+1)就行了

代码：

int n, k;
int Rank[MAXN], tmp[MAXN];
int sa[MAXN], lcp[MAXN];

bool compare_sa(int i, int j) {
    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 construct_sa(string S, int *sa) {
    n = S.length();
    for (int i = 0; i <= n; i++) {
        sa[i] = i;
        Rank[i] = i < n ? S[i] : -1;
    }
    for (k = 1; k <= n; k *= 2) {
        sort(sa, sa + n + 1, compare_sa);
        tmp[sa[0]] = 0;
        for (int i = 1; i <= n; i++) 
            tmp[sa[i]] = tmp[sa[i - 1]] + (compare_sa(sa[i - 1], sa[i]) ? 1 : 0);
        for (int i = 0; i <= n; i++) Rank[i] = tmp[i];
    }
}

void construct_lcp(string S, int *sa, int *lcp) {
    int n = S.length();
    for (int i = 0; i <= n; i++) Rank[sa[i]] = i;
    int h = 0;
    lcp[0] = 0;
    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 (S[j + h] != S[i + h]) break;
        lcp[Rank[i] - 1] = h;
    }
}

int main() {
    ios::sync_with_stdio(false), cin.tie(0);
    int T;
    cin >> T;
    while (T--) {
        string s;
        cin >> s;
        s += s + char('z' + 1);
        construct_sa(s, sa);
        rep(i, 0, n + 1) if (sa[i] < n / 2) {
            cout << sa[i] + 1 << endl;
            break;
        }
    }
    return 0;
}