后缀数组自学

sa[i]，表示在所有后缀中，排在第i位的后缀是哪个后缀。
rank[i]，表示以第i个字母为结尾的后缀在所有后缀中排第rank[i]个。
height 数组，表示height[i]=(suffix(SA[i-1])和suffix(SA[i])的最长公共前缀)，也就是排名相邻的两个后缀的最长公共前缀的长度 。

#include<iostream>
#include <cstdio>
#include <iostream>
#include <cstdlib>
#include <algorithm>
#include <cstring>
using namespace std;
const int maxn=200010;

int wa[maxn],wb[maxn],wsf[maxn],wv[maxn];
int rank[maxn],height[maxn],s[maxn],sa[maxn];
void make_sa(int n,int m){
    int i,j,p,*x=wa,*y=wb;
    for(i=0;i<m;i++)  wsf[i]=0;
    for(i=0;i<n;i++)  wsf[x[i]=s[i]]++;
    for(i=1;i<m;i++)  wsf[i]+=wsf[i-1];
    for(i=n-1;i>=0;i--)  sa[--wsf[x[i]]]=i;
    for(p=1,j=1;p<n;j<<=1,m=p){
        for(p=0,i=n-j;i<n;i++)  y[p++]=i;
        for(i=0;i<n;i++) {
            if(sa[i]>=j)  y[p++]=sa[i]-j;
        }
        for(i=0;i<n;i++)  wv[i]=x[y[i]];
        for(i=0;i<m;i++)  wsf[i]=0;
        for(i=0;i<n;i++)  wsf[wv[i]]++;
        for(i=1;i<m;i++)  wsf[i]+=wsf[i-1];
        for(i=n-1;i>=0;i--)  sa[--wsf[wv[i]]]=y[i];
        swap (x,y);
        x[sa[0]]=0;
        for(p=1,i=1;i<n;i++){
            x[sa[i]]=(y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+j]==y[sa[i]+j]) ? p-1 : p++;
        }
    }
}
void make_height(int n){
    int i,j,k=0;
    for(i=1;i<=n;i++) rank[sa[i]]=i;
    for(i=0;i<n;i++){
        if(k)k--;
        else k=0;
        j=sa[rank[i]-1];
        while(s[i+k]==s[j+k]) k++;
        height[rank[i]]=k;
    }
}
int main (){
    make_sa (srelen (s)+1,s中不同字符种类);
    make_height (strlen (s));
}