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

1717: [Usaco2006 Dec]Milk Patterns 产奶的模式

Time Limit: 5 Sec   Memory Limit: 64 MB
Submit: 993   Solved: 549
[ Submit][ Status][ Discuss]

Description

农夫John发现他的奶牛产奶的质量一直在变动。经过细致的调查，他发现：虽然他不能预见明天产奶的质量，但连续的若干天的质量有很多重叠。我们称之为一个“模式”。 John的牛奶按质量可以被赋予一个0到1000000之间的数。并且John记录了N(1<=N<=20000)天的牛奶质量值。他想知道最长的出现了至少K(2<=K<=N)次的模式的长度。比如1 2 3 2 3 2 3 1 中 2 3 2 3出现了两次。当K=2时，这个长度为4。

Input

* Line 1: 两个整数 N,K。

* Lines 2..N+1: 每行一个整数表示当天的质量值。

Output

* Line 1: 一个整数：N天中最长的出现了至少K次的模式的长度

Sample Input

8 2
1
2
3
2
3
2
3
1

Sample Output

4

HINT

Source

Gold

昨晚学后缀数组的时候前辈提到这道题，对后缀数组这种烧脑的玩意极其厌恶的我瞬间yy了hash的解法，今天交一发单hash及unsigned自然溢出试试，这么水的hash写法都过了真是开心。

等学会后缀数组再回来用后缀数组A掉这题

hash题解：

二分长度len，判断时O(n)枚举长度为len的子串计算hash值存入map，判断是否有出现≥k次的子串。

代码：

/**************************************************************
    Problem: 1717
    User: **********
    Language: C++
    Result: Accepted
    Time:180 ms
    Memory:4560 kb
****************************************************************/
 
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<climits>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define N 100010
using namespace std;
typedef unsigned long long ll;
map<ll,int>M;
const ll sed=19963;
int n,k,a[N],ans;
ll hash[N],pw[N];
int num[N],sav[N],tot;
ll geth(int x,int y){return hash[y]-hash[x-1]*pw[y-x+1];}
int find(int x)
{
    int l=1,r=tot,ret;
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(sav[mid]<=x)
            l=mid+1,ret=mid;
        else r=mid-1;
    }return ret;
}
bool check(int len)
{
    M.clear();
    for(int i=1;i+len-1<=n;i++)
    {
        int x=i,y=i+len-1;
        ll tmp=geth(x,y);M[tmp]++;
        if(M[tmp]==k)return true;
    }return false;
}
void solve()
{
    int l=1,r=n;
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(check(mid))
            l=mid+1,ans=mid;
        else r=mid-1;
    }
}
int main()
{
    scanf("%d%d",&n,&k);pw[0]=1;
    for(int i=1;i<=n;i++)pw[i]=pw[i-1]*sed;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        num[i]=a[i];
    }sort(num+1,num+1+n);
    for(int i=1;i<=n;i++)
        if(num[i]!=num[i-1])
            sav[++tot]=num[i];
    for(int i=1;i<=n;i++)
    {
        int t=find(a[i]);
        hash[i]=hash[i-1]*sed+t;
    }solve();printf("%d\n",ans);
}