HDU 5442 Favorite Donut

Favorite Donut

Time Limit: 1500/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 79    Accepted Submission(s): 17

Problem Description

Lulu has a sweet tooth. Her favorite food is ring donut. Everyday she buys a ring donut from the same bakery. A ring donut is consists of n parts. Every part has its own sugariness that can be expressed by a letter from a to z (from low to high), and a ring donut can be expressed by a string whose i-th character represents the sugariness of the i−th part in clockwise order. Note that z is the sweetest, and two parts are equally sweet if they have the same sugariness.

Once Lulu eats a part of the donut, she must continue to eat its uneaten adjacent part until all parts are eaten. Therefore, she has to eat either clockwise or counter-clockwise after her first bite, and there are 2n ways to eat the ring donut of n parts. For example, Lulu has 6 ways to eat a ring donut abc: abc,bca,cab,acb,bac,cba. Lulu likes eating the sweetest part first, so she actually prefer the way of the greatest lexicographic order. If there are two or more lexicographic maxima, then she will prefer the way whose starting part has the minimum index in clockwise order. If two ways start at the same part, then she will prefer eating the donut in clockwise order. Please compute the way to eat the donut she likes most.

Input
First line contain one integer T,T≤20, which means the number of test case.

For each test case, the first line contains one integer n,n≤20000, which represents how many parts the ring donut has. The next line contains a string consisted of n lowercase alphabets representing the ring donut.

Output
You should print one line for each test case, consisted of two integers, which represents the starting point (from 1 to n) and the direction (0 for clockwise and 1 for counterclockwise).

Sample Input

2

4

abab

4

aaab

 

Sample Output

2 0

4 0

 

Source

2015 ACM/ICPC Asia Regional Changchun Online

解题：后缀自动机

#include <bits/stdc++.h>
using namespace std;
const int maxn = 100010;
struct SAM {
    struct node {
        int son[26],f,len;
        void init() {
            f = -1;
            len = 0;
            memset(son,-1,sizeof son);
        }
    } sn[maxn<<1];
    int tot,last;
    void init() {
        tot = last = 0;
        sn[tot++].init();
    }
    int newnode() {
        sn[tot].init();
        return tot++;
    }
    void extend(int c) {
        int np = newnode(),p = last;
        sn[np].len = sn[last].len + 1;
        while(p != -1 && sn[p].son[c] == -1) {
            sn[p].son[c] = np;
            p = sn[p].f;
        }
        if(p == -1) sn[np].f = 0;
        else {
            int q = sn[p].son[c];
            if(sn[p].len + 1 == sn[q].len) sn[np].f = q;
            else {
                int nq = newnode();
                sn[nq] = sn[q];
                sn[nq].len = sn[p].len + 1;
                sn[np].f = sn[q].f = nq;
                while(p != -1 && sn[p].son[c] == q) {
                    sn[p].son[c] = nq;
                    p = sn[p].f;
                }
            }
        }
        last = np;
    }
} sam;
char str[maxn];
int fail[maxn];
void getFail(const char *word) {
    fail[0] = -1;
    fail[1] = 0;
    for(int i = 0,j = -1; word[i]; ++i) {
        while(j != -1 && word[i] != word[j]) j = fail[j];
        fail[i+1] = ++j;
    }
}
char text[maxn];
int mymin(int x,int y) {
    return x > y?y:x;
}
int mymax(int x,int y) {
    return x > y?x:y;
}
int FindIndex(const char *word,int (*f)(int,int),int idx,int len) {
    for(int i = 0, j = 0; text[i] ; ++i) {
        while(j > -1 && word[j] != text[i]) j = fail[j];
        if(!word[++j]) {
            if(i - len + 2 <= len)idx = f(idx,i - len + 2);
            j = fail[j];
        }
    }
    return idx;
}
int main() {
    int kase,len;
    scanf("%d",&kase);
    while(kase--) {
        sam.init();
        scanf("%d",&len);
        scanf("%s",str);
        int p = 0;
        for(int i = 0; i < (len<<1); ++i)
            sam.extend(str[i%len]-'a');
        string ss = "",sb = "";
        for(int i = 0; i < len; ++i) {
            for(int j = 25; j >= 0; --j) {
                if(sam.sn[p].son[j] != -1) {
                    p = sam.sn[p].son[j];
                    ss += char('a' + j);
                    break;
                }
            }
        }
        int clockwise = sam.sn[p].len - len + 1;
        getFail(ss.c_str());
        strcpy(text,str);
        strcpy(text + len,str);
        clockwise = FindIndex(ss.c_str(),mymin,len*2,len);
        reverse(str,str + len);
        sam.init();
        for(int i = p = 0; i < (len<<1); ++i) {
            sam.extend(str[i%len]-'a');
        }
        for(int i = 0; i < len; ++i) {
            for(int j = 25; j >= 0; --j) {
                if(sam.sn[p].son[j] != -1) {
                    p = sam.sn[p].son[j];
                    sb += char('a' + j);
                    break;
                }
            }
        }
        getFail(sb.c_str());
        strcpy(text,str);
        strcpy(text + len,str);
        int anticlockwise = len - FindIndex(sb.c_str(),mymax,0,len) + 1;
        if(ss == sb) {
            if(clockwise == anticlockwise) {
                printf("%d %d
",clockwise,0);
            } else if(clockwise < anticlockwise) {
                printf("%d %d
",clockwise,0);
            } else printf("%d %d
",anticlockwise,1);
        } else if(ss < sb) printf("%d %d
",anticlockwise,1);
        else printf("%d %d
",clockwise,0);
    }
    return 0;
}