1717: [Usaco2006 Dec]Milk Patterns 产奶的模式
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
1
2
3
2
3
2
3
1
Sample Output
4
HINT
Source
二分答案+height数组判定最长公共前缀。。。
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<iostream> 2 #include<cstdlib> 3 #include<cmath> 4 #include<cstring> 5 #include<cstdio> 6 #include<algorithm> 7 #include<string> 8 #include<map> 9 #include<queue> 10 #include<vector> 11 #include<set> 12 #define inf 1000000000 13 #define maxn 20000+1000 14 #define maxm 1000000+100 15 #define eps 1e-10 16 #define ll long long 17 #define for0(i,n) for(int i=0;i<=(n);i++) 18 #define for1(i,n) for(int i=1;i<=(n);i++) 19 #define for2(i,x,y) for(int i=(x);i<=(y);i++) 20 #define for3(i,x,y) for(int i=(x);i>=(y);i--) 21 #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go) 22 using namespace std; 23 string s,ss; 24 int n,k,a[maxn],sa[maxn],t[maxn],t2[maxn],c[maxm],h[maxn],rk[maxn]; 25 int read(){ 26 int x=0,f=1;char ch=getchar(); 27 while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} 28 while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();} 29 return x*f; 30 } 31 void getsa(int n,int m) 32 { 33 int *x=t,*y=t2; 34 for0(i,m)c[i]=0; 35 for0(i,n)c[x[i]=a[i]]++; 36 for1(i,m)c[i]+=c[i-1]; 37 for3(i,n,0)sa[--c[x[i]]]=i; 38 for(int k=1;k<=n+1;k<<=1) 39 { 40 int p=0; 41 for2(i,n-k+1,n)y[p++]=i; 42 for0(i,n)if(sa[i]>=k)y[p++]=sa[i]-k; 43 for0(i,m)c[i]=0; 44 for0(i,n)c[x[y[i]]]++; 45 for1(i,m)c[i]+=c[i-1]; 46 for3(i,n,0)sa[--c[x[y[i]]]]=y[i]; 47 swap(x,y);p=1;x[sa[0]]=0; 48 for1(i,n)x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?p:++p; 49 if(p>=n)break; 50 m=p; 51 } 52 for1(i,n)rk[sa[i]]=i; 53 for(int i=0,k=0,j;i<n;h[rk[i++]]=k) 54 for(k?k--:0,j=sa[rk[i]-1];a[i+k]==a[j+k];k++); 55 } 56 int main(){ 57 //freopen("input.txt","r",stdin); 58 //freopen("output.txt","w",stdout); 59 n=read();k=read(); 60 for0(i,n-1)a[i]=read()+1; 61 a[n]=0; 62 getsa(n,1000010); 63 int l=1,r=n; 64 while(l<=r){ 65 int mid=(l+r)>>1,tot=1; 66 for1(i,n){ 67 if(h[i]>=mid){ 68 tot++; 69 if(tot>=k){l=mid+1;break;} 70 } 71 else tot=1; 72 } 73 if(tot<k)r=mid-1; 74 } 75 printf("%d",r); 76 return 0; 77 }