UVA 719 / POJ 1509 Glass Beads (最小表示法/后缀自动机)

题目大意:

给出一个长度为N的字符串，求其字典序最小的循环同构。

N<=10W。

算法讨论：

算法一、最小表示法。定义题。

算法二、后缀自动机。

Codes:

 #include <iostream>
 #include <algorithm>
 #include <cstdlib>
 #include <cstdio>
 #include <cstring>
 #include <cassert>
 using namespace std;
 const int L =  + ;

 int ans = ;
 char str[L<<];

 struct State{
   int len, pre;
   int next[];

   State(){
     len = pre = ;
     memset(next, , sizeof next);
   }
   void clear(){
     len = pre = ;
     memset(next, , sizeof next);
   }
 }st[L<<];

 struct SuffixAutomaton{
   int sz, last;

   void Init(){
     last = ;
     sz = ;
     for(int i = ; i < (L<<); ++ i)
       st[i].clear();
     st[].len = ; st[].pre = -;
     sz ++;
   }
   void Extend(int c){
     int cur = sz ++;
     st[cur].len = st[last].len + ;
     int p;

     for(p = last; p != - && !st[p].next[c]; p = st[p].pre)
       st[p].next[c] = cur;

     if(p == -) st[cur].pre = ;
     else{
       int q = st[p].next[c];
       if(st[q].len == st[p].len + ) st[cur].pre = q;
       else{
         int cle = sz ++;
         st[cle].len = st[p].len + ;
         st[cle].pre = st[q].pre;
         for(int i = ; i < ; ++ i) st[cle].next[i] = st[q].next[i];
         for(; p != - && st[p].next[c] == q; p = st[p].pre)
           st[p].next[c] = cle;
         st[q].pre = st[cur].pre = cle;
       }
     }
     last = cur;
   }
 }SAM;

 void Input(){
   scanf("%s", str);
 }

 void Output(){
   printf("%d\n", ans);
 }

 void Solve(){
   int len = strlen(str);
   SAM.Init();
   for(int i = ; i < len; ++ i)
     str[i + len] = str[i];
   for(int i = ; i < (len<<); ++ i){
     SAM.Extend(str[i] - 'a');
   }

   int p = ;
   for(int i = ; i < len; ++ i){
     for(int j = ; j < ; ++ j){
       if(st[p].next[j]){
         p = st[p].next[j];
         break;
       }
     }
   }

   ans = st[p].len - len + ;
 }

 int main(){
   int t;
   scanf("%d", &t);
   while(t --){
     Input();
     Solve();
     Output();
   }
   return ;
 }

后缀自动机