51nod 1732 婚姻介绍所 (后缀数组、RMQ）

题目

题解

学过后缀数组的话，读完题基本就可以直接写了。（送分题？）

先求出后缀数组和高度数组，然后ST预处理lcp值的范围最小值，最后O（1）查询即可。

debug：
这个送分（命）题，我写了两个神奇的bug。
第一个是每次倍增排序后确定排名时cmp内应该比较的是sa[i-1]和sa[i].因为你要比较下第i名与第i-1名是否排名相同，肯定是比较第i名的后缀而不是i后缀。
第二个bug就更奇葩了，查询区间最值时以2为底取对数，我写成了log（L）？？？

AC代码

//#include <bits/stdc++.h>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int maxn=1e3+7;

int k,sa[maxn],lcp[maxn],n,rk[maxn],tmp[maxn];
char s[maxn];
bool cmp(const int &i,const int &j)
{
    if(rk[i]!=rk[j]) return rk[i]<rk[j];
    int fi=i+k<n?rk[i+k]:-1;
    int fj=j+k<n?rk[j+k]:-1;
    return fi<fj;
}
void get_sa()
{
    for(int i=0;i<n;i++)
    {
        sa[i]=i;
        rk[i]=s[i];
    }
    for(k=1;k<=n;k*=2)
    {
        sort(sa,sa+n,cmp);
        tmp[sa[0]]=0;
        for(int i=1;i<n;i++)
            tmp[sa[i]]=tmp[sa[i-1]]+(cmp(sa[i-1],sa[i])?1:0);
        for(int i=0;i<n;i++) rk[i]=tmp[i];
    }
}
void get_lcp()
{
    for(int i=0;i<n;i++) rk[sa[i]]=i;

    int h=0;
    for(int i=0;i<n;i++)
    {
        if(rk[i]==n-1)
        {
            lcp[rk[i]]=0;
            h=0;
            continue;
        }
        int j=sa[rk[i]+1];
        if(h) h--;
        for(;i+h<n && j+h<n;h++)
            if(s[i+h]!=s[j+h]) break;
        lcp[rk[i]]=h;
    }
}
int dp[maxn][20];
void rmq_init()
{
    lcp[n-1]=0;
    memset(dp,0x3f,sizeof(dp));
    for(int i=0;i<n;i++) dp[i][0]=lcp[i];
    for(int j=1;(1<<j)<=n;j++)
        for(int i=0;i+(1<<j)<n;i++)
        dp[i][j]=min(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
}
int query(int x,int y)
{
    int L=y-x+1;
    int k=(int)log2(L);
    return min(dp[x][k],dp[y-(1<<k)+1][k]);
}

int main()
{
    int q;
    while(~scanf("%d",&n))
    {
        scanf("%s",s);
        scanf("%d",&q);
        get_sa();
        get_lcp();
        rmq_init();
        while(q--)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            if(rk[x]>rk[y])
            {
                int z=x;
                x=y;
                y=z;
            }
            if(x==y)
            {
                printf("%d\n",n-x);
                continue;
            }

            int ans=query(rk[x],rk[y]-1);
            printf("%d\n",ans);
        }
    }
    return 0;
}