LCP 模板

LCP

Description

给定串 \(S\) .

\(m\) 组询问 \((X, Y, L, R)\): 求 \(S[X,Y]\) 与 \(S[L,R]\) 的最长公共前缀.

Input

第一行一个串 \(S\).

第二行一个数 \(m\).

接下来 \(m\) 行, 每行 \(4\) 个数 \(X, Y, L, R\).

Output

输出共 \(m\) 行, 第 \(i\) 行为第 \(i\) 个询问的答案.

HINT

\(1 <= |S| <= 100000\), 且 \(S\) 由小写字母构成.

\(1 <= m <= 100000\)

\(1 <= L <= R <= |S|\)

\(1 <= X <= Y <= |S|\)

---

LCP标准模板，放一下...

Code:

#include <cstdio>
#include <cstring>
#include <algorithm>
using std::min;
const int N=1e5+10;
char s[N];
int tax[N],Rank[N],sa[N],sec[N];
int n,m,q,st[N][20],h[N],height[N],Log[N];
void Rsort()
{
    for(int i=1;i<=m;i++) tax[i]=0;
    for(int i=1;i<=n;i++) ++tax[Rank[i]];
    for(int i=2;i<=m;i++) tax[i]+=tax[i-1];
    for(int i=n;i;i--) sa[tax[Rank[sec[i]]]--]=sec[i];
}
bool cmp(int x,int y,int len){return sec[x]==sec[y]&&sec[x+len]==sec[y+len];}
void SuffixSort()
{
    for(int i=1;i<=n;i++) Rank[i]=s[i]-'a'+1,sec[i]=i;
    m=26,Rsort();
    for(int p=0,w=1;p<n;w<<=1,m=p)
    {
        p=0;for(int i=n-w+1;i<=n;i++) sec[++p]=i;
        for(int i=1;i<=n;i++) if(sa[i]>w) sec[++p]=sa[i]-w;
        Rsort(),std::swap(Rank,sec),Rank[sa[p=1]]=1;
        for(int i=2;i<=n;i++) Rank[sa[i]]=cmp(sa[i],sa[i-1],w)?p:++p;
    }
    Log[0]=-1;
    for(int k,i=1;i<=n;height[Rank[i]]=h[i],Log[i]=Log[i>>1]+1,st[i][0]=i,++i)
        for(k=sa[Rank[i]-1],h[i]=h[i-1]?h[i-1]-1:0;s[i+h[i]]==s[k+h[i]];++h[i]);
    for(int j=1;j<=17;j++)
    {
        for(int i=1;i<=n-(1<<j)+1;i++)
        {
            int x=st[i][j-1],y=st[i+(1<<j-1)][j-1];
            st[i][j]=height[x]<height[y]?x:y;
        }
    }
}
int query(int l,int r)
{
    if(l>r) std::swap(l,r);
    if(++l==r) return height[st[l][0]];
    int d=Log[r+1-l],x=st[l][d],y=st[r-(1<<d)+1][d];
    return height[x]<height[y]?height[x]:height[y];
}
int main()
{
    scanf("%s",s+1),n=strlen(s+1);
    SuffixSort();
    scanf("%d",&q);
    for(int x,y,l,r,i=1;i<=q;i++)
    {
        scanf("%d%d%d%d",&x,&y,&l,&r);
        if(x==l) printf("%d\n",min(r,y)+1-l);
        else printf("%d\n",min(query(Rank[x],Rank[l]),min(r+1-l,y+1-x)));
    }
    return 0;
}

---

2019.1.1