【题意】
给定一个长度为n的整数序列,求其中至少出现k次的子序列长度最长为多长
【输入】
第一行n和k
接下来n个数字描述序列
【输出】
一个数字,表示至少出现过k次的子序列最长长度
也是09年论文《后缀数组——处理字符串的有力工具》上的例题
不过不同之处是二分答案之时,判断方式是一组中是否有k个后缀
/*
* POJ_3261.cpp
*
* Created on: 2013年11月25日
* Author: Administrator
*/
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn = 20010;
struct number{
int x;
int pos;
}num[maxn];
struct node{
int now;
int next;
}d[maxn];
int val[maxn][2];
int c[maxn];
int rank[maxn];
int sa[maxn];
int pos[maxn];
int x[maxn];
int n;
int k;
int h[maxn];
int height[maxn];
bool cmp(const number& a,const number& b){
return a.x < b.x;
}
void add_value(int u, int v, int i) {
d[i].next = c[u];
c[u] = i;
d[i].now = v;
}
void radix_sort(int l, int r) {
for(int k = 1 ; k >= 0 ; --k){
memset(c, 0, sizeof(c));
for(int i = r ; i >= l ; --i){
add_value(val[pos[i]][k], pos[i], i);
}
int t = 0;
for(int i = 0 ; i <= 20000 ; ++i){
for (int j = c[i]; j; j = d[j].next){
pos[++t] = d[j].now;
}
}
}
int t = 0;
for(int i = 1 ; i <= n ; ++i){
if (val[pos[i]][0] != val[pos[i - 1]][0]|| val[pos[i]][1] != val[pos[i - 1]][1]){
t++;
}
rank[pos[i]] = t;
}
}
bool exist(int len) {
int now = 0;
int s = 0;
for(int i = 1 ; i <= n ; ++i){//枚举名次数组...
if (height[i] < len) {
s = max(s,now);//结束当前组
now = 1;//now恢复为1
} else {
now++;
}
}
s = max(s,now);
if(s >= k){
return 1;
}
return 0;
}
void get_suffix_array() {
int t = 1;
while (t / 2 <= n) {
for(int i = 1 ; i <= n ; ++i){
val[i][0] = rank[i];
val[i][1] = (((i + t / 2 <= n) ? rank[i + t / 2] : 0));
pos[i] = i;
}
radix_sort(1, n);
t *= 2;
}
for(int i = 1 ; i <= n ; ++i){
sa[rank[i]] = i;
}
}
void get_common_prefix() {
memset(h, 0, sizeof(h));
for(int i = 1 ; i <= n ; ++i){
if (rank[i] == 1){
h[i] = 0;
}
else {
int now = 0;
if (i > 1 && h[i - 1] > 1){
now = h[i - 1] - 1;
}
while (now + i <= n && now + sa[rank[i] - 1] <= n&& x[now + i] == x[now + sa[rank[i] - 1]]){
now++;
}
h[i] = now;
}
}
for(int i = 1 ; i <= n ; ++i){
height[rank[i]] = h[i];
}
}
int binary_search(int l, int r) {
while (l <= r) {
int mid = (l + r) / 2;
if (exist(mid)){
l = mid + 1;
}else{
r = mid - 1;
}
}
return r;
}
void discretization(){//分组
sort(num+1,num+1+n,cmp);
int t = 0;
num[0].x = -1;
for(int i = 1 ; i <= n ; ++i){
if(num[i].x != num[i-1].x){
t++;
}
rank[num[i].pos] = t;
}
}
void solve() {
discretization();
get_suffix_array();
get_common_prefix();
int ans = binary_search(0,n);
printf("%d\n",ans);
}
int main(){
while(scanf("%d%d",&n,&k)!=EOF){
int i;
for(i = 1 ; i <= n ; ++i){
scanf("%d",&x[i]);
}
for(i = 1 ; i <= n ; ++i){
num[i].x = x[i];
num[i].pos = i;
}
solve();
}
return 0;
}