Description
佳媛姐姐过生日的时候,她的小伙伴从某东上买了一个生日礼物。生日礼物放在一个神奇的箱子中。箱子外边写了
一个长为n的字符串s,和m个问题。佳媛姐姐必须正确回答这m个问题,才能打开箱子拿到礼物,升职加薪,出任CE
O,嫁给高富帅,走上人生巅峰。每个问题均有a,b,c,d四个参数,问你子串s[a..b]的所有子串和s[c..d]的最长公
共前缀的长度的最大值是多少?佳媛姐姐并不擅长做这样的问题,所以她向你求助,你该如何帮助她呢?
Solution
二分答案
那么只需要判定 (s[c,c+mid-1]) 是否在 ([a+mid-1,b]) 出现过
找到串 (s[c,c+mid-1]) 对应的节点,查一下子树内的 (pos) 集合是否存在 ([a+mid-1,b])
找到对应节点,通过前缀 (s[1,c+mid-1]) 倍增跳就可以了,(pos) 集合用线段树维护一下就行了
#include<bits/stdc++.h>
using namespace std;
template<class T>void gi(T &x){
int f;char c;
for(f=1,c=getchar();c<'0'||c>'9';c=getchar())if(c=='-')f=-1;
for(x=0;c<='9'&&c>='0';c=getchar())x=x*10+(c&15);x*=f;
}
const int N=2e5+10;
int n,Q,ch[N][26],len[N],fa[N],cur=1,cnt=1,pos[N],tt=0,a,b,c,d,f[N][20];
char s[N];int ls[N*30],rs[N*30],head[N],nxt[N],to[N],num=0,rt[N];
inline void link(int x,int y){nxt[++num]=head[x];to[num]=y;head[x]=num;}
inline void ins(int c){
int p=cur;cur=++cnt;len[cur]=len[p]+1;
for(;p && !ch[p][c];p=fa[p])ch[p][c]=cur;
if(!p)fa[cur]=1;
else{
int q=ch[p][c];
if(len[p]+1==len[q])fa[cur]=q;
else{
int nt=++cnt;len[nt]=len[p]+1;
memcpy(ch[nt],ch[q],sizeof(ch[q]));
fa[nt]=fa[q];fa[q]=fa[cur]=nt;
for(;p && ch[p][c]==q;p=fa[p])ch[p][c]=nt;
}
}
}
inline void ins(int &x,int l,int r,int sa){
if(!x)x=++tt;
if(l==r)return ;
int mid=(l+r)>>1;
if(sa<=mid)ins(ls[x],l,mid,sa);
else ins(rs[x],mid+1,r,sa);
}
inline int merge(int x,int y){
if(!x||!y)return x+y;
int o=++tt;
ls[o]=merge(ls[x],ls[y]);
rs[o]=merge(rs[x],rs[y]);
return o;
}
inline void dfs(int x){
for(int i=head[x];i;i=nxt[i])dfs(to[i]),rt[x]=merge(rt[x],rt[to[i]]);
}
inline bool qry(int x,int l,int r,int sa,int se){
if(!x)return 0;
if(sa<=l && r<=se)return 1;
int mid=(l+r)>>1;
if(se<=mid)return qry(ls[x],l,mid,sa,se);
if(sa>mid)return qry(rs[x],mid+1,r,sa,se);
return qry(ls[x],l,mid,sa,mid)|qry(rs[x],mid+1,r,mid+1,se);
}
inline bool check(int mid,int r){
int x=pos[r];
for(int i=19;i>=0;i--)
if(f[x][i] && len[f[x][i]]>=mid)x=f[x][i];
return qry(rt[x],1,n,a+mid-1,b);
}
int main(){
freopen("pp.in","r",stdin);
freopen("pp.out","w",stdout);
cin>>n>>Q;
scanf("%s",s+1);
for(int i=1;i<=n;i++)ins(s[i]-'a'),ins(rt[cur],1,n,i),pos[i]=cur;
for(int i=2;i<=cnt;i++)link(fa[i],i);
for(int i=1;i<=cnt;i++)f[i][0]=fa[i];
for(int j=1;j<20;j++)
for(int i=1;i<=cnt;i++)f[i][j]=f[f[i][j-1]][j-1];
dfs(1);
while(Q--){
gi(a);gi(b);gi(c);gi(d);
int l=1,r=min(b-a+1,d-c+1),mid,ret=0;
while(l<=r){
mid=(l+r)>>1;
if(check(mid,c+mid-1))ret=mid,l=mid+1;
else r=mid-1;
}
printf("%d
",ret);
}
return 0;
}