【模板】後綴數組

模板

#include<bits/stdc++.h>
#define N 1000005
using namespace std;
int t1[N],t2[N],sa[N],h[N],rk[N],c[N],a[N],g[36][N],lg[N];
inline void SA(int n,int m)
{
	//DA->SA
    int *x=t1,*y=t2,p=0,f=0;
    for(int i=1;i<=m;i++) c[i]=0;
    for(int i=1;i<=n;i++) c[x[i]=a[i]]++;
    for(int i=1;i<=m;i++) c[i]+=c[i-1];
    for(int i=n;i;i--) sa[c[x[i]]--]=i;
    for(int i=1;i<=n&&p<=n;i<<=1)
	{
        p=0;
        for(int j=n-i+1;j<=n;j++) y[++p]=j;
        for(int j=1;j<=n;j++) if(sa[j]>i)y[++p]=sa[j]-i;
        for(int j=1;j<=m;j++) c[j]=0;
        for(int j=1;j<=n;j++) c[x[y[j]]]++;
        for(int j=1;j<=m;j++) c[j]+=c[j-1];
        for(int j=n;j;j--) sa[c[x[y[j]]]--]=y[j];
        swap(x,y);
		x[sa[1]]=1;
		p=2;
        for(int j=2;j<=n;j++)
			x[sa[j]]=y[sa[j]]==y[sa[j-1]]&&y[sa[j]+i]==y[sa[j-1]+i]?p-1:p++;
        m=p;
    }
	//Rank
    for(int i=1;i<=n;i++)
		rk[sa[i]]=i;
	//Height
    for(int i=1;i<=n;i++)
	{
        int j=sa[rk[i]-1];
        if(f) f--;
		while(a[i+f]==a[j+f])
			f++;
        h[rk[i]]=g[0][rk[i]]=f;
    }
	//ST
	for(int i=1;(1<<i)<=n;i++)
		for(int j=0;j+(1<<i)-1<=n;j++)
			g[i][j]=min(g[i-1][j],g[i-1][j+(1<<(i-1))]);
	for(int i=2;i<1000000;i++)
		lg[i]=lg[i>>1]+1;
}
inline int lcp(int x,int y)
{
	if(x>y)
		swap(x,y);
	int x1=lg[y-x];
	return min(g[x1][x+1],g[x1][y-(1<<x1)+1]);
}
int main()
{

}