【后缀数组】后缀排序

最近几天有很多题目没有写，尽量补上吧。

hz2016评测《《点击访问

caioj《《点击访问

后缀数组》后缀自动机

树状数组》线段树

其实两种都是用一个特殊的简单结构，维护一部分高级算法的查询结构。后缀自动机顾名思义，就是处理后缀的啦，这道题是模板题，把数组的所有后缀按照大到小排序。具体还是看代码注释吧，结合代码会比较好理解。

#include<map>
#include<queue>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Maxchar 1524288 //1000000 +2^19
#define Maxs 30
#define mes(x,y) memset(x,y,sizeof(x));
#define mpy(x,y) memcpy(x,y,sizeof(x))
#define INF 2147483647
using namespace std;
int a[Maxchar+1],tt[Maxchar+1];
char s[Maxchar+1];
int Rank[Maxchar+1],sa1[Maxchar+1],Rank2[Maxchar+1],sa2[Maxchar+1],Rsort[Maxchar+1];
//Rank名次数组: 你排第几？sa后缀数组：排第几的是谁？(记录是的起始位置) 后面的数字表示第几关键字
void get_sa(int n,int m){
memcpy(Rank,a,sizeof(Rank));
//预处理第一关键字的排名：对长度为1的字符串进行基数排序，求出sa
memset(Rsort,0,sizeof(Rsort));
for(int i=1;i<=n;i++)Rsort[Rank[i]]++;
for(int i=1;i<=m;i++)Rsort[i]+=Rsort[i-1];
for(int i=n;i>=1;i--)sa1[Rsort[Rank[i]]--]=i;
int ln=1,p=0;//ln为当前处理好的子串的长度,p表示有多少不相同的子串
while(p<n){//现在处理的是长度为ln*2的子串，开始将两个关键字合并排序
//为啥不写ln<=n？因为有可能不需要排到ln=n就排完了，而当p==n就说明全部字符串不同，后面再也不会改变
int k=0;//处理第二关键字的排名：
for(int i=n-ln+1;i<=n;i++)sa2[++k]=i;//先把要补0的后缀处理好，他们肯定在前面
for(int i=1;i<=n;i++)if(sa1[i]-ln>0)sa2[++k]=sa1[i]-ln;
//sa1[i]是sa1[i]-ln的第二关键字，因为我们枚举的是第二关键字的排名
//sa2记录sa1[i]-ln的第二关键字的排名，指向sa1[i]-ln
memset(Rsort,0,sizeof(Rsort));
for(int i=1;i<=n;i++)Rsort[Rank[i]]++;
for(int i=1;i<=m;i++)Rsort[i]+=Rsort[i-1];
for(int i=n;i>=1;i--)sa1[Rsort[Rank[sa2[i]]]--]=sa2[i];
for(int i=1;i<=n;i++)tt[i]=Rank[i];//由于Rank会改变，判断不了一、二关键字是否相等
//用sa得到新的Rank数组，为什么预处理的时候不能求出Rank？因为原来的Rank没有两个关键字
p=1;Rank[sa1[1]]=1;
for(int i=2;i<=n;i++){
if(tt[sa1[i]]!=tt[sa1[i-1]]||tt[sa1[i]+ln]!=tt[sa1[i-1]+ln] )p++;
//一二关键字其中一个不同，前后的字符串不同
Rank[sa1[i]]=p;
}
m=p;ln*=2;
}
}
int main(){
scanf("%s",s+1);
int len=strlen(s+1);
for(int i=1;i<=len;i++)a[i]=s[i]-'a'+1;
get_sa(len,Maxs+1);
for(int i=1;i<len;i++)printf("%d ",sa1[i]);printf("%d\n",sa1[len]);
return 0;
}