hdu 4513 吉哥系列故事——完美队形II(manacher)

Problem 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

 

题意:找到从左到中间不降序列的最长回文长度

思路:manacher在判断的时候多判断是否降序即可

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<stack>
#include<bitset>
#include<cstdlib>
#include<cmath>
#include<set>
#include<list>
#include<deque>
#include<map>
#include<queue>
#define ll long long int
using namespace std;
inline ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
inline ll lcm(ll a,ll b){return a/gcd(a,b)*b;}
int moth[13]={0,31,28,31,30,31,30,31,31,30,31,30,31};
int dir[4][2]={1,0 ,0,1 ,-1,0 ,0,-1};
int dirs[8][2]={1,0 ,0,1 ,-1,0 ,0,-1, -1,-1 ,-1,1 ,1,-1 ,1,1};
const int inf=0x3f3f3f3f;
const ll mod=1e9+7;
int p[200007];
int a[100007];
int temp[200007];
int n;
void manacher(){
    int cnt=2;
    temp[0]=-2;
    temp[1]=-1;
    for(int i=0;i<n;i++){
        temp[cnt++]=a[i];
        temp[cnt++]=-1;
    }
    int po=0; int mx=0;
    for(int i=0;i<cnt;i++){
        p[i]=mx>i?min(mx-i,p[2*po-i]):1;
        if(temp[i]==-1){
            if(temp[i+p[i]]==temp[i-p[i]])
                p[i]++;    
            while(temp[i+p[i]]==temp[i-p[i]]){
                if(temp[i+p[i]]==-1) p[i]++;
                else if(temp[i+p[i]]<=temp[i+p[i]-2]) p[i]++;
                else break;
            }
        }else{
            while(temp[i+p[i]]==temp[i-p[i]]){
                if(temp[i+p[i]]==-1) p[i]++;
                else if(temp[i+p[i]]<=temp[i+p[i]-2]) p[i]++;
                else break;
            }
        }
        if(i+p[i]>mx){
            mx=i+p[i];
            po=i;
        }
    }
}
int main(){
    ios::sync_with_stdio(false);
    int t;
    cin>>t;
    while(t--){
        cin>>n;
        for(int i=0;i<n;i++)
            cin>>a[i];
        manacher();
        int ans=0;
        for(int i=0;i<2*n+2;i++){
        //    cout<<p[i]<<endl;
            ans=max(ans,p[i]-1);
        }
        cout<<ans<<endl;
    }
    return 0;
}