bzoj 4556 字符串

后缀数组,暴力硬跑
贼快

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define MAXN 100005
using namespace std;
int buc[MAXN],wa[MAXN],wb[MAXN];
int r[MAXN],sa[MAXN],rank[MAXN],height[MAXN];
void getheight(int n){
    int i,j,k=0;
    for(i=0;i<n;i++)rank[sa[i]]=i;
    for(i=0;i<n;height[rank[i++]]=k)
        for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
    return ;
}
bool cmp(int *c,int a,int b,int d)
    {return c[a]==c[b]&&c[a+d]==c[b+d];}
void da(int n,int m=130){
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0;i<m;i++)buc[i]=0;
    for(i=0;i<n;i++)buc[x[i]=r[i]]++;
    for(i=1;i<m;i++)buc[i]+=buc[i-1];
    for(i=n-1;~i;i--)sa[--buc[x[i]]]=i;
    for(j=1,p=1;p<n;j*=2,m=p){
        for(i=n-j,p=0;i<n;i++)y[p++]=i;
        for(i=0;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;
        for(i=0;i<m;i++)buc[i]=0;
        for(i=0;i<n;i++)buc[x[y[i]]]++;
        for(i=1;i<m;i++)buc[i]+=buc[i-1];
        for(i=n-1;~i;i--)sa[--buc[x[y[i]]]]=y[i];
        for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
    }
    getheight(n);
    return ;
}
int n,q;
char s[MAXN];
int main(){
    scanf("%d%d",&n,&q);
    scanf("%s",s);
    for(int i=0;i<n;i++)r[i]=s[i];
    da(n+1);
    int ans,minn,pos,now,a,b,c,d;
    while(q--){
        scanf("%d%d%d%d",&a,&b,&c,&d);
        a--;b--;c--;d--;
        ans=0; minn=min(d-c+1,b-a+1); 
        pos=rank[c]; 
        if(sa[pos]>=a&&sa[pos]<=b)ans=max(ans,min(minn,min(b-sa[pos]+1,d-c+1)));
        minn=min(minn,height[pos]);
        for(int i=pos-1;i>=0;i--){
            if(minn<=ans)break;
            now=sa[i];
            if(now>=a&&now<=b)ans=max(ans,min(minn,min(b-sa[i]+1,d-c+1)));
            minn=min(minn,height[i]);
        }
        minn=min(d-c+1,b-a+1);
        for(int i=pos+1;i<=n;i++){
            minn=min(minn,height[i]);
            if(minn<=ans)break;
            now=sa[i];
            if(now>=a&&now<=b){ans=max(ans,min(minn,min(b-sa[i]+1,d-c+1)));}
        }
        printf("%d\n",ans);
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值