后缀数组

一点关于后缀数组的
乱七八糟的理解

#include<cstdio>
#include<cstring>
using namespace std;

const int maxn = 100010;

char r[maxn];
int wa[maxn];
int wb[maxn];
int wv[maxn];
int ws[maxn];//桶 
int S[maxn];

int cmp(int *r,int a,int b,int l){
    return r[a]==r[b]&&r[a+1]==r[b+1];
    //规定r[n-1]=0的好处
    //如果r[a]==r[b] 那么以a为起点和以b为起点的字符串肯定不会包含r[n-1] （因为只有这一个的值为0）
    //所以不用做特判，不用担心数组越界 
}

void da(int *sa,int n,int m){//r数组存放待排序的字符串
//sa为后缀数组 n为字符串长度 m为排名序号的范围 
    int i,j,p,*x=wa,*y=wb,*t;
    //???为了函数操作的方便，约定除r[n-1]外所有的r[i]都大于0 r[n-1]=0;  
    for(int i=0;i<m;i++) ws[i]=0;//清零 我大概只能看懂这个了233333
    //使用 基数排序 对长度为1的字符串进行排序
    for(int i=0;i<n;i++) ws[x[i]=r[i]]++;//x数组保存rank值
    for(int i=1;i<m;i++) ws[i]+=ws[i-1];
    for(int i=n-1;i>=0;i--) sa[--ws[x[i]]]=i;
    for(int i=0;i<n;i++) printf("%d ",sa[i]);
    //sa[i]表示排名为i的串首字母在原串中的位置
    //因为还没有排完序，所以当前x[i]表示的还是i位置所引领的串的首字母（s[i]），所以i位置所引领的串的排名自然就是“--ws[x[i]]”
    for(j=1,p=1;p<n;j*=2,m=p){//利用倍增思想继续完成
    //j表示的是当前要比较的字符串的长度 
        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;//y保存对第二关键字排序的结果 
        //小优化 
        //基数排序的结果可利用上一次求得的sa直接算出 
        for(i=0;i<n;i++) wv[i]=x[y[i]];
        for(i=0;i<m;i++) ws[i]=0;
        for(i=0;i<n;i++) ws[wv[i]]++;
        for(i=1;i<m;i++) ws[i]+=ws[i-1];
        for(i=n-1;i>=0;i--) sa[--ws[wv[i]]]=y[i];
        //对第一关键字进行排序 
        for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
           x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
        //由于rank值可能相同，所以需要比较两个字符串是否完全相同
        //【节省空间】用y数组保存rank值
        //小优化：复制整个数组可以用交换指针的值代替

    }
    return ;
}

int main(){
    scanf("%s",r);
    int len=strlen(r);
    da(S,len,127);
    return 0;
}