题面:Luogu2852
题目就是让我们求在串s中出现至少k次的最长子串c
首先Height这个性质还是要用上:lcp(x,y)=min(Height[rank[x]]~Height[rank[y]])
不过这题的rank倒没多大用。。。
然后我们可以二分这个长度l(答案啦),然后利用这个性质贪心
这样就可以求出串s中lcp长度大于等于l的最大连续长度,与k比较一下就好了
就是跑得慢了点。。。
具体二分过程看注释好了
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
int n,m,q,Rank[200001],Rank1[200001],sa[200001],h[200001];
int ton[1000001],a[200001];
inline void Sort(){
for(int i=0;i<=m;i++)ton[i]=0;
for(int i=1;i<=n;i++)ton[Rank[Rank1[i]]]++;
for(int i=1;i<=m;i++)ton[i]+=ton[i-1];
for(int i=n;i;i--)sa[ton[Rank[Rank1[i]]]--]=Rank1[i];
}
inline bool check(int x){//求串s中lcp长度大于等于x的最大连续长度
int cnt=0;
for(int i=1;i<=n;i++)if(h[i]//此时lcp已小于x
if(cnt+1>=q)return 1;//已经满足条件
cnt=0;//否则清空再来
}else cnt++;
return cnt+1>=q;
}
int main()
{
scanf("%d%d",&n,&q);
for(int i=1;i<=n;i++)scanf("%d",&a[i]);
for(int i=1;i<=n;i++)Rank[i]=a[i],Rank1[i]=i;
m=1000000;Sort();
for(int i=1,j,p=1;p2,m=p){
for(p=0,j=n-i+1;j<=n;j++)Rank1[++p]=j;
for(j=1;j<=n;j++)if(sa[j]>i)Rank1[++p]=sa[j]-i;
Sort();swap(Rank,Rank1);Rank[sa[1]]=p=1;
for(j=2;j<=n;j++)Rank[sa[j]]=(Rank1[sa[j]]==Rank1[sa[j-1]]&&Rank1[sa[j]+i]==Rank1[sa[j-1]+i])?p:++p;
}
for(int i=1,k=0,j=0;i<=n;h[Rank[i++]]=k)
for(k=k?k-1:k,j=sa[Rank[i]-1];a[i+k]==a[j+k];k++);
int l=0,r=n,ans;//二分长度
while(l<=r){
int mid=l+r>>1;
if(check(mid))ans=mid,l=mid+1;
else r=mid-1;
}
printf("%d",ans);
return 0;
}