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

中文题面，大概算得上模板题...不过在模板的while部分加上限制条件即可..因为是整数数组，都是正数，所以manacher模板里的#，？,\0什么的可以修改为-1 -2 -3...题目让判断回文从左到中间非递减，那么如果是以i为中心，那么b[i-1]和b[i+1]如果回文就是相等的，并且假定b[i-x]非负，那么b[i-x+2]>=b[i-x]对于所有x都符合情况。所以加上这个限制条件即可...

#include"cstdio"
#include"iostream"
#include"cstring"
#include"algorithm"
#define N 110005
using namespace std;
int a[N],b[N*2];
int rad[N*2];  // rad[i]代表以i为中心的最大回文半径
int main()
{
    int t;
    cin>>t;
    for(int cas=1;cas<=t;cas++)
    {       //以下标1读入方便使用
        int maxl,maxid,id;
        int i,len;
        scanf("%d",&len);
        for(int i=1;i<=len;i++)
            scanf("%d",&a[i]);
        for(i=1;i<=len;i++)
        {         // 插入'#' 使得所有回文串为奇数
            b[i*2]=a[i];
            b[i*2+1]=-1;
        }
        len=2*i;b[0]=-2; b[1]=-1; b[len]=-3;   // 防越界  头尾区别
        maxid=id=0;
        maxl=0;
        for(i=1;ii)  rad[i]=min(rad[2*id-i],maxid-i);  // 利用对称性优化复杂度
            else  rad[i]=1;
            while(b[i-rad[i]]==b[i+rad[i]])
            {       //对于当前点 搜索半径
                if(b[i-rad[i]]>0&&b[i-rad[i]]>b[i-rad[i]+2])
                    break;
                rad[i]++;
            }
            if(rad[i]+i>maxid)
            {
                maxid=rad[i]+i;
                id=i;
            }
            if(rad[i]>maxl) maxl=rad[i];
        }
        printf("%d\n",maxl-1);
    }
    return 0;
}