后缀数组(๑Ő௰Ő๑)ﾉ“

---

后缀数组与LCP

后缀数组

首先，将一个字符串的所有后缀按字典序排序。
后缀数组 $sa[i]$：表示所有后缀在排完序后，排名为 $i$ 的后缀在原串中的位置。
名次数组 $rank[i]$：表示所有后缀在排序完后，原字符串中第 $i$ 名现在的排名。
显然，$sa$ 与 $rank$ 是互逆的，即$$rank[sa[i]]=i,sa[rank[i]]=i;$$

倍增求sa

我们求后缀的排名，可以每次求出区间 $[i,i$+$2\ast w$−$1]$ 的排名，$w$=$1,2,4,8...$求 $[i,i$+$2\ast w$−$1]$ 的排名，我们可以将它分为 $[i,i$+$w$-$1]$ 和 $[i$+$w,i$+$2\ast w$−$1]$ 两部分。
所有可以以 $[i,i$+$w$-$1]$ 的排名为第一关键字，$[i$+$w,i$+$2\ast w$−$1]$的排名为第二关键字，来对所有的 $i$ 排序。
于是我们就可以用个 $sort$ 在 $O(nlo{g}^{2}n)$ 内进行后缀排序了。

CODE

int rnk[N],sa[N],h[N];
struct p{
   
	int a,b,id;
}p[N];
bool cmp(const P&a,const P&b){
   
	return a.a==b.a?a.b<b.b:a.a<b.a;
}
int init_sa(char *s){
   
	int n=std::strlen(s);
	for(int i=0;i<n;i++) rnk[i]=s[i];
	for(int i=1;i<n;i<<=1){
   
		for(int j=0;j<n;j++)
			p[j]=P{
   rnk[j],j+1>=n?-1:rnk[j+1],j};
		std::sort(p,p+n,cmp);
		int c;
		for(int j=0;j<n;j++){
   
			rnk[p[j].id]=c;
			c+=!(p[j].a==p[j+1].a&&p[j].b==p[j+1].b);
		}
	}
}

LCP(最长公共前缀)

LCP，即最长公共前缀，Longgest Common Prefix。
记 $lcp(i,j)$ 表示 $s{a}_i$ 与 $s{a}_j$ 的最长公共前缀。
则对任意