后缀数组的倍增算法

一、基数排序。

题目描述：给你n个字符串，把它们排序。  n<=100000  （不要用sort水过，除非你在NOI的考场上）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<ctime>
#include<algorithm>
#define MAXN 100010
using namespace std;
int n,rank[MAXN],p[MAXN],sa[MAXN],cnt[MAXN];
char s[MAXN][15];
inline int read()
{
    int x=0,f=1;  char ch=getchar();
    while(!isdigit(ch))  {if(ch=='-')  f=-1;  ch=getchar();}
    while(isdigit(ch))  {x=x*10+ch-'0';  ch=getchar();}
    return x*f;
}
void Init()
{
    n=read();
    for(int i=1;i<=n;i++)
    {
        char ch[15];
        scanf("%s",ch);
        for(int j=1;j<=10;j++)
            s[i][j]=j>strlen(ch)?'a'-1:ch[j-1];
    }
}
void Radix_sort()
{
    for(int i=1;i<=n;i++)  sa[i]=i;
    for(int l=10;l>0;l--)
    {
        for(int i=1;i<=n;i++)  rank[i]=s[i][l]-'a'+1,p[i]=sa[i];
        for(int i=0;i<=26;i++)  cnt[i]=0;
        for(int i=1;i<=n;i++)  cnt[rank[i]]++;
        for(int i=1;i<=26;i++)  cnt[i]+=cnt[i-1];
        for(int i=n;i>=1;i--)  sa[cnt[rank[p[i]]]--]=p[i];
    }
}
void print()
{
    for(int i=1;i<=n;i++)
        printf("%d ",sa[i]);
}
int main()
{
    Init();
    Radix_sort();
    print();
    return 0;
}

二、求后缀数组。

int N,ans,rank[MAXN],sa[MAXN],p[MAXN],cnt[MAXN],tmp[MAXN],height[MAXN];
char s[MAXN];
bool equ(int x,int y,int l) {return rank[x]==rank[y]&&rank[x+l]==rank[y+l];}
void doubling()
{
    for(int i=1;i<=N;i++)  rank[i]=s[i],sa[i]=i;
    for(int pos=0,sig=255,l=0,i;pos<N;sig=pos)
    {
        for(i=N-l+1,pos=0;i<=N;i++)  p[++pos]=i;
        for(i=1;i<=N;i++)  if(sa[i]>l)  p[++pos]=sa[i]-l;
        for(i=1;i<=sig;i++)  cnt[i]=0;
        for(i=1;i<=N;i++)  cnt[rank[i]]++;
        for(i=1;i<=sig;i++)  cnt[i]+=cnt[i-1];
        for(i=N;i;i--)  sa[cnt[rank[p[i]]]--]=p[i];
        for(pos=0,i=1;i<=N;i++)  
            tmp[sa[i]]=equ(sa[i],sa[i-1],l)?pos:++pos;
        for(int i=1;i<=N;i++)  rank[i]=tmp[i];
        l=!l?1:l<<1;
    }
}

三、height数组

void get_height()
{
    for(int i=1,j=0,k;i<=N;i++)
    {
        if(!(k=sa[rank[i]-1]))  {j=0;  continue;}
        if(j)  j--;
        while(s[i+j]==s[k+j])  j++;
        height[rank[i]]=j;
        ans=max(ans,j);
    }
}