bzoj 4556 [Tjoi2016&Heoi2016]字符串

4556: [Tjoi2016&Heoi2016]字符串
Time Limit: 20 Sec Memory Limit: 128 MB
Submit: 918 Solved: 365
[Submit][Status][Discuss]
Description

佳媛姐姐过生日的时候，她的小伙伴从某东上买了一个生日礼物。生日礼物放在一个神奇的箱子中。箱子外边写了

一个长为n的字符串s，和m个问题。佳媛姐姐必须正确回答这m个问题，才能打开箱子拿到礼物，升职加薪，出任CE

O，嫁给高富帅，走上人生巅峰。每个问题均有a,b,c,d四个参数，问你子串s[a..b]的所有子串和s[c..d]的最长公

共前缀的长度的最大值是多少？佳媛姐姐并不擅长做这样的问题，所以她向你求助，你该如何帮助她呢？

Input

输入的第一行有两个正整数n,m，分别表示字符串的长度和询问的个数。接下来一行是一个长为n的字符串。接下来

m行，每行有4个数a,b,c,d，表示询问s[a..b]的所有子串和s[c..d]的最长公共前缀的最大值。1<=n,m<=100,000,

字符串中仅有小写英文字母，a<=b,c<=d,1<=a,b,c,d<=n

Output

对于每一次询问，输出答案。

Sample Input

5 5

aaaaa

1 1 1 5

1 5 1 1

2 3 2 3

2 4 2 3

2 3 2 4

Sample Output

1

1

2

2

2

HINT

Source

---

【分析】
半暴力的方法
扫描height数组。
但是应该是不能造出来数据卡住这种方法的…不会证明。

---

【代码】

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#define ll long long#define M(a) memset(a,0,sizeof a)#define fo(i,j,k) for(i=j;i<=k;i++)using namespace std;const int mxn=100005;char s[mxn];int n,m,Q,len;inline int read(){    int x=0;char ch=getchar();    while(ch<'0'||ch>'9') ch=getchar();    while(ch>='0'&&ch<='9') x=(x<<1)+(x<<3)+ch-'0',ch=getchar();    return x;}struct suffix  //大家好我是黑箱 {    int rank[mxn],height[mxn];    int a[mxn],b[mxn],x[mxn],y[mxn],sa[mxn];    inline void get()    {        int i,j;        scanf("%s",s+1);        fo(i,1,len) a[i]=s[i]-'a'+1;        work();    }    inline bool comp(int i,int j,int l)    {        return y[i]==y[j]&&(i+l>len?-1:y[i+l])==(j+l>len?-1:y[j+l]);    }    inline void work()    {        int i,j,k,p;m=30;        fo(i,0,m) b[i]=0;        fo(i,1,len) b[x[i]=a[i]]++;        fo(i,1,m) b[i]+=b[i-1];        for(i=len;i>=1;i--) sa[b[x[i]]--]=i;        for(k=1;k<=len;k<<=1)        {            p=0;            fo(i,len-k+1,len) y[++p]=i;            fo(i,1,len) if(sa[i]>k) y[++p]=sa[i]-k;            fo(i,0,m) b[i]=0;            fo(i,1,len) b[x[y[i]]]++;            fo(i,1,m) b[i]+=b[i-1];            for(i=len;i>=1;i--) sa[b[x[y[i]]]--]=y[i];            swap(x,y),p=2,x[sa[1]]=1;            fo(i,2,len)              x[sa[i]]=comp(sa[i-1],sa[i],k)?p-1:p++;            if(p>len) break;            m=p;        }        p=k=0;        fo(i,1,len) rank[sa[i]]=i;        for(i=1;i<=len;height[rank[i++]]=k)          for(k?k--:0,j=sa[rank[i]-1];a[i+k]==a[j+k];k++);    }}S;int main(){    int i,j,a,b,c,d;    len=read(),Q=read();    S.get();    while(Q--)    {        a=read(),b=read(),c=read(),d=read();        int tmp=S.rank[c],ans=0,mn=len;        if(a<=c && c<=b) ans=min(b-c+1,d-c+1);        fo(i,tmp+1,len)        {            mn=min(mn,S.height[i]);            if(mn<=ans) break;            if(a<=S.sa[i] && S.sa[i]<=b)              ans=max(ans,min(mn,b-S.sa[i]+1));        }        mn=len;        for(i=tmp;i>=2;i--)        {            mn=min(mn,S.height[i]);            if(mn<=ans) break;            if(a<=S.sa[i-1] && S.sa[i-1]<=b)              ans=max(ans,min(mn,b-S.sa[i-1]+1));        }        printf("%d\n",min(ans,d-c+1));    }    return 0;}