后缀数组

后缀数组是处理字符串的有力工具。后缀数组是后缀树的一个非常精巧的替代品，它比后缀树容易编程实现，能够实现后缀树的很多功能而时间复杂度也并不逊色，而且它比后缀树所占用的内存空间小很多。可以说，在信息学竞赛中后缀数组比后缀树要更为实用。（以上内容摘自度娘~~~）
不久前写过后缀数组，今天来复习一下……

//bzoj 1031
#include <cstdio>
#include <algorithm>
#include <iostream>
#include <string>
#define maxn 200005
using namespace std;
string ans,s;
int n,c[maxn],sa[maxn],t2[maxn],t1[maxn],h[maxn],rank[maxn];
int *x=t1,*y=t2;

void get_sa()
{
    int m=257;
    for (int i=0;i<m;i++) c[i]=0;
    for (int i=0;i<n;i++) c[x[i]=sa[i]]++;
    for (int i=1;i<m;i++) c[i]+=c[i-1];
    for (int i=n-1;i>=0;i--) sa[--c[x[i]]]=i;
    for (int k=1;k<=n;k<<=1)
    {
        int p=0;
        for (int i=n-k;i<n;i++) y[p++]=i;//第二关键字排序 
        for (int i=0;i<n;i++) if (sa[i]>=k) y[p++]=sa[i]-k;//sa[i]>=k则有第二关键字 ，否则没有y 
        for (int i=0;i<m;i++) c[i]=0;//第一关键字排序
        for (int i=0;i<n;i++) c[x[y[i]]]++;
        for (int i=1;i<m;i++) c[i]+=c[i-1];
        for (int i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
        swap(x,y);
        p=1;x[sa[0]]=0;
        for (int i=1;i<n;i++)
          x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
        if (p>=n) break;
        m=p;
    }
}

void get_h()
{
    int k=0,j;
    for (int i=0;i<n;i++) rank[sa[i]]=i;
    for (int i=0;i<n;h[rank[i++]]=k)
      for (k?k--:0,j=sa[rank[i]-1];s[i+k]==s[j+k];k++);
}

int main()
{
    cin>>s;
    s=s+s;
    n=s.length();
    for (int i=0;i<n;i++) sa[i]=s[i];
    get_sa();
    get_h(); 
    //for (int i=0;i<n;i++) printf("%d ",sa[i]);
    for (int i=0;i<n;i++)
      if (sa[i]<n/2) ans+=s[sa[i]+n/2-1];
    cout<<ans;
    return 0;
}