后缀数组——Luogu2852 [USACO06DEC]牛奶模式Milk Patterns

题面: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;
}

关于这题我们的讲课人的课件给出的题解完全看不懂。。。
能够解答的同学请发下面讨论,十分感谢!
后缀数组——Luogu2852 [USACO06DEC]牛奶模式Milk Patterns_第1张图片

你可能感兴趣的:(贪心,二分/三分,后缀数组/后缀自动机)