Stammering Aliens

本文介绍了一种使用后缀数组和最长公共前缀(LCP)技术来查找字符串中出现次数大于等于k的最长子串的方法。通过二分查找确定满足条件的最长子串长度,并记录最右侧起始位置。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题意:
输入一个k,和一个字符串,询问是否存在出现次数>=k的子串,若存在,输出长度最长的子串和这些子串中起始下标最靠右的(下标从0开始),否则输出"none";

例如:
3
baaaababababbababbab
其中babab出现了3次,长度为5,起始下标分别出现在5,7,12,输出最靠右的也就是12


题解:
先用SA对所有后缀排序,由于一段区间的LCP表示的子串一定出现在每个后缀里,因此出现次数就可以解决了,也就是一段连续区间内LCP>=x的个数,然后考虑怎么得到最长的子串,我们可以二分,二分长度mid,若一段区间内LCP>=mid的个数>=k,那么就存在这个长度的子串出现次数>=k,至于靠右的位置,check的时候记录max(sa)即可


AC代码:

#pragma GCC optimize(2)
#include<bits/stdc++.h>
#include<ext/rope>
using namespace std;
using namespace __gnu_cxx;
#define LL long long
const int MAXN = 40000+50;
const int MAXM = 2e5+50;
const int MOD = 1e9+7;
const int INF = 0x3f3f3f3f;
const double PI = acos(-1);
int n,k;
char s[MAXN];
int rk[MAXN],sa[MAXN],h[MAXN],c[MAXN],y[MAXN];
inline void SA(int n,int m){
    for(int i=0;i<m;i++) c[i]=0;
    for(int i=0;i<n;i++) c[rk[i]=s[i]]++;
    for(int i=1;i<m;i++) c[i]+=c[i-1];
    for(int i=n-1;i>=0;i--) sa[--c[rk[i]]]=i;
    for(int k=1;k<=n;k<<=1){
        int num=0;
        for(int i=n-k;i<n;i++) y[num++]=i;
        for(int i=0;i<n;i++) if(sa[i]>=k) y[num++]=sa[i]-k;
        for(int i=0;i<m;i++) c[i]=0;
        for(int i=0;i<n;i++) c[rk[i]]++;
        for(int i=1;i<m;i++) c[i]+=c[i-1];
        for(int i=n-1;i>=0;i--) sa[--c[rk[y[i]]]]=y[i];
        swap(rk,y);
        rk[sa[0]]=0,num=1;
        for(int i=1;i<n;i++)
            rk[sa[i]]=(y[sa[i]]==y[sa[i-1]] && y[sa[i]+k]==y[sa[i-1]+k] ? num-1 : num++);
        if(num>=n) break;
        m=num;
    }
}
inline void H(int n){
    for(int i=1;i<=n;i++) rk[sa[i]]=i;
    for(int i=0,k=0;i<n;i++){
        if(k) k--; else k=0;
        int j=sa[rk[i]-1];
        while(s[i+k]==s[j+k]) k++;
        h[rk[i]]=k;
    }
}
inline int check(int mid){
    int mx=-1; int cnt=0,ans=-1;
    for(int i=1;i<=n;i++){
        if(h[i]>=mid){
            cnt++; ans=max(ans,sa[i]); ans=max(ans,sa[i-1]);
            if(cnt>=k) mx=max(mx,ans);
        }else{
            cnt=1; ans=-1;
        }
    }
    return mx;
}
signed main(){
#ifndef ONLINE_JUDGE
    freopen("C:\\Users\\Administrator\\Desktop\\in.txt","r",stdin);
#endif // ONLINE_JUDGE
    while(~scanf("%d",&k) && k){
        scanf("%s",s); n=strlen(s); s[n]=0;
        if(k==1){ printf("%d 0\n",n);continue; }
        SA(n+1,150); H(n);
        int l=1,r=n,res=-1,mx=-1;
        while(l<=r){
            int mid=(l+r)>>1;
            int ans=check(mid);
            if(ans!=-1) l=mid+1,res=mid,mx=ans;
            else r=mid-1;
        }
        if(res==-1) puts("none"); else printf("%d %d\n",res,mx);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值