后缀数组

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

#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;}