【BZOJ1717】产奶的模式 后缀数组

原创 2016年06月01日 20:33:15

无脑后缀数组大法好QAQ

还是在height数组上做文章,等价于找序列中【连续K-1个height值中的最小值】的最大值,滑动窗口维护更新即可。

第一次发现pair如此好♂用

/**************************************************************
    Problem: 1717
    User: cqyzhb
    Language: C++
    Result: Accepted
    Time:20 ms
    Memory:2564 kb
****************************************************************/
 
#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
#define MAXN 40005
void _read(int &x)
{
    x=0; char ch=getchar(); bool flag=false;
    while(ch<'0' || ch>'9'){if(ch=='-')flag=true; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();} if(flag)x=-x; return ;
}
int N,K,A[MAXN],M,B[MAXN];
char s[MAXN];
void Init()
{
    _read(N); _read(K); K--;
    for(int i=1;i<=N;_read(A[i]),B[i]=A[i],i++);
    sort(B+1,B+1+N);
    M=unique(B+1,B+1+N)-B-1;
    for(int i=1;i<=N;i++)A[i]=lower_bound(B+1,B+1+M,A[i])-B;
    return ;
}
int wa[MAXN],wb[MAXN],wc[MAXN],sa[MAXN],h[MAXN],rank[MAXN];
void make_sa()
{
    int *t,*x=wa,*y=wb,i,k,cnt;
    for(i=0;i<=M;i++)wc[i]=0;
    for(i=1;i<=N;i++)wc[x[i]=A[i]]++;
    for(i=1;i<=M;i++)wc[i]+=wc[i-1];
    for(i=N;i>=1;i--)sa[wc[x[i]]--]=i;
    for(k=1;k<=N;k=k*2)
    {
        cnt=0;
        for(i=N-k+1;i<=N;i++)y[++cnt]=i;
        for(i=1;i<=N;i++)if(sa[i]>k)y[++cnt]=sa[i]-k;
        for(i=0;i<=M;i++)wc[i]=0;
        for(i=1;i<=N;i++)wc[x[y[i]]]++;
        for(i=1;i<=M;i++)wc[i]+=wc[i-1];
        for(i=N;i>=1;i--)sa[wc[x[y[i]]]--]=y[i];
        for(t=x,x=y,y=t,cnt=1,x[sa[1]]=1,i=2;i<=N;i++)
        x[sa[i]]= (y[sa[i-1]]==y[sa[i]])&&(y[sa[i-1]+k]==y[sa[i]+k]) ? cnt : ++cnt ;
        M=cnt; if(cnt==N)break;
    }
    for(i=1;i<=N;i++)rank[sa[i]]=i;
    k=0;
    for(int i=1;i<=N;h[rank[i]]=k,i++)
    {
        if(k)k--;for(int j=sa[rank[i]-1];A[i+k]==A[j+k];k++);
    }
    return ;
}
priority_queue<pair<int,int>,vector<pair<int,int> > >pq;
int main()
{
//  freopen("in.txt","r",stdin);
    Init();
    make_sa();
    int ans=0;
    for(int i=1;i<K;i++)pq.push(make_pair(-h[i],i));
    for(int i=K;i<=N;i++)
    {
        pq.push(make_pair(-h[i],i));
        while(pq.top().second <= i-K)pq.pop();
        ans=max(ans,-pq.top().first );
    }
    printf("%d\n",ans);
    return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

[BZOJ1717][Usaco2006 Dec]Milk Patterns 产奶的模式(后缀数组)

错过的年华在北漠开出斑斓的紫薇花,却荒芜了轮回的春夏。

BZOJ1717: [Usaco2006 Dec]Milk Patterns 产奶的模式 后缀数组

1717: [Usaco2006 Dec]Milk Patterns 产奶的模式 Time Limit: 5 Sec  Memory Limit: 64 MB Submit: 1064  Solv...
  • Oakley_
  • Oakley_
  • 2017年01月12日 17:05
  • 316

bzoj1717[Usaco2006 Dec]Milk Patterns 产奶的模式 (二分+hash)题解①

1717: [Usaco2006 Dec]Milk Patterns 产奶的模式 Time Limit: 5 Sec  Memory Limit: 64 MB Submit: 993  Solved:...

BZOJ1717——[Usaco2006 Dec]Milk Patterns 产奶的模式

1、题意:给个数列板字符串,求出现了至少K(2≤K≤N)K(2\leq K\leq N)次的字符串的长度 2、分析:字符串题,考虑heightheight数组的性质,利用它我们可以求出一些字符串之间...

后缀数组(bzoj 1031: [JSOI2007]字符加密Cipher)

后缀数组主要功能: 长度为n的字符串总共有n个后缀,求这n个后缀的字典序 实现方法:倍增+基数排序,过程就是下面那张表 求log(n)次rank数组,每次的rank数组都可以通过上次的rank数组...

[BZOJ2534]Uva10829L-gap字符串(后缀数组+st表)

题目描述传送门题目大意:求字符串s中有多少子串,满足形如ABA形式,其中A是非空字符串,且B的长度正好为L题解这道题和股市的预测实际上时一样的…不过现在忘得快干净了… B的长度已知是L,首先枚举A的...

【BZOJ4310】跳蚤,后缀数组+ST表求LCP+二分答案

bzoj 2119: 股市的预测 (后缀数组+差分+st表)

2119: 股市的预测 Time Limit: 10 Sec  Memory Limit: 259 MB Submit: 283  Solved: 135 [Submit][Status][Discu...

BZOJ 2251: [2010Beijing Wc]外星联络|后缀数组

后缀数组求出height数组 然后扫一边 暴力统计出解就可以了 #include #include #include #include #include #include #include #i...
  • ws_yzy
  • ws_yzy
  • 2016年01月12日 15:30
  • 434

【BZOJ4698】【Sdoi2008】Sandy的卡片(kmp 后缀数组)

DescriptionSandy和Sue的热衷于收集干脆面中的卡片。然而,Sue收集卡片是因为卡片上漂亮的人物形象,而Sandy则是为了积攒卡片兑换超炫的人物模型。每一张卡片都由一些数字进行标记,第i...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【BZOJ1717】产奶的模式 后缀数组
举报原因:
原因补充:

(最多只允许输入30个字)