后缀数组

解释及作用

对于一个字符串，将其所有后缀排序，用后缀中的第一个字符代表这个后缀，那么排序后的数组就是后缀数组，例如，对于字符串ababs，它的后缀数组是1 3 2 4 5，因为
1.ababs–开始位置为1，排名第1
2.babs—开始位置为2，排名第3
3.abs—–开始位置为3，排名第2
4.bs——开始位置为4，排名第4
5.s——–开始位置为5，排名第5
对于这个性质,我们可以用来利用二分进行字符串匹配,复杂度为(匹配串长度*log被匹配串长度),因而相比与kmp,它在被匹配串长度较长时有优势。

实现方法

如果直接枚举后缀进行排序，那么复杂度为O(n^2*logn),利用倍增的方法可以优化此算法，那么复杂度为O(n *(logn)^2).
具体实现就是知道长度为l时，以每个字符为首的长度为l的字符串的排名，就可以推出以各个字符为首的长度为2*l的字符串的排名，比较方法就是先比较前面的排名，再比较后面的。

代码

#include<iostream> 
#include<cstdio>
#include<algorithm>
#define N 1001000
using namespace std;

int len,rank[N],st[N],tmp[N],n;
string str;

inline bool cmp(int u,int v)
{
    if(rank[u]!=rank[v]) return rank[u]<rank[v];
    int nu,nv;
    nu= (u+len>=n)?-1:rank[u+len];
    nv= (v+len>=n)?-1:rank[v+len];
    return nu<nv;
}

int main()
{
    int i,j;
    cin>>str;
    n=str.size();
    for(i=0;i<n;i++) st[i]=i,rank[i]=str[i];
    for(len=1;len<=n;len<<=1)
    {
        sort(st,st+n,cmp);
        tmp[st[0]]=1;
        for(i=1;i<n;i++)
        {
            tmp[st[i]]=tmp[st[i-1]]+cmp(st[i-1],st[i]);
        }
        for(i=0;i<n;i++) rank[i]=tmp[i];
    }
    for(i=0;i<n;i++) printf("%d ",st[i]+1);
}