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 = 10000 + 5;

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

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

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

struct SuffixAutomaton{
  int sz, last;

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

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

    if(p == -1) st[cur].pre = 0;
    else{
      int q = st[p].next[c];
      if(st[q].len == st[p].len + 1) st[cur].pre = q;
      else{
        int cle = sz ++;
        st[cle].len = st[p].len + 1;
        st[cle].pre = st[q].pre;
        for(int i = 0; i < 26; ++ i) st[cle].next[i] = st[q].next[i];
        for(; p != -1 && 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
", ans);
}

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

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

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

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