求解后缀数组的最长前缀

用RMQ中的ST算法求出Height数组的区间最小值，然后进行区间查询，先计算出查询值的

rank值，然后再查询。

代码：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#define  LL long long
#define  ULL unsigned long long
using namespace std;
const int MAXN=100010;
int wa[MAXN],wb[MAXN],wv[MAXN],Ws[MAXN];
int Min[MAXN][25];
int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b]&&r[a+l]==r[b+l];
}
void da(const char r[],int sa[],int n,int m)
{
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0; i<m; i++)
        Ws[i]=0;
    for(i=0; i<n; i++)
        Ws[x[i]=r[i]]++;
    for(i=1; i<m; i++)
        Ws[i]+=Ws[i-1];
    for(i=n-1; i>=0; i--)
        sa[--Ws[x[i]]]=i;
    for(j=1,p=1; p<n; j*=2,m=p)
    {
        for(p=0,i=n-j; i<n; i++)
            y[p++]=i;
        for(i=0; i<n; i++)
            if(sa[i]>=j)
                y[p++]=sa[i]-j;
        for(i=0; i<n; i++)
            wv[i]=x[y[i]];
        for(i=0; i<m; i++)
            Ws[i]=0;
        for(i=0; i<n; i++)
            Ws[wv[i]]++;
        for(i=1; i<m; i++)
            Ws[i]+=Ws[i-1];
        for(i=n-1; i>=0; i--)
            sa[--Ws[wv[i]]]=y[i];
        for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1; i<n; i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
    }
    return;
}
int sa[MAXN],Rank[MAXN],height[MAXN];
void calheight(const char *r,int *sa,int n)
{
    int i,j,k=0;
    for(i=1; i<=n; i++)
        Rank[sa[i]]=i;
    for(i=0; i<n; height[Rank[i++]]=k)
        for(k?k--:0,j=sa[Rank[i]-1]; r[i+k]==r[j+k]; k++);
    for(int i=n; i>=1; --i)
        ++sa[i],Rank[i]=Rank[i-1];
}
void RMQ(int n)
{
    for(int i=1;i<=n;i++)
        Min[i][0]=height[i];
        Min[1][0]=MAXN;
    for(int i=1;(1<<i)<=n;i++)
    {
        for(int j=1;j+(1<<i)-1<=n;j++)
        {
            Min[j][i]=min(Min[j][i-1],Min[j+(1<<(i-1))][i-1]);
        }
    }
}
int Query(int l,int r)
{
    l=Rank[l];
    r=Rank[r];
    if(l>r)
        swap(l,r);
    l=l+1;
    int k=(int)(log((double)(r-l+1))/(log((double)2)));
    return min(Min[l][k],Min[r-(1<<k)+1][k]);
}
char str[MAXN];
int main()
{
    while(scanf("%s",str)!=EOF)
    {
        int len=strlen(str);
        int m;
        da(str,sa,len+1,130);
        calheight(str,sa,len);
        RMQ(len);
        scanf("%d",&m);
        int l,r;
        while(m--)
        {
            scanf("%d%d",&l,&r);
            int ans=Query(l,r);
            printf("%d\n",ans);
        }
    }
    return 0;
}