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 = +;
 const int MAXM = ;
 const double eps = 1e-;

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

 int manacherOdd()
 {
     mem0(rad);
     for(int i=,j=,k;i<=N;)
     {
         while (high[i-j-]==high[i+j+] && high[i-j]>=high[i-j-]) j++;  //扫描得出rad值
         rad[i]=j;
         for (k=; 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,);  //更新下一个rad值的初始值
         //printf("%d\n", i);
     }
     int ans = ;
     for(int i=;i<=N;i++)
     {
         ans = max(ans, rad[i]*+);
     }
     return ans;
 }

 int manacherEven()
 {
     mem0(rad);
     int j = , k;
     for(int i=;i<N;)
         if(high[i]==high[i+])
         {
             while(high[i-j-] == high[i+j+] && high[i-j]>=high[i-j-]) j++;
             rad[i] = j+;
             for(k=;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, );
             //printf("%d\n", i);
         }
         else i++;
     int ans = ;
     for(int i=;i<=N;i++)
     {
         ans = max(ans, rad[i]*);
     }
     return ans;
 }

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