bzoj1717[Usaco2006 Dec]Milk Patterns 产奶的模式*

bzoj1717[Usaco2006 Dec]Milk Patterns 产奶的模式

题意：

John记录了n天的牛奶质量值。他想知道最长的出现了至少k次的模式（即一个连续子串）的长度。n≤20000。

题解：

求一个总哈希值，二分模式长度，枚举每个模式开始端点，获得它的哈希值，然后排序比较有没有至少k个相等。

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#define inc(i,j,k) for(int i=j;i<=k;i++)
#define maxn 20010
#define ll unsigned long long
#define yyl 2333
using namespace std;

inline int read(){
    char ch=getchar(); int f=1,x=0;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return f*x;
}
ll hash[maxn],mi[maxn],calc[maxn],pat[maxn]; int n,k,l,r,ans,tot;
int main(){
    n=read(); k=read(); inc(i,1,n)pat[i]=read(); mi[0]=1;
    inc(i,1,n)hash[i]=hash[i-1]*yyl+pat[i],mi[i]=mi[i-1]*yyl; l=1; r=n;
    while(l<=r){
        int mid=(l+r)>>1; tot=0; inc(i,mid,n)calc[++tot]=hash[i]-hash[i-mid]*mi[mid];
        sort(calc+1,calc+tot+1); int x=0,y=0;
        inc(i,1,tot){if(!x||calc[i]!=calc[x])x=i; y=max(y,i-x+1);}
        if(y>=k)ans=mid,l=mid+1;else r=mid-1;
    }
    printf("%d",ans); return 0;
}

20161102