【模板】后缀数组

串处理问题神器，RMQ预处理后能O(1)回答俩后缀的LCP。

另外有一种常见处理方法：二分LCP长度H，将sa分块，则每一块内的后缀同且仅同当前块内的另一后缀LCP>=H。

还有一种常见思考方法：考虑贡献值。

Code：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int SN=100010;

int n,s[SN],sa[SN],c[SN],rank[SN],hei[SN],h[SN],lg[SN],rmq[SN][30],ta0[SN],ta1[SN];

bool Cmp(int i, int j) { return (i<=n && j<=n && s[i]==s[j]); }
void SufArr()
{
     int i,j,k,p,m=SN, *x=ta0,*y=ta1;
     for(i=1; i<=m; i++) c[i]=0;
     for(i=1; i<=n; i++) c[x[i]=s[i]]++;
     for(i=1; i<=m; i++) c[i]+=c[i-1];
     for(i=n; i>=1; i--) sa[c[x[i]]--]=i;
     for(k=1,p=0; ; k<<=1,p=0)
     {
	  for(i=n-k+1; i<=n; i++) y[++p]=i;
	  for(i=1; i<=n; i++) if(sa[i]-k>0) y[++p]=sa[i]-k;
	  for(i=1; i<=m; i++) c[i]=0;
	  for(i=1; i<=n; i++) c[x[y[i]]]++;
	  for(i=1; i<=m; i++) c[i]+=c[i-1];
	  for(i=n; i>=1; i--) sa[c[x[y[i]]]--]=y[i];
	  swap(x,y);
	  for(i=2,x[sa[1]]=p=1; i<=n; i++)
	       x[sa[i]]=(y[sa[i]]==y[sa[i-1]] && y[sa[i]+k]==y[sa[i-1]+k])? p: ++p;
	  if(p==n) break;  else m=p;
     }
     for(i=1; i<=n; i++) rank[sa[i]]=i;

     for(i=1,j=0; i<=n; i++) lg[i]=( (1<<j) >= (i+1)/2 )? j: ++j;
     for(i=1,h[0]=1; i<=n; i++)
	  if(rank[i]==1) h[i]=hei[rank[i]]=0;
	  else { for(k=Max(h[i-1],1); Cmp(i+k-1, sa[rank[i]-1]+k-1); k++); h[i]=hei[rank[i]]=k-1; }
}
void RMQ()
{
     int i,k,up=Lg(SN);
     for(i=1; i<=n; i++) rmq[i][0]=hei[i];
     for(k=1; k<=up; k++)
	  for(i=1; i+(1<<k)-1<=n; i++)
	       rmq[i][k]=Min(rmq[i][k-1], rmq[i+(1<<(k-1))][k-1]);
}
int Ask(int l, int r)
{
     int len,tmp;
     l=rank[l], r=rank[r]; if(l>r) tmp=l,l=r,r=tmp; l++;
     len=r+1-l, tmp=lg[len]-1, tmp=(tmp)?tmp:1;
     return Min(rmq[l][tmp], rmq[r-(1<<tmp)+1][tmp]);
}

int main()
{
     freopen("oj.in", "r", stdin);
     freopen("oj.out", "w", stdout);
     return 0;
}