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

这个比最长回文子串就多了一个条件，就是回文字串（这里相当于人的高度）由两端向中间递增。

才刚刚看了看manacher，在用模板A了一道题后，还没有完全理解manacher，然后就准备把这道题也直接带模板的。

却发现这个递增的要求实在是麻烦，然后实在是没办法了，就把manacher算法拆开了，写两遍，按照子串的中间是一个数还是两个数来写。

其他的什么地方都没改，只是扩展的时候把题目要求的条件加上就好了。

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <ctime>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
#define INF 0x3f3f3f3f
#define inf (-((LL)1<<40))
#define lson k<<1, L, mid
#define rson k<<1|1, mid+1, R
#define mem0(a) memset(a,0,sizeof(a))
#define mem1(a) memset(a,-1,sizeof(a))
#define mem(a, b) memset(a, b, sizeof(a))
#define FOPENIN(IN) freopen(IN, "r", stdin)
#define FOPENOUT(OUT) freopen(OUT, "w", stdout)

template<class T> T CMP_MIN(T a, T b) { return a < b; }
template<class T> T CMP_MAX(T a, T b) { return a > b; }
template<class T> T MAX(T a, T b) { return a > b ? a : b; }
template<class T> T MIN(T a, T b) { return a < b ? a : b; }
template<class T> T GCD(T a, T b) { return b ? GCD(b, a%b) : a; }
template<class T> T LCM(T a, T b) { return a / GCD(a,b) * b;    }

//typedef __int64 LL;
typedef long long LL;
const int MAXN = 4000000+10;
const int MAXM = 2000010;
const double eps = 1e-12;

int high[MAXN];
int rad[MAXN];
int T, N;

int manacherOdd()
{
    mem0(rad);
    for(int i=1,j=0,k;i<=N;)
    {
        while (high[i-j-1]==high[i+j+1] && high[i-j]>=high[i-j-1]) j++;  //扫描得出rad值
        rad[i]=j;
        for (k=1; k<=j && rad[i-k]!=rad[i]-k; k++) rad[i+k]=min(rad[i-k],rad[i]-k);  //k指针扫描
        i+=k;  //i跳到下一个需要计算rad值的位置
        j=max(j-k,0);  //更新下一个rad值的初始值
        //printf("%d
", i);
    }
    int ans = 0;
    for(int i=1;i<=N;i++)
    {
        ans = max(ans, rad[i]*2+1);
    }
    return ans;
}

int manacherEven()
{
    mem0(rad);
    int j = 0, k;
    for(int i=1;i<N;)
        if(high[i]==high[i+1])
        {
            while(high[i-j-1] == high[i+j+2] && high[i-j]>=high[i-j-1]) j++;
            rad[i] = j+1;
            for(k=1;k<=j && rad[i-k]!=rad[i]-k; k++) rad[i+k]=min(rad[i-k],rad[i]-k);
            i += k;
            j = max(j-k, 0);
            //printf("%d
", i);
        }
        else i++;
    int ans = 0;
    for(int i=1;i<=N;i++)
    {
        ans = max(ans, rad[i]*2);
    }
    return ans;
}

int main()
{
    while(~scanf("%d", &T))while(T--)
    {
        scanf("%d", &N);
        for(int i=1;i<=N;i++) scanf("%d", &high[i]);
        high[0] = high[N+1] = INF;
        int ans = max( manacherOdd(), manacherEven());
        printf("%d
", ans);
    }
    return 0;
}