bzoj 1717 [USACO06DEC] Milk Patterns (后缀数组+二分)

题目大意:求可重叠的相同子串数量至少是K的子串最长长度  

洛谷传送门

依然是后缀数组+二分,先用后缀数组处理出height

每次二分出一个长度x,然后去验证,在排序的后缀串集合里,有没有连续数量多于K个串的长度>=x,

但据说有一种高端做法是把二分换成单调队列,能减少常数,可惜我并没有看懂......

原题好像是哈希的骚操作,但网上的题解好像都是后缀数组......

比上一道男人八题简单多了,我原来的错代码竟然卡过去了70分......

#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 20010
#define inf 0x3f3f3f3f
#define maxn 1000010
using namespace std;

int n,len,tot,K;
int a[N],tr[N],sa[N],rk[N],hs[N],h[N];
struct node{int id,val,w;}d[N];
int cmp1(node s1,node s2){return s1.val<s2.val;}
int cmp2(node s1,node s2){return s1.id<s2.id;}
int gint()
{
    int rett=0,fh=1;char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')fh=-1;c=getchar();}
    while(c>='0'&&c<='9'){rett=(rett<<3)+(rett<<1)+c-'0';c=getchar();}
    return rett*fh;
}
bool check(int k,int x,int y){
    if(x+k>len||y+k>len) return 0;
    else return (rk[x]==rk[y]&&rk[x+k]==rk[y+k])?1:0;
}
void get_sa()
{
    int cnt=0,i;
    for(i=1;i<=len;i++) hs[a[i]]++;
    for(i=1;i<=tot;i++) if(hs[i]) tr[i]=++cnt;
    for(i=1;i<=tot;i++) hs[i]+=hs[i-1];
    for(i=1;i<=len;i++) rk[i]=tr[a[i]],sa[hs[a[i]]--]=i;
    for(int k=1;cnt<len;k<<=1)
    {
        for(i=1;i<=cnt;i++) hs[i]=0;
        for(i=1;i<=len;i++) hs[rk[i]]++;
        for(i=1;i<=cnt;i++) hs[i]+=hs[i-1];
        for(i=len;i>=1;i--) if(sa[i]>k) tr[sa[i]-k]=hs[rk[sa[i]-k]]--;
        for(i=1;i<=k;i++) tr[len-i+1]=hs[rk[len-i+1]]--;
        for(i=1;i<=len;i++) sa[tr[i]]=i;
        for(i=1,cnt=0;i<=len;i++) tr[sa[i]]=check(k,sa[i],sa[i-1])?cnt:++cnt;
        for(i=1;i<=len;i++) rk[i]=tr[i];
    }
}
void get_height()
{ 
    for(int i=1;i<=n;i++){
        if(rk[i]==1) continue;
        for(int j=max(1,h[rk[i-1]]-1);;j++)
            if(a[i+j-1]==a[sa[rk[i]-1]+j-1]) h[rk[i]]=j;
            else break;
    }
}
void descrete()
{
    sort(d+1,d+n+1,cmp1);
    d[0].val=-1;
    for(int i=1;i<=n;i++)
        if(d[i].val==d[i-1].val)
            d[i].w=d[i-1].w;
        else d[i].w=++tot;
    sort(d+1,d+n+1,cmp2);
    for(int i=1;i<=n;i++)
        a[i]=d[i].w;
}
int check(int ans)
{
    int i,cnt=0;
    for(i=1;i<=n;)
    {
        if(h[i]<ans) {i++;continue;}
        for(cnt=0;h[i]>=ans&&i<=n;i++) 
            cnt++;
        if(cnt>=K-1) return 1;
    }return 0;
}

int main()
{
    scanf("%d%d",&n,&K);len=n;
    for(int i=1;i<=n;i++) 
        d[i].val=gint(),d[i].id=i;
    descrete();
    get_sa();
    get_height();
    int l=0,r=n,ans=0;
    while(l<=r){
        int mid=(l+r)>>1;
        if(check(mid)) ans=mid,l=mid+1;
        else r=mid-1;
    }printf("%d\n",ans);
    return 0;
}


 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值