[模板]后缀数组

int s[MAXN+10];//char
int sa[MAXN+10],t[MAXN+10],t2[MAXN+10],c[MAXN+10];
void build_sa(int m)
{
int i,*x=t,*y=t2;
for(i=0;i<m;i++) c[i]=0;
for(i=0;i<n;i++) c[x[i]=s[i]]++;
for(i=1;i<m;i++) c[i]+=c[i-1];
for(i=n-1;i>=0;i--) sa[--c[x[i]]]=i;

for(int k=1;k<=n;k<<=1)
{
int p=0;
for(i=n-k;i<n;i++) y[p++]=i;
for(i=0;i<n;i++) if(sa[i]>=k) y[p++]=sa[i]-k;

for(i=0;i<m;i++) c[i]=0;
for(i=0;i<n;i++) c[x[y[i]]]++;
for(i=1;i<m;i++) c[i]+=c[i-1];
for(i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];

swap(x,y);
p=1;
x[sa[0]]=0;
for(i=1;i<n;i++)
if(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k])
x[sa[i]]=p-1;
else x[sa[i]]=p++;
if(p>=n) break;
m=p;
}
}
int Rank[MAXN+10],h[MAXN+10];
void getheight()
{
int i,k=0;
for(i=0;i<n;i++) Rank[sa[i]]=i;
for(i=0;i<n;i++)
{
if(k) k--;
if(!Rank[i]) continue;
int j=sa[Rank[i]-1];
while(s[i+k]==s[j+k]) k++;
h[Rank[i]]=k;
}
}

int LOG(int x)
{
int res=0;
while(x)
{
x>>=1;
res++;
}
return res-1;
}
int sparse_table[MAXN+10][MAXLOG];
void pre()
{
int i,j;
for(i=0;i<n;i++)
sparse_table[i][0]=h[i];
for(j=1;j<=MAXLOG;j++)
{
int up=n-(1<<j);
for(i=0;i<up;i++)
sparse_table[i][j]=min(sparse_table[i][j-1],sparse_table[i+(1<<(j-1))][j-1]);
}
}
int LCP(int x,int y)
{
x=Rank[x];
y=Rank[y];
Log=LOG(y-x+1);
return min(sparse_table[x][Log],sparse_table[y-(1<<Log)+1][Log]);
}