HDU 4513 吉哥系列故事――完美队形II （manacher算法）

 

吉哥系列故事――完美队形II

hdu-4513

Description

　　吉哥又想出了一个新的完美队形游戏！ 
　　假设有n个人按顺序站在他的面前，他们的身高分别是h[1], h[2] ... h[n]，吉哥希望从中挑出一些人，让这些人形成一个新的队形，新的队形若满足以下三点要求，则就是新的完美队形： 

　　1、挑出的人保持原队形的相对顺序不变，且必须都是在原队形中连续的； 
　　2、左右对称，假设有m个人形成新的队形，则第1个人和第m个人身高相同，第2个人和第m-1个人身高相同，依此类推，当然如果m是奇数，中间那个人可以任意； 
　　3、从左到中间那个人，身高需保证不下降，如果用H表示新队形的高度，则H[1] <= H[2] <= H[3] .... <= H[mid]。 

　　现在吉哥想知道：最多能选出多少人组成新的完美队形呢？

 

Input

　　输入数据第一行包含一个整数T，表示总共有T组测试数据(T <= 20)； 
　　每组数据首先是一个整数n(1 <= n <= 100000)，表示原先队形的人数，接下来一行输入n个整数，表示原队形从左到右站的人的身高（50 <= h <= 250，不排除特别矮小和高大的）。

 

Output

　　请输出能组成完美队形的最多人数，每组输出占一行。

 

Sample Input

2 3 51 52 51 4 51 52 52 51

 

Sample Output

3 4

 

既然人的身高在50以上，那么就用0和-1来分隔数组搞manacher就好了，在算法的判断条件里加入有关大小的判断，变换一下回文串中的最小值，当然，这个最小值不能是用来分隔的数字。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <set>
using namespace std;

const int maxn = 1e5 + 7;
int n;
int s[maxn];
int ss[maxn<<1];
int p[maxn<<1];

void make() {
    ss[0] = -1;
    int pos = 0;
    for(int i = 0; i < n; i++) {
        ss[++pos] = 0;
        ss[++pos] = s[i];
    }ss[++pos] = 0;
    n = n * 2 + 2;
    //for(int i = 0; i < n; i++) printf("%d ", ss[i]); puts("");
}

void solve() {
    int maxx = 0, id;
    for(int i = 1; i < n; i++) {
        if(maxx > i)
            p[i] = min(p[id * 2 -  i], p[id] + id - i);
        else p[i] = 1;
        int Min;
        if(ss[i] > 0) Min = ss[i];
        else Min = ss[i + 1];
        while(ss[i - p[i]] == ss[i + p[i]] && ss[i - p[i]] <= Min) {
            //printf("are you here?\n");
            if(ss[i - p[i]] > 0)
                Min = ss[i - p[i]];
            if(i + p[i] > maxx) {
                maxx = i + p[i];
                id = i;
            }
            p[i]++;
        }
    }
    int ans = 0;
    for(int i = 1; i < n; i++) ans = max(ans, p[i]);
    printf("%d\n", ans - 1);
}

int main() {
    int t;
    scanf("%d", &t);
    while(t--) {
        scanf("%d", &n);
        for(int i = 0; i < n; i++) {
            scanf("%d", &s[i]);
        }
        make();
        solve();
    }
}