后缀数组总结及模板

http://www.cnblogs.com/zwfymqz/p/8413523.html

https://xminh.github.io/2018/02/27/后缀数组-最详细(maybe)讲解.html

 

 

【模板】后缀排序

题目背景

这是一道模板题。

题目描述

读入一个长度为 n 的由大小写英文字母或数字组成的字符串，请把这个字符串的所有非空后缀按字典序从小到大排序，然后按顺序输出后缀的第一个字符在原串中的位置。位置编号为 1 到 n。

输入输出格式

输入格式：

一行一个长度为 n 的仅包含大小写英文字母或数字的字符串。

输出格式：

一行，共n个整数，表示答案。

输入输出样例

输入样例#1： 

ababa

输出样例#1： 

5 3 1 4 2

说明

n <= 10^6

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=1e6+15;
char s[N]; int n,m=75,rank[N],sa[N],tp[N],tax[N],h[N];
 
inline void qsort()
 {for(int i=1;i<=m;i++) tax[i]=0;
  for(int i=1;i<=n;i++) tax[rank[i]]++;
  for(int i=1;i<=m;i++) tax[i]+=tax[i-1];
  for(int i=n;i>=1;i--) sa[tax[rank[tp[i]]]--]=tp[i];
 }
 
inline void suffix()
 {for(int i=1;i<=n;i++) rank[i]=s[i]-'0'+1,tp[i]=i;
  qsort();
  for(int w=1,p=0;p<n;m=p,w<<=1)
   {p=0;
    for(int i=1;i<=w;i++) tp[++p]=n-w+i;
    for(int i=1;i<=n;i++) if(sa[i]>w) tp[++p]=sa[i]-w;
    
    qsort();
    swap(tp,rank);
    rank[sa[1]]=p=1;
    for(int i=2;i<=n;i++)
     rank[sa[i]]=(tp[sa[i]]==tp[sa[i-1]]&&tp[sa[i]+w]==tp[sa[i-1]+w])?p:++p;
    
   }
 } 
inline void geth()
 {int j,k=0;
  for(int i=1;i<=n;i++)
   {if(k) k--;
    j=sa[rank[i]-1];
    while(s[i+k]==s[j+k])k++;
    h[rank[i]]=k;
   }
 } 
int main()
{		
 scanf("%s",s+1); n=strlen(s+1);	
	
 suffix();	
 
 for(int i=1;i<n;i++) printf("%d ",sa[i]);	printf("%d",sa[n]);
	
return 0;
}