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;
}